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Executive  Summary 


The  main  thrust  of  this  effort  was  centered  around  the  actual  dy- 
namic, rather  than  static,  thermal  behavior  and  response  of  the  fabric 
of  an  experimental  masonry  building.  Presently,  most  thermal  design 
procedures  for  a building  involve  the  assumptions  that  steady-state 
temperatures  exist  indoors  and  outdoors  and  that  the  mass  of  the  build- 
ing can  be  neglected.  This  report  describes  the  dynamic  (non-steady- 
state)  thermal  behavior  of  an  experimental  building  as  it  is  affected 
by  changing  outdoor  air  temperatures. 

A full-scale  building  was  erected  in  a high-bay  environmental 
chamber  and  the  exterior  surfaces  were  exposed  to  a diurnal  temperature 
cycle.  Several  features  of  the  building  were  changed  during  the  experi 
ment  to  note  the  effect  on  the  thermal  performance  of  the  building. 
These  features  were:  fenestration,  amount  and  location  of  insulation 

and  indoor  mass.  In  all  cases  measured  values  of  temperature  and  heat 
transfer  were  compared  with  the  corresponding  values  predicted  by  a 
computer  program  called  National  Bureau  of  Standards  Load  Determination 
(NBSLD). 
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The  experimental  structure  was  a one  room  house  which  was  20*  long, 
20*  wide,  and  10*  high.  The  walls  were  made  of  solid  concrete  cinder 
aggregate  blocks  with  fully  bedded  mortar  joints.  The  floor  consisted 
of  two  inch  thick  concrete  placed  over  two  inches  of  polystyrene 
board-type  insulation.  The  roof  was  made  from  five  reinforced  pre-cast 
concrete  slabs.  When  it  was  desired  to  simulate  an  indoor  mass,  2600 
pounds  of  concrete  block  were  stacked  on  the  floor.  When  insulation  was 
used,  2 inch  thick  polystyrene  board-type  insulation  was  spot-glued  to 
the  inside  or  the  outside  surfaces  of  the  building.  This  building  was 
instrumented  to  obtain  heat  transfer  data. 

Two  types  of  basic  tests  were  performed.  The  first  was  a floating 
test  in  which  no  heat  energy  was  added  to  or  taken  away  from  the  interior 
of  the  building.  For  this  test,  the  interior  thermal  environment  of  the 
building  was  allowed  to  respond  to  the  outdoor  air  temperature  cycle. 

For  most  of  the  tests  the  outdoor  temperature  was  varied  between  40  and 
100  F each  day.  The  second  type  of  test  was  a thermostated  test.  For 
these  tests  the  experimental  building  was  exposed  to  the  outdoor  air 
temperature  cycle,  while  the  indoor  air  temperature  was  maintained  within 
+ 1 °F  by  controlling  four  electric  fan  heaters  located  on  the  floor 
of  the  experimental  structure. 


iii 


It  was  found  that  the  combination  of  mass  in  the  walls  and  roof 
facing  the  interior  with  insulation  placed  on  the  outside  surfaces  of 
the  building  was  very  effective  in  reducing  and  controlling  the  variation 
of  the  indoor  air  temperature o This  desired  effect  was  predicted  by  the 
computer  program.  For  example,  when  the  inside  air  temperature  was  not 
controlled  and  the  building  was  floating  in  response  to  the  outside  air 
temperature  cycle  (about  60  °F  change)  the  indoor  air  temperature  change 
over  24  hours  was  about  + 1 °F.  In  addition,  comparing  cases  of  no  in- 
sulation, insulation  inside,  and  insulation  outside,  the  temperature 
differences  from  floor  to  ceiling  on  the  walls  and  of  the  indoor  air 
were  lowest  when  the  insulation  was  placed  on  the  outside  of  the  building. 

The  effect  of  an  indoor  mass  on  the  thermal  behavior  of  the  experi- 
mental structure  was  small.  For  tests  in  which  the  variation  of  the  in- 
side temperature  was  small  (such  as  the  thermostated  tests),  the  effect 
of  an  indoor  mass  was  practically  negligible. 

The  NBSLD  computer  program  was  experimentally  validated  for  pre- 
dicting the  daily  indoor  air  temperature  profile  as  it  is  influenced  by 
known  outdoor  temperature  conditions  and  the  effect  of  the  mass  and  thermal 
resistance  of  this  experimental  building.  Furthermore,  when  the  inside 
air  temperature  was  thermostated,  this  program  predicted  the  peak  and 
daily  average  heating  loads  and  may  therefore  be  used  to  size  equipment 
needed  to  condition  the  interior  of  a building  and  to  predict  energy 
requirements.  It  was  shown  that  steady-state  methods  of  heating  load 
calculation  could  result  in  oversizing  heating  equipment  by  30%  or  more. 

The  NBSLD  dynamic  method  takes  into  account  heat  storage  effects  and 
therefore  predicts  the  peak  heating  load  more  realistically.  The  maximum 
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difference  between  the  computer  calculated  peak  heating  load  and  measured 
values  was  six  percent  and  the  average  difference  was  3,5  percent  for  the 
five  tests. 
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1.  Introduction 


To  provide  a functional  and  habitable  indoor  environment  for  a 
building  requires  careful  consideration  of  the  properties  and  performance 
of  the  materials  that  cover  the  building  frame  together  with  careful  de- 
sign, specification  and  installation  of  its  mechanical  and  electrical 
systems.  The  building  materials  and  systems  taken  together  are  a major 
part  of  the  cost  of  a new  building  and  the  fuel  consumption  also  con- 
stitutes a substantial  long-term  expense  in  the  operation  of  a building. 

The  indoor  thermal  environment  of  a building  is  influenced  by  the 
weather,  by  the  thermal  behavior  of  the  walls,  roof  and  floors,  by  heat- 
producing  occupant-related  activities,  and  especially  by  the  mechanical, 
electrical  and  service  systems  that  must  function  to  provide  control  of 
the  heating  and  cooling  devices  that  serve  to  make  living  spaces  habitable. 

This  study  explores  the  actual  dynamic  or  time-variable  flow  of  heat 
into  and  out  of  the  fabric  of  a building  and  the  resulting  temperature 
patterns  of  the  indoor  air  and  the  structure  itself.  Present  practices 
are  based  largely  on  steady-state  assumptions  and  techniques.  The  actual 
performance  is  dynamic  because  of  the  changing  patterns  of  weather  and 
climate.  Therefore,  analysis  and  predictions  of  hourly,  daily  and  seasonal 
system  performance  should  be  based  on  dynamic  considerations.  The  theory 
and  basic  mathematics  for  the  dynamics  of  such  a system  were  first  ex- 
plained by  Fourier,  about  1820,  but  the  complexity  of  calculation  and 
the  time  and  expense  involved  has  deterred  architects  and  engineers  from 
using  such  sophisticated  procedures  to  design  and  evaluate  buildings. 
Simplified  steady-state  approaches  have  been  and  are  still  used  in  com- 
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bination  with  engineering  judgment. 

Design  calculations  for  the  heating  and  cooling  loads  for  buildings 
have  been  performed  virtually  by  a multiplicity  of  arithmetic  and  algebraic j 
computations.  It  was  not  practical  to  make  an  extensive  type  of  design 
analysis  and  the  loads  were  generally  determined  by  employing  simple 
equations  using  selected  fixed  winter  and  summer  design  temperatures. 
Experience  has  shown  that  systems  designed  on  this  basis  are  sometimes 
oversized  and  may  not  operate  at  full  load  and  optimum  efficiency. 

With  the  advent  of  high  speed  electronic  digital  computers  with  a 
large  memory  bank,  it  is  now  possible  to  make  a comprehensive  design 
analysis  which  includes  the  dynamic  performance  of  buildings  as  affected 
by  diurnal  and  seasonal  patterns  of  the  weather  and  the  time  dependent 
interactions  within  the  building  itself.  This  approach  allows  an  engi- 
neer to  rapidly  and  inexpensively  calculate:  (a)  energy  requirements 

with  consideration  of  operating  costs,  (b)  heating  and  cooling  load 
profiles  for  equipment  design  or  selection  and  operation,  (c)  the  in- 
formation that  will  permit  the  design  engineer  to  rapidly  evaluate  a 
large  number  of  options  in  the  design  process,  and  (d)  optimum  efficiency 
of  energy  utilization  which  is  becoming  increasingly  important  as  a 
national  concern. 

Computer  programs  usually  contain  approximations  that  require  ex- 
perimental validation  before  being  adopted  for  wide-scale  use.  In  ad- 
dition, the  performance  data  on  building  materials  and  elements,  design 
weather  data  and  boundary  conditions  at  surfaces  need  better  definition 
to  assure  accuracy  of  predicted  results. 
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It  is  the  objective  of  this  study  to  produce  a computer  program 
suited  to  the  variable  temperature  and  heat  flow  regimes  in  most  real 
situations  and  to  compare  results  as  predicted  by  this  program  with 
measurements  made  in  the  laboratory  on  full  scale  structures  that  are 
subjected  to  changing  simulated  weather  patterns.  Further,  it  was 
hypothesized  that  building  walls,  roofs  and  floors  can  be  better  de- 
signed to  take  advantage  of  thermal  lags  that  occur  due  to  the  mass  of 
the  building  and  thereby  allow  a reduction  of  the  installed  capacity 
of  mechanical  equipment  for  heating  and  cooling  while  still  maintaining 
performance  satisfactory  for  human  comfort  and  health.  For  example, 
it  was  hypothesized  that  if  the  masonry  of  a building  is  located  on  the 
indoor  side  of  walls  and  roofs  with  thermal  insulation  on  the  outside 
the  stability  of  indoor  temperature  changes  should  be  improved  with  less 
gross  energy  expended  for  maintaining  a selected  indoor  temperature 
level.  Also,  locating  masonry  on  the  inside  of  the  walls  with  insula- 
tion on  the  outside  provides  other  potential  advantages  such  as : a 

reduction  of  cracking  and  spalling  because  the  masonry  remains  unex- 
posed to  weather  and  at  essentially  a constant  temperature  and  moisture 
content;  the  use  of  strong  durable  indoor  surfaces  should  allow  a re- 
duction in  the  costs  of  maintenance  and  redecorating;  a possible  im- 
provement in  acoustic  performance;  and  a greater  resistance  of  the 
building  to  an  interior  fire,  or  its  rapid  spread.  When  compared  with 
the  usual  construction  of  walls  with  masonry  outside  and  insulation 
inside,  the  proposed  inverted  system  with  insulation  outside  has 
elicited  considerable  interest. 
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A concerted  effort  towards  the  experimental  verification  of  com- 
puter calculation  methods  and  the  technical  merits  of  the  inverted  sys- 
tem is  needed.  At  the  National  Bureau  of  Standards  the  initial  experi- 
mental phases  in  this  regard  included  laboratory  testing  in  a high-bay 
environmental  chamber  employing  a prototype  building  where  the  time 
varying  external  environment  could  be  controlled,  reproduced  and  varia- 
tions in  important  parameters  could  be  studied. 

This  report  presents  a computer  program  for  prediction  of  dynamic 
thermal  and  energy  loads  of  buildings,  the  experimental  results  obtained 
from  laboratory  measurements  made  on  a prototype  building  and  the  com- 
parison of  experimental  results  with  those  calculated  by  the  computer 
program.  In  conjunction  with  the  experimental  phases  involving  the 
dynamic  thermal  performance  of  a prototype  building,  two  other  significant 
experiments  were  performed  on  the  building.  The  first  experiment  was 
concerned  with  the  air  infiltration  rate  of  the  building.  The  method, 
procedure  and  results  are  contained  in  Appendix  A of  this  report.  The 
second  experiment  involved  a series  of  noise  transmission  measurements 
made  on  the  building.  The  method,  procedure  and  results  are  contained 
in  Appendix  B of  this  report.  Other  observations  included  monitoring  of 
the  moisture  content  of  the  prototype  building  and  the  movement  of  the 
walls  of  the  building  under  the  influence  of  the  changes  in  simulated 
outdoor  air  temperatures.  The  moisture  content  reached  low  equilibrium 
values  early  in  the  program  and  remained  stable  thereafter.  Wall  move- 
ment was  little  and  about  what  would  be  expected  using  predictive  engi- 
neering calculations.  No  surface  or  through-the-wall  cracks  in  masonry 
were  observed  at  any  time  in  the  program. 
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2.  Prediction  and  Evaluation  Analyses 


In  order  to  evaluate  the  dynamic,  rather  than  steady-state  thermal 
behavior  and  response  of  the  fabric  of  a building  as  affected  by  diurnal 
and  seasonal  variations  of  weather  and  the  time  dependent  interactions 
within  the  building,  it  was  necessary  to  make  a comprehensive  mathemati- 
cal analysis  of  the  various  heat  transfer  problems  and  translate  the 
derived  expressions  into  computer  programs.  It  was  found  that  the  heat 
conduction  portion  of  the  overall  problem  could  not  be  satisfied  by 
purely  rigorous  mathematical  solutions  to  the  applicable  partial  dif- 
ferential equations  because  some  of  the  boundary  conditions  at  solid 
surfaces  cannot  be  represented  in  a rigorous  form  in  a reasonable  man- 
ner. For  these  reasons,  the  Response  Factor  method  was  employed  for 
those  portions  of  the  problem  involving  heat  conduction,  because  it 
allows  a time  variation  of  boundary  conditions  and  can  readily  be  re- 
lated to  similar  and  other  modes  of  heat  flow,  such  as  radiation  and 
time  varying  changes  in  the  nature  of  convection  heat  flow. 
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Basically,  the  Response  Factor  method  predicts  one-dimensional  heat 
flow  by  utilizing  the  superposition  principle  in  such  a manner  that  the 
overall  thermal  response  of  a solid  at  a selected  time  is  the  sum  of  the 
responses  caused  by  many  individual  temperatures  or  heat  flux  pulses 
during  preceding  time  steps.  Thereby,  transient  boundary  conditions  are 
simulated  by  a train  of  pulses.  By  summing  up  the  fluxes  or  temperatures 
caused  by  each  pulse,  the  total  heat  flux  or  temperature  at  a given  time 
can  be  determined.  The  differential  equations  of  heat  conduction  for 
multi-layer  systems  of  a building  are  solved  in  this  method  by  employing 
matrix  equations  of  the  Laplace  transforms.  The  matrix  algebra,  super- 
position principle,  and  inversion  of  the  Laplace  transforms  are  shown 

•>v  j 

and  discussed  by  Kusuda~  . Experience  has  shown  that  when  this  method 
is  compared  with  a rigorous  analytical  solution  under  simplified  condi- 
tions, the  agreement  is  very  good,  except  for  the  case  where  sudden 
changes  or  amplitude  peaks  of  a weather  cycle  are  encountered.  This  is 
probably  due  to  the  time  steps  employed  and  is  not  considered  to  be  a 
serious  drawback. 


Thermal  Response  Factors  for  Multi-layer  Structures  of  Various 
Heat  Conduction  Systems,  American  Society  of  Heating,  Refrigerating 
and  Air-Conditioning  Engineers  (ASHRAE)  Transactions,  1969,  pp . 


246-271. 
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Appendix  C contains  the  complete  computer  program,  NBSLD,  Computer 
Programs  to  Obtain  Heating  and  Cooling  Loads  and  to  Estimate  Room  Air 
Temperature  Change  Using  Thermal  Response  Factors.  For  the  purpose 
of  predicting  performance  in  the  experiment  certain  subroutines  of 
NBSLD  were  not  needed.  Appendix  D is  the  computer  program  as  adapted 
from  NBSLD  for  use  in  this  report  for  comparing  predicted  results  with 
experimentally  measured  results.  Appendix  E gives  a sample  set  of  in- 
put and  a print-out  of  corresponding  computer  results  as  used  with  the 
program  of  Appendix  D. 

For  this  thermal  analysis,  the  following  assumptions  were  made: 

1.  The  conduction  heat  transfer  through  all  the  components 
of  the  experimental  structure  was  assumed  to  be  one-di- 
mensional . 

2.  All  building  materials  were  assumed  to  be  homogeneous 
having  constant  physical  and  thermal  properties  over 
the  operating  temperature  range  of  the  tests. 

3.  For  the  tests  considered  in  this  report  the  heat- 
transfer  coefficients  for  the  inside  and  outside 
surfaces  of  the  experimental  structure  were  assumed 
to  be  constant . 

4.  Heat  and  mass  transfer  of  water  in  vapor  or  liquid  form 
or  the  latent  heats  of  condensation  and  evaporation  were 
not  considered  in  the  analysis.  For  most  tests,  the  dew 
point  temperature  of  the  outside  air  was  maintained  below 
that  for  any  temperature  occurring  in  daily  cycle. 
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5.  Infiltration  of  air  from  the  outside  to  the  inside  and 


from  the  inside  to  the  outside  was  considered  to  be  a 
constant  for  a particular  test.  Two  tests  were  performed 
for  determining  the  air  infiltration  rates  of  the  build- 
ing, one  with  and  the  other  without  windows  installed. 

The  description  and  results  for  the  air  infiltration 
tests  is  in  Appendix  A. 

3.  Description  of  Building 

The  building  was  constructed  in  a high-bay  environmental  labora- 
tory of  approximately  70,000  cubic  feet  in  volume.  A photograph  of 
the  experimental  structure  located  in  the  environmental  chamber  is 
shown  in  figure  1.  In  this  laboratory  the  temperature  and  relative 
humidity  can  be  controlled  over  the  ranges  -50  to  150  °F  and  15  to  85 
percent,  respectively.  Temperatures  and  relative  humidities  can  be 
changed  as  a function  of  time  using  cam-operated  controllers.  The 
floor  of  the  laboratory  is  undisturbed  earth  suitable  for  placing 
building  foundations. 
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The  outside  plan  dimensions  of  the  building  were  20*  x 20*  with 
10  foot  high  walls.  The  flat  roof  consisted  of  five  20*  long  by  4*  wide 
and  4 inch  thick  steel  reinforced  concrete  roof  slabs  as  shown  in  figure 
2,  The  walls  were  made  of  nominal  8"  high  by  8"  wide  and  16"  long  solid- 
cinder  aggregate  concrete  blocks  joined  with  fully  bedded  mortar  joints. 
The  blocks  were  of  a nominal  100  pound  per  cubic  foot  density.  Eight 
concrete  lintels  were  installed  at  appropriate  locations;  one  above  each 
of  the  seven  windows  that  were  40"  high  and  32"  wide  and  one  above  the 
solid  wood  door  measuring  79"  high  x 32"  wide  x 2"  thick.  Window  open- 
ings were  filled  with  blocks  for  the  first  two  tests  (see  figure  2).  The 
blocks  were  removed  and  the  windows  installed  for  the  remaining  tests. 

The  exposed  glass  area  was  about  8 percent  of  the  exposed  wall  area  or 
about  18  percent  of  the  floor  area.  Figure  3 shows  the  configuration. 

A detailed  illustration  of  the  floor  and  the  footing  supporting 
the  walls  is  given  in  figure  4.  Below  the  ground  level,  four  inch 
thick  polystyrene  insulation  was  placed  on  the  outside  and  a one  inch 
thickness  on  the  inside  of  the  concrete  blocks  to  a depth  of  16  inches. 
Below  the  16  inch  depth  a one  inch  thickness  was  placed  on  the  outside 
of  the  footing.  The  floor  was  made  of  two  inches  of  polystyrene  insu- 
lation placed  on  the  earth  with  a two  inch  thick  concrete  slab  on  top  of 
the  insulation.  Considerable  insulation  was  purposedly  placed  below  grade 
to  reduce  the  known  long-term  influence  of  heat  flow  to  the  earth  from  the 
building  and  to  minimize  the  time  necessary  for  experimental  test. 
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Cracks  at  the  roof-wall  interface  and  between  the  roof  slabs  were 


caulked  with  a polysulfide  sealant.  When  the  windows  were  installed, 
all  cracks  including  those  at  the  glass-wood  frame  interface  were  also 
caulked  with  the  same  sealant.  Windows  were  as  shown  in  figure  4. 

Commercial  expanded  polystyrene  board-type  insulation  2 inches 
thick,  when  used,  was  spot  glued  to  either  the  inside  or  outside  sur- 
faces and  all  cracks  were  tape  sealed.  The  identical  insulation  was 
used  inside  and  outside.  An  internal  mass  consisted  of  2600  pounds 
of  solid  concrete  blocks  stacked  on  the  floor  as  shown  in  figure  3 was 
used  to  simulate  the  heat  capacity  effect  of  interior  partitions, 
furniture,  etc. 


4.  Instrumentation  and  Transducers 

Temperatures  were  measured  using  24  gage  copper-constantan  thermo- 
couples. The  dots  on  figure  5 indicate  thermocouple  locations.  The 
five  vertical  planes  A,  B,  C,  D and  E,  as  shown  on  the  plan  view  of 
figure  6,  each  contained  the  same  thermocouple  configuration  given 
in  figure  5,  except  for  the  indoor  air  thermocouples  which  were  located 
only  in  the  vertical  plane  B.  Four  thermocouples  were  placed  in  the 
air  one  foot  from  the  outside  surfaces.  One  of  these  was  located  at 
the  center  of  the  roof  and  the  other  three  were  located  at  the  mid- 
height of  the  three  walls  denoted  by  vertical  planes  B,  D and  E of  fig- 
ure 6 . 
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Six  heat  flow  meters  were  placed  on  inside  surfaces,  five  of  them 
in  vertical  plane  B of  figure  6.  One  was  placed  at  the  center  of  the 
floor  and  a second  meter  was  placed  on  the  floor  at  a distance  two  feet 
in  from  the  wall.  Two  meters  were  placed  on  the  ceiling  opposite  those 
on  the  floor.  The  fifth  meter  was  placed  on  wall  at  mid-height.  The 
sixth  meter  was  placed  mid-height  on  the  wall  of  vertical  plane  D. 

The  heat  flow  meters  were  circular  disks  2.0  in.  in  diameter  and  0.13 
in.  thick,  made  of  tan  polyvinylchloride  filler  material,  each  having 
an  embedded  spiral  of  he li cal ly -wound  wire  comprising  a large  number 
of  thermo junctions  in  series  (with  internal  resistance  range  of  135  to 
170  Q)  distributed  over  a circular  area  15/8  inches  in  diameter  located 
centrally  in  the  disk.  Two  wires  attached  in  each  meter  acted  as  leads 
for  the  series  thermopile  of  the  meter.  The  meters  were  calibrated  in 
an  8 in.  guarded  hot  plate  apparatus  conforming  with  the  requirements 
of  Standard  Method  of  Test  ASTM  C177. 

All  thermocouple  and  heat  flow  meter  leads  were  connected  to  ther- 
mally isolated  terminal  strips  at  the  center  of  the  room  from  which  cop- 
per leads  went  to  a data  acquisition  system.  The  terminal  strips  were 
mounted  on  a one -quarter  inch  thick  aluminum  plate  which  in  turn  was 
surrounded  by  three  inches  of  polyurethane  insulation.  All  lead  wires 
were  surrounded  by  three  inches  of  the  same  insulation  for  a distance 
of  seven  inches.  This  assembly  is  termed  a zone  box.  Four  additional 
thermocouple  leads  were  connected  to  the  terminal  strips  at  ends  of  the 
zone  box  and  their  junctions  were  placed  in  an  ice  point  reference  ex- 
ternal to  the  building.  The  readings  from  these  four  thermocouples 
gave  the  temperature  of  the  zone  box  as  a reference  temperature  for 
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the  other  thermocouple  leads. 

Copper  leads  from  the  zone  box  were  connected  to  terminals  of  the 
data  acquisition  system  which  converted  the  analog  signals  to  digital 
information  which  in  turn  was  recorded  on  punched  cards. 

Electric  power,  when  supplied,  to  the  building  was  measured  using  a 
calibrated  single-phase  watthour  meter  equipped  with  an  impulse  generator 
The  impulse  generator  is  a photo-electric  device  which  counts  the  revo- 
lutions of  the  disk  inside  the  watthour  meter.  A digital  signal  (number 
of  revolutions  of  the  disk)  was  fed  into  the  data  acquisition  system 
which  in  turn  recorded  the  digital  signal  on  punched  cards  at  selected 
time  intervals. 


5.  Experimental  Procedure 

Figure  7 is  a representative  sample  of  the  outside  air  temperature 
wave-form  imposed  on  the  structure  for  each  24-hour  time  period.  The 
limits  40  °F  to  100  °F  were  selected  for  experimental  convenience  and 
because  their  average  would  be  approximately  a normal  room  temperature. 
The  curve  is  the  average  of  the  four  individual  temperatures  indicated 
by  thermocouples  in  the  air  one  foot  from  the  exterior  surface  of  the 
structure.  The  maximum  difference  in  temperature  between  any  of  these 
four  locations  was  always  less  than  4 °F.  The  outside  dew  point  tem- 
perature was  maintained  constant  at  approximately  5 °F  below  the  lowest 
temperature  of  a cycle.  The  temperature  cycle  of  figure  7 was  selected 
as  a simulated  sol-air  temperature  pattern  as  given  in  Table  25,  page 
490  of  the  "Handbook  of  Fundamentals",  published  by  the  American  Society 
of  Heating,  Refrigerating  and  Air-Conditioning  Engineers,  1967.  Sol- 
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air  temperatures  were  area  averaged  for  orientations  north,  east,  south, 
west  and  horizontal.  The  temperature  variation  as  indicated  on  figure 
7 was  maintained  for  a period  of  from  three  to  four  days  before  a final 
set  of  data  was  taken.  This  conditioning  period  was  deemed  to  be 
necessary  and  sufficient  to  eliminate  transient  heat  flows  thereby 
giving  only  those  heat  flows  that  would  occur  in  a steady-periodic 
condition. 

A complete  set  of  data  for  each  test  consisted  of  recording  the 
digital  output  from  analog  signals  of  171  sensing  elements  (thermocouples, 
etc.)  every  30  minutes  for  a 24  hour  period.  The  recorded  data  on 
punched  cards  was  fed  into  the  computer  programmed  to  process  the  data 
into  temperatures,  heat  flows,  etc.  The  converted  data  were  then 
transferred  to  magnetic  tape  for  use  in  analyses,  and  plotting  as 
temperature  and  heat  flow  patterns. 

The  results  from  ten  tests  given  in  this  report  are  derived  from 
the  five  floating  tests  and  five  thermostated  tests  summarized  in 
Tables  1 and  2» 


a«  Floating  Tests 

Floating  tests  are  defined  as  those  tests  where  no  heat  energy 
was  added  or  taken  away  from  the  interior  air  of  the  experimental 
structure  by  mechanical  equipment.  The  temperature  of  the  interior 
air  was  allowed  to  "float"  or  respond  to  changes  in  the  outside  air 
temperature.  Five  floating  tests  were  conducted  with  variations  in 
test  conditions  as  shown  in  Table  1. 
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Table  1 


Floating  Tests 


Test  No. 

Insulation 

Windows 

Internal  ! 

1 

None 

None 

None 

2 

None 

None 

Mass* 

3 

None 

Single  Pane 

None 

4 

Inside 

Single  Pane 

Mass* 

5 

Outside 

Single  Pane 

Mass* 

2600  lbs  of 

concrete  blocks 

• 

b.  Thermostated  Tests 

Thermostating  tests  are  defined  as  those  tests  where  heat  energy 
was  added  to  the  interior  air  of  the  experimental  structure  by  four 
electric  heaters  under  thermostatic  control.  The  variations  in  test 
conditions  are  shown  in  Table  2 along  with  the  average  inside  air  tem- 
perature maintained  and  its  root  mean  square  deviation. 


Table  2 

Thermostated  Tests 


Inside  Air 

Test  No, Insulation  Windows Internal  Mass Temp , 


Single 
Single 
Single 
Double 
Double 

•k 


Pane 

None 

78.9 

+ 

1.2 

Pane 

Mass* 

76.9 

+ 

0.8 

Pane 

Mass* 

77.6 

+ 

0.6 

Pane 

Mass* 

77.6 

+ 

0.6 

Pane 

Mass* 

74.2 

+ 

0.8 

6 None 

7 Inside 

8 Outside 

9 Outside 

10  Outside 


2600  lbs  of  concrete  blocks 


The  sensing  element  for  thermostating  was  a thermocouple  placed  in  the 
middle  of  the  room  at  mid-height.  It  controlled  the  operation  of  four 
fan  heaters  placed  as  shown  in  figure  3 in  an  on-off  type  of  control 
with  a differential  of  approximately  +2  ®F.  Each  drum-type  fan  heater, 
as  shown  in  figure  8,  consisted  of  a 600  watt  cone  heater  and  a blower 
which  takes  air  from  the  floor  level  passes  it  through  the  heater  cham- 
ber and  into  the  room  through  peripheral  holes  near  the  top  of  the  drum. 
For  test  10,  the  daily  temperature  cycle  for  the  outside  air  ranged  from 
10  to  70  °F,  but  the  cycle  was  identical  in  shape  to  that  given  in 
figure  7. 


6.  Results  and  Discussion 


The  thermal  and  physical  properties  of  the  materials  comprising 
the  building  which  are  necessary  for  use  in  the  computer  program  are 
given  in  Table  3. 


Table  3 


Thermal  and  Physical  Properties 


Thickness 

Thermal 

Conductivity 

Density 

Specific  Heat 

in. 

Btu  hr~^  ff^  F“^ 

lbs  ft“^ 

Btu  lb"l  F’l 

Concrete 

7.5 

.29 

100 

.18 

Block 

Roof  Slab 

4. 

.80 

150 

.2 

Polystyrene 

2. 

.018 

2.5 

.27 

Insulation 

Concrete 

2. 

.80 

150 

.2 

Floor 

Earth 

.5 

120 

.2 
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Measurements  of  thermal  conductivity,  thickness  and  density  were  made  on 
oven-dried  samples  of  the  concrete  block  and  polystyrene  insulation  in 
accordance  with  the  hot  plate  method  given  in  ASTM  C177.  All  other 
properties  were  obtained  from  available  literature. 

The  coefficients  of  heat  transfer  at  the  inside  and  outside  sur- 
faces were  the  most  difficult  of  the  numerous  parameters  to  define  for 
this  experimental  work.  Values  given  in  literature  are  usually  determined 
from  steady-state  conditions  whereas  the  test  conditions  were  dynamic 
and  the  coefficients  vary  with  orientation  of  surfaces,  direction  of 
heat  flow,  temperature  of  surface  and  the  air  motion  over  the  surface. 

During  one  of  the  tests,  an  attempt  was  made  to  measure  air  veloci- 
ties at  inside  and  outside  surfaces  of  the  building  with  a vane  anemometer. 
The  air  velocities  were  not  sufficient  to  rotate  the  vanes  indicating 
that  the  velocities  were  somewhat  less  than  50  fpm  and  that  conditions 
at  the  surfaces  could  be  considered  as  natural  convection.  Under  natural 
convection  conditions  the  convection  component  of  the  heat  transfer  co- 
efficients is  defined  in  literature  as  being  proportional  to  the  one- 
third  power  of  the  absolute  temperature  difference  between  a surface  and 
the  adjacent  air  if  in  the  turbulent  range.  This  relationship  may  apply 
to  vertical  surfaces,  heated  horizontal  surfaces  facing  upward  and  cooled 
horizontal  surfaces  facing  downward.  For  horizontal  surfaces  either 
heated  facing  downward  or  cooled  facing  upward,  the  adjacent  air  is  con- 
sidered to  be  in  the  laminar  range  and  the  convection  component  becomes 
very  small  and  the  radiation  component  of  the  heat  transfer  coefficient 
becomes  dominant . 
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From  the  above  considerations,  values  for  the  coefficients  of  heat 
transfer  at  the  various  surfaces  were  selected  and  used  in  the  computer 
program  as  constants  for  the  time  period  of  a test.  The  coefficients 

used  for  the  inside  surfaces  at  the  ceiling,  walls,  and  floor  were  1.08, 

- 1 -2  - 1 

1,1,  and  1.08  Btu  hr  ft  F , respectively.  The  heat  transfer  coef- 

- 1 -2  - 1 

ficient  for  the  outside  surfaces  was  selected  to  be  1.47  Btu  hr  ft  F . 
In  general  these  values  are  based  on  a value  of  0.9  for  the  radiation 
component  of  heat  transfer  and  time  averaged  temperature  differences  be“ 
tween  surfaces  and  adjacent  air  of  1 F and  14  F for  the  inside  and  the 

outside,  respectively.  For  test  10,  where  the  outside  temperature  was 

-1  -2  -1 

considerably  lower,  the  coefficient  selected  was  3 Btu  hr  ft  F . 
Reasonable  variations  in  these  coefficients  show  a negligible  effect  on 
results  from  the  computer  program. 

For  the  computer  program,  the  heat  capacity  effects  of  the  door  and 
windows  were  assumed  to  be  negligibly  small  and  only  the  thermal  re- 
sistance of  these  components  was  used.  For  the  door  and  single  and 

double  pane  windows  the  overall  coefficients  of  heat  transfer  were 

-1  -2  -1 

calculated  to  be  0.25,  0.45,  and  0.39  Btu  hr  ft  F , respectively, 
for  the  conditions  of  tests  1 through  9.  For  the  double  pane  windows 
of  test  10  the  selected  coefficient  was  0.46. 
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For  heat  flow  to  or  from  the  floor,  the  underlying  earth  was  con- 
sidered to  be  a one-dimensional  semi-infinite  medium  for  the  Response 
Factor  program,  and  the  average  of  temperatures  measured  at  the  one-foot 
depth  in  the  earth  was  used  as  the  earth  temperature  at  a depth  con- 
siderably removed  from  the  floor.  For  the  duration  of  the  tests  this 
was  deemed  an  adequate  assumption  because  the  root  mean  square  devia- 
tion of  the  earth  temperatures  at  the  one-foot  level  was  less  than  0.2  °F 
for  all  tests  where  the  diurnal  outside  air  temperature  varied  from  40 
to  100  °F  and  less  than  0.3  ®F  for  the  10  to  70  ®F  cycle.  For  the  40  to 
100  °F  tests,  the  average  temperature  at  the  top  of  the  footing  (figure 
5)  was  about  r5  F lower  than  the  earth  temperature  at  the  one-foot  level, 
and  for  the  10  to  70  °F  test  was  about  2 F lower.  This  indicates  that 
some  of  the  heat  is  flowing  from  the  earth  underlying  the  floor  toward 

the  footing  and  was  not  accounted  for  in  the  one-dimensional  heat  transfer 
approach  of  the  Response  Factor  program.  The  error  due  to  this  heat  flow 
is  believed  to  be  very  small  in  relation  to  other  heat  flows.  A mathe- 
matical analysis  was  performed  for  the  heat  flow  at  the  ground  level  in 
the  wall  section  below  ground  level  to  the  top  of  the  footing  (figure  4) 
using  the  temperature  variations  with  time  from  thermostated  test  6.  The 
computed  heat  flows  showed  that  heat  was  flowing  into  and  out  of  this 
section  with  time,  but  the  magnitude  of  these  heat  flows  was  small  in 
relation  to  other  heat  flows. 
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Air  infiltration  rates  were  determined  by  a tracer  gas  method  using 
helium  as  the  tracer  gas.  (See  Appendix  A.)  For  the  building  without 
and  with  windows,  measured  values  were  0.06  and  0.38  air  changes  per  • 
hour,  respectively.  Since  there  was  little  air  movement  at  the  inside 
and  outside  surfaces,  the  thermal  head  (the  difference  in  temperature 
between  the  inside  and  outside  air)  is  the  predominant  driving  force  for 
air  infiltration.  The  above  values  are  considered  maximum  rates  for  air 
infiltration,  because  the  tests  were  performed  when  the  thermal  head  was 
the  greatest.  It  would  be  expected  that  the  air  infiltration  rate  would 
hj  proportional  to  the  thermal  head.  For  this  reason,  average  air  infil- 
tration rates  were  selected  as  being  2 cfm  for  tests  with  no  windows  and 
10  cfm  for  test  with  windows.  The  tests  of  Appendix  A were  performed 
on  the  building  without  thermal  insulation.  Placing  insulation  on  either 
the  inside  or  outside  surfaces  would  increase  the  resistance  to  air  in- 
filtration. For  this  reason,  a rate  of  5 cfm  was  used  for  tests  with 
insulation. 
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Noise  reduction  measurements  were  made  on  the  prototype  building 
as  given  in  Appendix  B.  The  results  for  conditions  of  no  windows,  single 
pane  windows  only,  single-pane  windows  with  insulation  inside  and  single- 
pane windows  with  insulation  outside  are  shown  in  figure  2 of  Appendix 
Bo  As  indicated  and  expected  the  noise  reduction  was  greatest  without 
windows  and  some  improvement  is  shown  when  insulation  was  applied.  Com- 
parison of  noise  reduction  measurements  for  the  tests  with  insulation  on 
the  inside  and  on  the  outside  indicate  that  insulation  on  the  inside  had 
better  characteristics  because  of  the  higher  noise  reduction  values  in 
the  range  from  500  to  2500  Hz.  This  range  is  considered  to  contain  the 
most  objectionable  portion  of  the  audible  frequency  spectrum. 

As  mentioned  in  the  introduction  an  equilibrium  moisture  content 
of  the  block  was  rapidly  achieved  and  the  influence  of  moisture  in  these 
tests  is  considered  to  be  negligible.  The  moisture  content  of  the  block 
was  monitored  by  observing  the  change  of  weight  of  a single  oven  dried 
concrete  block  placed  in  the  environmental  chamber  and  the  test  room 
throughout  the  tests « The  equilibrium  moisture  content  of  the  single 
block  was  47o  by  weight.  Similarly,  vertical  and  horizontal  thermal  ex- 
pansion of  the  concrete  block  wall  attained  an  equilibrium  range  and 
was  considered  to  be  desirably  low  especially  since  no  surface  or  through 
the  wall  cracks  were  visible. 
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a.  Floating  Tests 


For  the  floating  tests  numbered  1 through  5,  Table  1,  the  measured 
and  computer  calculated  inside  air  temperatures  are  plotted  in  figures 
9 through  13,  respectively,  each  with  its  measured  outdoor  air  temperatures « 
The  curve  of  measured  indoor  air  temperature  is  the  arithmetic  average  of 
the  six  indoor  air  thermocouples  as  shown  in  figure  5.  The  vertical  dis- 
tribution of  temperature  within  the  room  will  be  treated  later  in  this 
discussion.  There  was  generally  good  agreement  between  the  measured  and 
predicted  average  inside  air  temperatures  in  all  cases,  although  there 
is  a trend  for  the  predicted  values  to  have  slightly  higher  maximum  values 
and  lower  minimum  values  during  the  24-hour  cycle.  This  indicates  that 
the  mass  of  the  building  dampens  temperature  changes  more  than  is  accounted 
for  in  the  predictive  computer  program.  This  may  be  due  to  several  factors 
such  as  the  theoretical  model  that  was  used  in  the  computer  programs  neg- 
lects the  additional  thermal  inertia  introduced  at  the  corners  of  the 
building  and  neglects  slight  changes  in  material  physical  properties 
during  exposure  as  compared  with  measured  dry  values. 
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Comparing  the  indoor  air  temperature  curves  of  figures  9 through 
13,  it  can  be  seen  that  placing  insulation  on  inside  and  outside  building 
surfaces  had  a marked  influence  on  the  inside  air  temperature  profiles. 
Compare  figure  12  with  figure  11,  and  figure  13  with  figure  11.  The 
temperature  deviations  from  the  daily  average  inside  air  temperature  for 
the  building  with  windows  are  plotted  in  figure  14  for  the  cases  of  no 
insulation,  insulation  on  the  inside  building  surface,  and  insulation  on 
the  outside  building  surface  (corresponding  to  figures  11,  12  and  13). 
Adding  insulation  on  the  inside  surface  of  the  building  reduced  the  peak 
to  peak  variations  of  the  inside  air  temperature  from  about  10„5  to  5.5  °F. 
The  effect  of  the  insulation  then  was  to  damp  out  the  cyclic  fluctuations 
of  inside  air  temperature  with  windows  installed o Furthermore,  when 
insulation  was  placed  on  the  outside  surfaces,  the  peak  to  peak  variation 
was  reduced  to  about  2 °F.  This  experimental  finding  is  considered  to 
be  significant  because  no  heat  energy  was  purposedly  added  to  or  taken 
away  from  the  indoor  air  during  the  tests  and  the  performance  results 
illustrate  that  considerable  control  of  the  indoor  air  temperature  can 
be  exercised  by  simply  placing  the  mass  of  the  walls  and  roof  facing 
indoors  with  insulation  facing  the  outdoors. 
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To  investigate  the  effect  of  an  interior  mass  on  the  inside  air  tem- 
perature for  a floating  test,  a comparison  was  made  between  tests  1 and 
2,  (figures  9 and  10).  Temperature  deviations  from  the  measured  mean  in- 
side air  temperature  for  these  two  tests  are  plotted  in  figure  15.  It 
can  be  seen  that  for  these  cases  with  no  insulation,  the  presence  of  an 
internal  mass  slightly  damps  the  inside  air  temperature  cycle.  This  ef- 
fect was  also  predicted  by  the  Response  Factor  program.  For  floating 
tests  with  either  insulation  on  the  inside  or  outside  surfaces  (figures 
10  and  11),  the  effect  of  an  internal  mass  is  reduced  to  negligible  pro- 
portions. This  is  because  the  heat  absorption  and  rejection  by  the  in- 
ternal mass  is  very  small  when  the  cyclic  fluctuations  of  the  inside  air 
temperature  are  small. 

To  examine  the  effect  of  windows  on  the  inside  air  temperature,  a com- 
parison was  made  between  tests  numbered  1 and  3.  The  measured  temperature 
deviations  from  the  mean  inside  air  temperature  for  these  two  tests  are 
plotted  in  figure  16.  From  figure  16  it  may  be  seen  that  for  the  two  cases 
without  insulation  the  effect  of  adding  windows  had  little  effect  on  the 
cyclic  fluctuations  of  the  inside  air  temperature.  The  percent  glass  to 
wall  area  was  8.4.  For  cases  with  insulation,  one  would  expect  the  addi- 
tion of  windows  would  have  a more  pronounced  effect  on  the  cyclic  fluctu- 
ations of  the  inside  air  temperature,  since  the  heat  flow  through  windows 
would  be  a larger  percentage  of  the  total  heat  flow.  Direct  experimental 
comparison  is  not  possible  because  measurements  were  not  made  on  the 
structure  with  insulation  either  inside  or  outside  without  windows.  For 
practical  purposes  an  improvement  in  the  indoor  temperature  profile  as  shown 
in  figure  13  by  elimination  of  windows  is  considered  to  be  negligible. 
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Figures  17  and  18  show  the  inside  and  outside  wall  surface  tempera- 
ture variations  for  test  1;  no  insulation,  no  windows  and  no  internal 
mass.  Each  curve  represents  the  average  temperature  of  five  thermo- 
couples located  at  the  same  height  above  the  floor  and  at  the  wall 
positions  as  shown  in  figure  6.  From  these  graphs  it  can  be  seen  that 
the  inside  and  outside  wall  surface  temperatures  for  this  floating  test 
differ  from  each  other  within  a 2 ®F  band  except  for  the  average 
temperatures  at  the  0 and  10  foot  levels.  This  suggests  that  the  as- 
sumption of  one-dimensional  heat  transfer  is  valid  over  a major  area 
of  the  wall  surface,  multi-dimensional  effects  being  confined  in  a re- 
gion near  the  junctures  of  the  floor  to  wall  and  the  roof  to  wall.  Fig- 
ures 17  and  18  also  show  by  comparison  the  effect  of  thermal  resistance 
and  mass  of  the  building,  i.e.,  at  the  10  foot  level  the  outside  surface 
changed  in  temperature  by  about  30  °F  while  the  inside  surface  at  the 
same  level  changed  by  about  16  °F.  Also,  the  highest  and  lowest  tem- 
peratures on  the  outside  surface  occurred  about  2 hours  sooner  than  the 
inside  surface.  The  use  of  thermal  insulation  resulted  in  a much  more 
uniform  inside  wall  temperature  distribution.  For  instance,  when  in- 
sulation was  placed  on  the  outside  surface  of  the  building  (floating 
test  numbered  5)  a maximum  inside  wall  surface  temperature  fluctuation 
of  2.3  °F  occurred  over  the  24-hour  cycle  at  the  juncture  of  the  wall 
and  the  ceiling.  In  addition,  at  any  instant  the  maximum  floor  to 
ceiling  temperature  difference  along  the  inside  wall  surface  was  1.8  °F, 
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Comparisons  between  the  measured  and  calculated  heat  fluxes  at  the 
inside  surfaces  of  major  building  components  for  floating  test  1 are 
shown  in  figures  (19)  through  (21)  where  negative  values  denote  heat 
flow  into  the  room.  Measured  heat  fluxes  shown  for  the  floor,  roof, 
and  the  wall  were  obtained  using  heat  flow  meters  located  at  the  center 
of  the  floor,  the  center  of  the  roof,  and  at  the  midpoint  of  wall  in 
plane  B,  respectively,  (see  figure  6).  Since  both  the  measured  and  cal- 
culated data  contain  many  small  fluctuations  due  to  local  variations  of 
the  inside  air  temperature,  it  was  necessary  to  apply  a harmonic  analysis 
to  each  set  of  heat  flux  data,  maintaining  only  the  first  eight  terms  to 
give  the  smoothed  curves  shown  in  the  graphs.  From  figures  19  and  20  it 
can  be  seen  hat  the  agreement  between  the  measured  and  calculated  heat 
fluxes  at  the  inside  surface  for  floor  and  the  roof  was  very  good. 

Figure  21  shows  fairly  large  deviations  for  the  smoothed  measured  wall 
heat  flux  from  the  calculated  values.  The  same  type  of  performance  char- 
acteristic was  obtained  in  other  tests  where  the  floor  and  ceiling  also 
showed  good  agreement.  The  calculated  heat  flux  values  were  completed 
assuming  a constant  film  resistance  for  all  heat  flow  conditions.  From 
previous  discussion  concerning  the  film  resistance,  it  will  not  be  a con- 
stant, but  will  be  a function  of  the  heat  flow  conditions  that  promote 
air  flow  adjacent  to  the  surfaces.  Heat  flow  meters  are  very  sensitive 
instruments  and  the  signals  from  them  can  vary  considerably  when  subjected 
to  the  turbulent  air  motion  along  a wall.  Readings  that  were  taken  at  one 
instant  and  at  half-hour  intervals  would  not  be  expected  to  give  a true 
representation  of  average  signal  from  the  meters  for  the  time  period  under 
consideration.  For  measuring  heat  flow  at  wall  surfaces  the  signal  from 
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the  meters  should  have  been  recorded  for  finite  time  periods  to  give  more 
representative  values. 

To  study  the  processes  which  combine  together  to  produce  the  thermal 
performance  of  the  air  inside  a building,  the  profiles  of  the  heat  flux 
at  the  separate  inside  surfaces  during  the  outside  air  temperature  cycle 
were  plotted.  Figure  22  shows  the  variations  of  the  heat  flux  at  the 
inside  surfaces  of  the  roof,  walls,  floor  and  window  for  the  case  of  no 
insulation  (floating  test  3).  The  heat  flux  profiles  appearing  in  these 
graphs  were  calculated  by  the  Response  Factor  computer  program.  Positive 
values  signify  heat  flow  in  a direction  from  the  inside  to  the  outside. 

The  net  heat  transfer  to  or  from  the  indoor  air  at  any  instant  of  time 
is  equal  to  the  algebraic  sum  of  the  products  of  the  heat  fluxes  at  the 
surfaces  and  their  respective  areas  plus  the  heat  exchange  resulting  from 
air  infiltration.  For  the  floating  tests  this  sum  should  be  equal  to 
zero.  The  heat  flux  at  the  inside  surface  is  affected  by  the  resistance 
to  heat  flow  and  the  thermal  heat  capacity  of  the  materials  across  which 
heat  must  flow  to  the  surface  as  well  as  the  dynamic  conditions  of  the 
temperature  of  the  outside  and  inside  air.  For  this  reason,  heat  is 
simultaneously  flowing  out  of  and  into  different  surfaces  of  the  room. 

The  heat  flow  at  the  windows  is  in  phase  with  the  temperature  potential 
created  by  the  difference  in  the  outside  and  inside  air  temperature  be- 
cause heat  storage  (mass)  of  the  windows  was  negligible.  The  roof  and 
walls  are  not  in  phase  with  this  potential  due  to  their  appreciable  heat 
storage  capacity  and  their  minimum  values  (maximum  heat  flows  into  the 
room)  lag  behind  that  for  the  windows  by  about  3 and  9 hours,  respectively. 
The  roof  was  approximately  one-half  the  thickness  of  the  walls,  and  a 
smaller  delay  time  to  reach  a maximum  or  minimum  was  expected.  Heat  flow 
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into  and  out  of  the  floor  was  approximately  in  phase  with  the  inside  air 
temperature  cycle  shown  in  figure  11.  This  was  as  expected  because  the 
ground  temperature  beneath  the  floor  was  relatively  constant  with  time. 

A similar  analysis  of  heat  flow  was  performed  for  the  case  of  in- 
sulation placed  on  the  outside  surfaces  (floating  test  5).  Figure  23 
shows  the  profiles  of  the  heat  flow  at  the  inside  surfaces  for  this  test 
condition.  With  the  peak  outside  air  temperature  at  the  fourteenth  hour, 
the  delay  times  for  maximum  heat  flows  into  the  room  were  12  and  5 hours 
for  the  walls  and  roof,  respectively.  The  effect  of  placing  insulation 
on  the  outside  surface  was  to  increase  the  delay  time  (9  and  3 versus  12 
and  5)  and  considerably  reduce  the  amplitude  of  the  heat  flux  profiles. 

Figure  24  is  a plot  of  deviations  of  the  inside  air  temperature  from 
the  instantaneous  average  of  the  six  air  thermocouple  locations  shown  in 
figure  5 over  a twenty-four  hour  period  for  the  case  of  no  insulation 
(floating  test  3).  As  in  all  previous  plots,  the  peak  outside  air  tem- 
perature occurred  at  the  fourteenth  hour.  Positive  deviations  signify 
that  the  air  temperature  at  that  location  was  higher  than  the  average 
inside  air  temperature.  On  a daily  average  the  air  adjacent  to  the 
ceiling  was  about  2 °F  warmer  than  the  air  layer  next  to  the  floor  with 
the  floor  being  as  much  as  3 °F  warmer  and  8.5  °F  colder  than  the  ceiling 
during  portions  of  the  cycle.  The  portion  of  the  cycle  with  the  largest 
floor  to  ceiling  temperature  difference  (about  hour  18)  shows  a good 
example  for  a heated  surface  facing  downward  (ceiling)  and  a cooled  sur- 
face facing  upward  (floor)  where  the  air  flows  adjacent  to  the  two  sur- 
faces were  in  the  laminar  range  thus  producing  little  mixing  of  air  and 
large  vertical  temperature  gradients.  Conversely,  the  portion  of  the 
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cycle  with  the  smaller  temperature  differences  (about  hour  5)  shows  an 
example  for  a cooled  surface  facing  downward  (ceiling)  and  a heated  sur- 
face facing  upward  (floor)  where  the  air  flows  adjacent  to  the  surfaces 
were  in  the  turbulent  region  producing  mixing  of  air  by  natural  convection 
and  smaller  vertical  temperature  gradients.  One  must  conclude  from  figures 
22  to  24  that  the  indoor  convection  pattern  is  continually  changing,  as 
well  as  surface  coefficients  of  heat  transfer.  The  same  observations  can 
be  made  from  the  plots  of  deviations  from  the  average  inside  air  tempera- 
tures given  in  figures  25  and  26  for  insulation  placed  on  the  inside  (test 
4)  and  the  outside  (test  5)  surfaces,  respectively.  In  these  two  cases 
the  vertical  temperature  gradients  are  considerably  dampened  due  to  the 
addition  of  insulation,  and  subsequent  reductions  in  variations  of  the 
surface  temperatures. 


b.  Thermostated  Tests 

For  the  thermostated  tests  the  inside  room  air  temperature  was  main- 
tained within  an  approximate  2 °F  band  by  controlling  the  heat  input  to 
the  experimental  structure.  The  room  air  temperature  was  obtained  by 
averaging  the  six  air  temperatures  (figure  5)  at  each  time  interval. 

Figures  (27)  through  (31)  are  graphs  for  tests  numbered  6 to  10 
which  compare  the  measured  power  supplied  to  the  electric  heaters  and 
the  heating  load  calculated  by  the  Response  Factor  program  over  the 
24-hour  outdoor  air  temperature  cycle  as  shown  in  each  figure.  The 
calculated  load  was  computed  by  summing  the  net  heat  flows  through 
each  building  component  and  heat  flow  due  to  air  infiltration  at  each 
time  interval.  Areas  used  for  computing  heat  flows  were  the  arithmetic 
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averages  of  the  inside  and  outside  areas  of  each  building  component » 

Since  both  the  measured  and  calculated  heating  load  data  contained  many 
small  fluctuations  due  to  variations  of  the  inside  air  temperature,  it 
was  necessary  to  apply  a harmonic  analysis  to  each  set  of  heat  load  data. 
Only  the  first  eight  terms  were  maintained  to  give  the  smoothed  curves 
shown  in  the  graphs. 

As  shown  on  figures  27  through  31  the  minimum  measured  and  calculated 
heating  load  usually  occurred  later  in  the  day  than  the  peak  outside  air 
temperature  (hour  14)  because  of  the  effect  of  the  mass  of  the  building 
and  insulation  retarding  heat  flow  through  building  components.  Comparing 
the  cases  without  and  with  insulation,  figure  27  with  figures  28,  29,  30 
and  31,  it  can  be  seen  that  the  effect  of  placing  insulation  on  either 
the  inside  or  outside  surfaces  of  the  building  was  to  substantially  reduce 
the  amount  of  heating  needed  to  maintain  a constant  inside  air  temperature. 
Generally  the  correlation  between  computer  prediction  and  the  measured 
heating  load  profiles  is  reasonably  good.  There  was  less  than  a six 
percent  difference  between  the  maximum  computed  and  measured  heating 
loads  for  all  cases.  The  average  difference  was  3.5  percent  for  the 
five  tests. 

For  test  10  (figure  31)  and  test  9 (figure  30)  the  building  was  i- 
dentical  but  for  test  10  the  outdoor  temperature  cycle  was  changed  from 
40  - 100  F to  10  - 70  F and  the  indoor  air  temperature  was  changed  from 
77.6  F to  73.8  F.  The  shape  of  the  heating  load  profiles  are  similar 
but  for  test  9 the  maximum  and  minimum  loads  were  about  2600  and  500 
Btu/hr,  respectively,  and  for  test  10  about  6300  and  3600  Btu/hr,  re- 
spectively. The  maximum  loads  for  both  tests  are  lower  than  the  values 
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that  would  be  estimated  on  the  basis  of  steady-state  procedures  as  is 
discussed  later  in  this  paper  in  detail. 

For  the  thermostated  tests  with  insulation  (tests  7 through  10), 
the  measured  heating  loads  lag  the  calculated  heating  loads  over  part 
of  the  24-hour  cycle.  Consistently,  the  phase  lag  occurred  on  the  pro- 
files in  the  time  period  between  the  maximum  and  minimum  loads.  Also, 
some  phase  lags  occurred  following  the  minimum  loads.  The  reasons  for 
these  phase  lags  are  not  obvious  because  the  phase  lags  varied  from  one 
test  to  the  other  and  the  lag  is  especially  evident  in  test  9,  figure  30. 
It  was  found  during  analysis  that  the  calculated  heating  load  was  in- 
fluenced by  whether  the  inside,  outside  or  average  area  was  used,  lack 
of  heat  flow  allowance  for  corners  and  the  building  foundation,  varia- 
tions of  inside  air  temperature,  and  heat  transfer  coefficients  at  the 
inside  and  outside  surfaces. 

To  illustrate  the  effect  of  windows  on  the  thermal  behavior  of  the 
experimental  house,  calculations  were  made  using  the  Response  Factor 
method  for  the  cases  of  7 single  pane  windows,  7 double  pane  windows, 
and  no  windows  with  insulation  on  the  outside  surfaces.  The  outside  air 
temperature  cycle  used  was  40  - 100  F and  the  inside  air  temperature  was 
77.6  °F.  Figure  32  shows  the  computed  heating  load  profiles  for  the 
above  cases.  The  peak  heating  loads  for  single  pane  windows  was  50?o 
higher  and  occurred  approximately  two  hours  earlier  than  the  case  without 
windows.  The  peak  heating  load  for  double  pane  windows  was  7%  lower  than 
single  pane  windows.  Some  validation  by  measurement  of  the  latter  can  be 
seen  by  comparing  the  peak  heating  loads  as  shown  in  figures  29  and  30, 
about  4%  difference. 
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Heating  Load  Predictions 


c . 


Steady-state  methods  are  usually  used  for  predicting  maximum  heating 
loads  from  which  the  size  of  heating  equipment  is  selected.  Sometimes 
this  process  results  in  oversizing  of  heating  equipment.  To  illustrate 
and  compare  the  steady-state  procedure  and  the  dynamic  procedure  as  given 
by  the  computer  programs  in  Appendices  C and  D,  Table  4 was  prepared. 

The  values  listed  in  the  column  under  Steady-State  Method  in  Table 
4 were  calculated  for  the  experimental  structure  as  used  in  tests  6 
through  10,  and  for  the  outside  air  temperature  cycles  used  during  the 
tests.  The  steady-state  maximum  heat  flow  rate  was  calculated  using 
the  following  formula: 


(T.  - T)+(T.  - T)rU  A + lo08V(T.  - T) 


where 


q = heating  load,  Btu  hr 

-1  -2 

U = coefficient  of  transmission  for  the  floor,  Btu  hr  ft  F 
F 


2 

= area  of  the  floor,  ft 


T^  = average  inside  air  temperature, 

T = average  ground  temperature,  F 
S 

T^  = outdoor  temperature,  F 

U = coefficient  of  transmission  for  the  nth  surface 
n 


A = area  of  the  nth  surface,  ft 


n 


V = air  infiltration  rate,  cfm 
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The  first  term  corresponds  to  the  heat  transferred  through  the  floor. 

The  second  term  is  for  heat  transferred  through  the  walls,  windows  and 
roof.  The  third  term  is  heat  transfer  due  to  air  infiltration.  When  the 
above  equation  was  used  to  predict  the  maximum  heating  load  the  minimum 
outdoor  temperature  was  used  for  T^.  When  the  above  equation  was  used 
to  calculate  the  daily  average  heating  load,  the  daily  mean  outdoor 
temperature  was  used  for  T^. 

The  peak  and  daily  average  heating  loads  as  calculated  using  the 
steady-state  and  Response  Factor  methods  are  presented  for  comparison 
with  measured  values  in  Table  4. 

The  maximum  heat  flow  rates  as  calculated  by  the  steady-state  method 
for  the  conditions  during  tests  6,  7,  8,  9,  and  10  were  31,  59,  65,  68 
and  30  percent,  respectively,  higher  than  the  measured  rates.  The  maxi- 
mum heat  flow  rates  as  predicted  by  the  Response  Factor  method  were  6 
percent  or  less  of  the  rates  measured  during  the  tests.  The  above  high 
percentages  indicate  that  when  steady-state  maximum  rates  are  used  to 
size  heating  equipment  without  taking  into  account  the  heat  capacity 
effects  of  the  building  considerable  oversizing  could  result. 

When  comparing  daily  average  heat  flow  rates  between  the  steady- 
state  method,  the  response  factor  method  and  measured  values.  Table  4 
shows  that  all  values  are  reasonably  close  to  each  other  for  a given 
test  number  (about  10%  or  less).  This  was  expected  because  a minimum 
quantity  of  heat  energy  is  necessary  to  maintain  the  indoor  air  tem- 
perature over  a period  of  24  hours. 
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Comparison  of  Maximum  and  Average  Heating  Loads 
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7.  Conclusions 


The  NBSLD  computer  program  was  experimentally  validated  for  pre- 
dicting the  daily  indoor  air  temperature  as  it  is  influenced  by  known 
outdoor  temperature  conditions  and  the  mass  and  thermal  resistance  of 
the  building.  Furthermore,  when  the  inside  air  temperature  was  thermo- 
stated,  this  program  predicted  the  peak  and  daily  average  heating  loads 
and  may  therefore  be  used  to  size  equipment  needed  to  condition  the  in- 
terior of  a building  and  to  predict  energy  requirements.  It  was  shown 
that  steady-state  methods  of  heating  load  calculations  could  result  in 
oversizing  heating  equipment  by  30%  or  more.  The  NBSLD  dynamic  method 
takes  into  account  heat  storage  effects  and  therefore  predicts  the  peak 
heating  load  more  realistically.  The  maximum  percent  difference  between 
the  computer  calculated  peak  heating  load  and  measured  values  was  six 
percent,  and  the  average  difference  was  3.5  percent  for  the  five  tests. 

The  combination  of  mass  in  the  walls  and  roof  facing  the  interior  with 
insulation  placed  on  the  outer  surfaces  of  the  building  was  very  effective 
in  reducing  and  controlling  the  variation  of  the  indoor  air  temperature. 
This  desired  effect  was  also  predicted  by  the  computer  program.  When  the 
inside  air  temperature  was  not  thermostated  and  the  building  floated  in 
response  the  outside  air  temperature  condition,  placing  insulation  on  the 
inside  building  surface  reduced  the  variation  of  the  inside  air  temperature 
from  10  1/2  F to  5 1/2  F.  Furthermore,  when  this  insulation  was  placed 
on  the  outside  surface  of  the  building,  the  peak  to  peak  variation  in  the 
inside  air  temperature  was  reduced  to  2 °F.  In  addition,  comparing  cases 
of  no  insulation,  insulation  inside,  and  insulation  outside,  the  tempera- 
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ture  distribution  from  floor  to  ceiling  on  walls  and  in  the  indoor  air 
was  a minimum  when  insulation  was  placed  on  the  outside  of  the  building. 

The  effect  of  an  internal  mass  on  the  thermal  behavior  of  the  ex- 
perimental structure  was  generally  small.  An  internal  mass  may  have  a 
greater  effect  in  a less  massive  building.  On  the  other  hand,  windows 
had  a significant  effect  on  the  computed  thermal  behavior  of  the  experi- 
mental structure.  For  instance,  the  peak  heating  load  for  the  experimental 
structure  with  windows  and  insulation  was  507o  higher  than  the  same  build- 
ing without  windows.  Use  of  storm  windows  reduced  the  peak  heating  load 
by  ^ percent . 


8.  Acknowledgments 

The  authors  wish  to  thank  H.  E.  Robinson  for  his  many  ideas,  con- 
ceptual approach  and  technical  suggestions  prior  to  and  during  the  ex- 
perimental phases  of  this  investigation.  Also  appreciation  is  expressed 
to  D.  R.  Showalter,  J.  D.  Allen,  J.  M.  Dungan  and  J.  W,  Grimes  for  their 
technical  assistance  in  the  installation  of  the  instrumentation  and  day 
to  day  operation  of  the  experiments.  The  very  helpful  suggestions  and 
comments  made  during  review  of  this  paper  by  P.  R.  Achenbach  are  gratefully 
acknowledged. 


35 


vD 

cn 


Figure  1 Experimental  structure  in  environmental  chamber 
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Figure  2 The  basic  experimental  structure 


Figure  3 The  experimental  structure  showing  windows,  internal  mass, 
heaters,  and  insulation  on  the  outside  of  the  building 
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Figure  4 


Floor,  footing  and  window  details 
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Figure  5 The  experimental  structure  showing  thermocouple  locations 
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Figure  6 


Plan  view  of  thermocouple  locations 
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Figure  / Outside  sol-air  temperature  cycle 


r 


ib) 


Figure  8 (a) 

(b) 


Fan  heater 

Fan  heater  with  top  cover  removed 
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Figure  9 Comparison  between  measured  and  calculated 
inside  air  temperatures  for  floating  test  1 
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Figure  10  Comparison  between  measured  and  calculated 
inside  air  tempjeratures  for  floating  test  2 
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TEMPERRTURES 


NO  INSULATION,  WINDOWS,  NO  INTERNAL  MASS. 


Figure  11 


Comparison  between  measured  and  calculated 
inside  air  temperatures  for  floating  test  3 


TEMPERnrURES 


INSULATION  INSIDE,  WINDOWS,  INTERNAL  MASS. 


Figure  12 


Comparison  between  measured  and  calculated 
inside  air  temperatures  for  floating  test  4 
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INSULATION  OUTSIDE,  WINDOWS,  INTERNAL  MASS 


Comparison  between  measured  and  calculated 
inside  air  temperatures  for  floating  test  5 
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Figure  14  Comparison  of  inside  air  temperature  deviations  from  daily  average  inside 

air  temperature  for  identical  houses  for  cases  of  no  insulation,  insulation 
inside,  and  insulation  outside 
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Figure  15  Comparison  of  the  inside  air  temperature  deviations  from  daily 
average  for  identical  houses  with  internal  mass  (test  2)  and 
without  internal  mass  (test  1) 
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Figure  16  Comparison  of  the  inside  air  temperature  deviations  from  daily  average 
inside  air  temperature  for  identical  uninsulated  houses  with  windows 
(test  3)  and  without  windows  (test  1) 
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Figure  17  Variations  of  inside  wall  surface  temperatures  for  test  1 
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Figure  18 


Variations  of  outside  wall  surface  temperatures  for  test  1 


Figure  19  Comparison  between  the  measured  and  calculated  heat  fluxes 
at  the  inside  surface  of  the  floor  for  test  1 
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Figure  21  Comparison  between  the  measured  and  calculated  heat  fluxes 
at  the  inside  surface  of  the  wall  for  test  1 
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Figure  22  Computed  variations  of  the  heat  flow  rates  at  the  inside  surfaces  of  the 
roof,  walls,  floor,  and  the  window  for  the  case  of  no  insulation  (test  3) 
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Figure  23  Computed  variations  of  the  heat  flow  rates  at  the  inside  surfaces  of  the  roof,  walls, 
floor,  and  the  windows  for  the  case  of  insulation  placed  on  the  outside  surfaces  of 
the  building  (test  5) 
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Figure  24  Deviations  of  the  inside  air  temperature  from  instantaneous 
average  of  the  six  indoor  air  thermocouples  for  the  case  of 
no  insulation  (test  3) 


TEMPERATURE  DEVIATIONS 


□ 

3"  above  floor,  2 

' from  wall 

+ 

3'*  above  floor,  center  of  room 

A 

3"  below  ceiling. 

2 * from  wall 

o 

3"  below  ceiling. 

center  of  room 

o 

5*3"  above  floor. 

2 * from  wall 

X 

5*3"  above  floor. 

center  of  room 

-2- 


-3- 


-4- 


-ot 1 1 1 1 1 1 1 

0 3 6 9 12  15  18  21  24 

TIME  - HOURS 

Figure  25  Deviations  of  the  inside  air  temperature  from  instantaneous  average 

of  the  six  indoor  air  thermocouples  for  the  case  of  insulation  placed 
on  the  outside  surfaces  of  the  building  (test  5) 
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TEnPERflTURE  DEVIPTI0N5 
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Figure  26  Deviations  of  the  inside  air  temperature  from  daily  average  of 
the  six  indoor  air  thermocouples  for  the  case  of  insulation 
placed  on  the  inside  surfaces  of  the  building  (test  4) 
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Figure  27  Comparison  between  the  measured  and  calculated  heating  loads  for  test 
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Figure  28  Comparison  between  the  measured  and  calculated  heating  loads  for  t 
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Figure  30  Comparison  between  the  measured  and  calculated  heating  loads  for  test 
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Figure  32  Comparison  of  the  heating  load  profiles  for  the  cases  of  no  windows, 
single  pane  windows,  and  double  pane  windows  for  the  same  building 
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Appendix  A 


Air  Infiltration  Measurements  on  the  NBS 
Prototype  Building 


Thomas  K.  Faison 

Environmental  Engineering  Section 
Sensory  Environment  Branch 
Building  Research  Division 
Institute  for  Applied  Technology 
National  Bureau  of  Standards 


1.  Introduction 


In  the  early  stages  of  the  project  on  thermal  performance  of  the 
experimental  structure,  measurements  were  made  to  determine  the  magnitude 
of  air  exchange  between  the  structure  and  the  surrounding  chamber  during 
the  process  of  cyclic  temperature  changes.  Since  wind  forces  were  negli- 
gible during  the  testing  period,  the  major  driving  force  influencing  the 
exchange  of  air  was  the  thermal  difference  between  the  air  inside  of  the 
structure  and  that  of  the  surrounding  air  in  the  chamber. 

2.  Analysis  and  Instriomentation 

The  instrumentation  used  in  the  determination  of  the  air  exchange 
rates  was  developed  at  the  National  Bureau  of  Standards”^,  and  the  process 
of  measurement  was  that  of  the  tracer  gas  method  using  helium  as  the 
tracer  gas. 

The  rate  of  change  in  concentration  of  a tracer  gas  caused  by  ex- 
change or  infiltration  of  outside  air  under  a steady-state  temperature 
difference  is  expressed  by  the  formula: 

-V  (dc/dt)  = Kc  (1) 


— Coblentz,  C.  W.,  and  Achenbach,  P.  R.,  "Design  and  Performance  of  a 
Portable  Infiltration  Meter",  Transactions,  American  Society  of 
Heating  and  Air  Conditioning  Engineers,  Vol.  63,  1957. 


la 


where  V - volume  of  enclosure 


c = concentration  of  tracer  gas  at  time  t 


K = average  volume  of  air  infiltration  per  unit  time  for  the 


time  interval 


t = time 


When  c = c^  at  time  = o,  the  solution  of  Equation  1 is  as  follows: 


c = c e 
o 


-Kt/V 


(2) 


or 


Kt/V  = log  (c  /c) 
e o 


(3) 


Equation  3 shows  that  the  number  of  air  changes  occurring  during  time  t 
is  equal  to  the  natural  logarithm  of  the  ratio  of  the  tracer  gas  concen 
trations  at  the  beginning  and  at  the  end  of  the  time  interval. 


Prior  to  the  test,  the  apparatus  was  calibrated  and  brought  into 
equilibrium  with  its  surroundings,  then  helium,  the  tracer  gas,  was  re- 
leased into  the  room.  As  the  helium  was  introduced  it  was  mixed  with 
the  room  air  by  means  of  a portable  fan  and  the  final  mixture  of  air  and 
helium  contained  about  l/27o  of  helium  by  volume. 

Four  helium  sensing  elements  were  distributed  within  the  space.  Each 
sensor  was  positioned  3 feet  above  floor  level  and  4 feet  from  an  outside 
wall  near  each  of  the  four  corners.  Air  temperature  measurements  of  the 
two  spaces  were  recorded  during  the  test. 


3.  Procedure  and  Results 


2a 


Initially  a test  was  made  to  determine  the  amount  of  air  exchange 
through  the  structure  with  the  surrounding  environmental  chamber  prior  to 
cutting  openings  for  the  glass  windows.  Later  additional  tests  were  made 
to  determine  the  rate  of  air  exchange  when  glass  windows  were  introduced 
into  the  structure.  The  windows  were  of  a fixed  type  and  were  caulked 
in  place.  The  door  was  closed  for  all  tests. 

Measurements  were  made  at  the  time  of  day  when  the  air  in  the  en- 
vironmental chamber  was  lowest  and  unchanging,  providing  a maximum  tem- 
perature difference  and  air  exchange  between  the  inside  and  outside. 
Measurements  of  air  exchange  were  made  when  the  tightly  fitting  weather- 
stripped  door  was  normally  closed  and  when  all  cracks  around  the  door 
were  taped. 

For  the  building  without  windows  the  measured  values  of  air  exchange 
were  0.03  and  0.06  air  changes  per  hour  for  the  conditions  of  the  taped 
and  untaped  door,  respectively.  These  air  exchange  rates  for  the  basic 
structure  are  very  small.  In  fact,  they  are  the  smallest  ever  measured 
at  NBS.  They  do  provide  a minimum  value  for  comparison  with  other  tests 
and  show  that  heat  gain  or  loss  to  the  structure  was  almost  solely  by 
heat  conduction  and  the  influence  of  air  leakage  for  the  test  without 
windows  was  practically  negligible. 
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After  the  windows  were  installed,  single  glass  only,  additional 
measurements  were  made  to  determine  the  exchange  rate  under  these  con- 
ditions. The  same  procedure  was  followed  and  approximately  the  same 
temperature  difference  was  observed.  Under  these  conditions,  but  with 
the  windows  installed,  the  door  not  taped  and  no  insulation  on  the  walls, 
the  measured  value  was  0.38  air  changes  per  hour,  a significant  increase 
over  the  first  tests  having  no  window  openings. 
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Appendix  B 


Noise  Transmission  Measurements  of  the  NBS 
Prototype  Building 


John  S . Forrer 

Applied  Acoustics  and  Illumination  Section 
Sensory  Environment  Branch 
Building  Research  Division 
Institute  for  Applied  Technology 
National  Bureau  of  Standards 
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1.  Objectives  of  Tests 


Measurements  were  made  of  the  attenuation  of  outdoor  noise  provided 
by  the  prototype  concrete  block  structure  constructed  in  the  NBS  high- 
bay  environmental  laboratory  in  order  to  establish  the  feasibility  of 
noise  reduction  testing  in  such  a space  and  to  determine  the  sound 
transmission  characteristics  corresponding  to  four  different  conditions 
of  the  structure. 


2.  Building  Variations  Tested 

The  building  construction  during  the  first  series  of  tests  was  a 
simple  concrete  block  cubicle  with  a 20*  x 20*  floor  plan  and  a 10  ft 
high  ceiling  (outside  dimensions).  The  walls  were  made  of  8*'  x 8"  x 16" 
solid  concrete  blocks.  A concrete  slab  floor  and  a flat  four-inch  thick 
pre-cast  concrete  slab  roof  completed  the  enclosure.  A two-inch  thick 
solid  wooden  door  (foam  rubber  gasketed)  provided  the  only  break  in  the 
otherwise  solid  shell  of  the  structure. 

The  test  structure  configurations  employed  during  the  noise  trans- 
mission tests  were  as  follows  : 

1.  Concrete  shell  with  a single  wooden  door  (described  above). 

2 c Seven  32*'  x 40"  x 3/32**  single-pane  windows  installed  as 
shown  in  Figure  1 (bottom  of  sills  40  in.  above  floor). 

3.  Two-inch  thick  rigid  polystyrene  thermal  insulation  applied 
to  the  inside  walls  and  ceiling. 

4.  Insulation  removed  from  inside  the  structure  and  similar 
material  used  to  cover  the  outside  walls  and  roof. 
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3.  Test  Procedures 


Figure  B-1  shows  the  location  of  each  of  the  five  microphones  of  the 
receiving  room  array  (inside  the  house)  and  the  six  microphones  of  the 
source  room  array  (outside  the  house).  The  microphone  systems  employed 
one-inch  pressure-type  condenser  microphone  cartridges  with  attached 
preamplifiers.  Each  array  was  powered  by  a six-channel  microphone 
energizer  and  multiplexer  which  scanned  the  microphone  array  at  a rate 
of  five  channels  per  second.  The  multiplexer  output  was  fed  into  a 
one-third  octave  band-pass  filter  set.  The  filtered  signal  was  measured 
by  means  of  a precision  sound  level  meter  or  a graphic  level  recorder 
(see  Table  B-1). 

Calibration  of  the  measurement  system  was  performed  using  a calibrated 
pistonphone — a precision  sound  source  which  produces  a sound  pressure 
level  of  124  + .2  dB  at  a frequency  of  250  Hz  at  the  microphone  diaphragm. 

The  signal  for  the  noise  transmission  tests  was  provided  by  four 
speakers  energized  with  pink  random  noise*.  These  speakers  were  located 
opposite  the  outside  corners  of  the  house  as  shown  in  Figure  B-1.  The 
noise  reduction  provided  by  the  house  at  each  test  frequency  was  determined 
by  subtracting  the  one-third  octave  band  sound  pressure  level  measured 
in  the  receiving  room  from  the  corresponding  level  measured  in  the  source 
room. 


Pink  random  noise  is  white  noise  passed  through  a network  which  weights 
at  -3  dB  per  octave. 
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4.  Results 


The  curves  plotted  in  Figure  B-2  present  the  measured  noise  reduc- 
tion provided  by  the  house  for  each  of  the  four  variations  in  construc- 
tion, As  shown,  the  use  of  windows  caused  an  average  loss  of  sound  iso- 
lation of  about  10  db  for  frequencies  above  200  Hz.  The  addition  of 
thermal  insulation  either  on  the  inside  or  the  outside  improved  the 
acoustic  performance  but  not  enough  to  overcome  the  loss  from  windows . 

Data  was  gathered  at  frequencies  below  500  hertz  but  the  short  in- 
tegration times  used  in  the  r.m.s.  detection  system,  along  with  diffi- 
culties encountered  in  achieving  a uniform  sound  field  in  the  test  space 
rendered  the  measurements  inconclusive  for  frequencies  below  500  hertz. 
Specifically,  measurements  of  the  sound  distribution  inside  and  around 
the  house  with  the  speakers  energized  revealed  differences  in  the  range 
of  4-12  db  for  frequencies  below  200  Hz  in  the  sound  pressure  levels 
measured  at  microphones  in  the  same  array  in  the  receiving  room  and  for 
frequencies  below  500  Hz  in  the  source  room.  Differences  of  this  magni- 
tude render  a spatial  average  achieved  by  a five  or  six  microphone  array 
of  little  value. 
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Table  B-1  Instrumentation  for  Noise  Reduction  Measurements* 

1.  Bruel  and  Kjaer  Model  4220  Pistonphone 

2.  Bruel  and  Kjaer  Model  4132  Pressure  Microphone 

3.  Bruel  and  Kjaer  Model  2619  FET  Preamplifiers 

4.  Bruel  and  Kjaer  Model  221  Microphone  Energizer  and  Multiplexer 

5.  Bruel  and  Kjaer  Model  1612  Band-pass  Filter  Set 

6.  Bruel  and  Kjaer  Model  2204  Precision  Sound  Level  Meter  (used  during 

design  stages  1,  2,  and  3). 

7.  Bruel  and  Kjaer  Model  2305  Graphic  Level  Recorder  (used  during  design 
stages  3 and  4). 

8.  Kudelski  (Nagra  III)  tape  recording  of  pink  noise  used  as  signal 
source  in  design  stages  1,  2,  and  3. 

9.  Bruel  and  Kjaer  Model  1024  Sine  Random  Generator  used  as  pink  noise 
source  in  design  stage  4. 


Commercial  instruments  are  identified  in  this  report  in  order  to 
adequately  specify  the  experimental  procedure.  In  no  case  does 
such  identification  imply  recommendations  or  endorsement  by  the 
National  Bureau  of  Standards,  nor  does  it  imply  that  the  equipment 
identified  is  necessarily  the  best  available  for  the  purpose. 
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Figure  B-1  Floor  plan  of  NBS  High-bay  Environmental  Laboratory  and 

prototype  concrete  block  building.  Microphone  and  speaker 
positions  used  for  the  noise  reduction  tests  are  indicated. 
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Figure  B-2  Noise  reduction  versus  frequency  for  various 
construction  modifications  of  the  concrete 
block  building  in  the  High-bay  Environmental 
Laboratory . 
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Appendix  C 


Computer  Programs  (NBSLD)  to  Obtain  Heating  and  Cooling  Loads  and 
to  Estimate  Room  Air  Temperature  Change  Using  Thermal  Response  Factors 


Tamami  Kusuda 

Environmental  Engineering  Section 
Sensory  Environment  Branch 
Building  Research  Division 
Institute  for  Applied  Technology 
National  Bureau  of  Standards 


Computer  Programs  (NBSLD)  to  Obtain  Heating  and  Cooling  Loads  and 
to  Estimate  Room  Air  Temperature  Change  Using  Thermal  Response  Factors 

1.  Introduction 

The  NBS  computer  programs  called  NBSLD  are  a group  of  routines  to 
permit  the  determination  of  heating  and  cooling  loads  of  a room  based 
upon  a calculation  methodology  proposed  by  the  American  Society  of 
Heating,  Refrigerating  and  Air-Conditioning  Engineers  (ASHRAE)  Task 
Group  on  Energy  Requirements. 

For  a given  24-hour  weather  pattern  the  program  calculates  heat 
exchange  due  to  solar  and  sky  thermal  radiation  through  windows,  heat 
conduction  through  walls  and  roofs,  heat  convection  due  to  air  infil- 
tration and  internal  heat  generation.  Heat  exchange  is  computed  for 
every  hour  and  later  converted  into  the  room  heating  or  cooling  load 
in  conjunction  with  weighting  factors.  Details  of  these  calculation 
procedures  and  the  theoretical  background  for  the  weighting  factors' 
application  are  given  hereo  They  are  available  in  the  1971  ASHRAE 
publication  entitled  "Procedures  for  Determining  Heating  and  Cooling 
Loads  for  Computerized  Calculation  of  Energy  Requirements".  This  pub- 
lication was  prepared  by  the  ASHRAE  Task  Group  on  Energy  Requirements 
with  the  assistance  of  the  National  Bureau  of  Standards  and  the  National 
Research  Council  of  Canada. 
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The  ASHRAE  Task  Group  procedure  incorporates  what  is  considered  to 
be  the  most  up-to-date  computation  methodology  for  evaluating  the  dynamic 
aspects  of  building  heat  conduction  by  the  response  factor  method.  Since 
the  algorithms  employed  in  this  procedure  are  new  and  rather  complex, 
their  use  has  been  limited. 

Presented  in  this  report  is  the  Fortran  listing  of  the  NBS  program 
of  the  ASHRAE  Task  Group  algorithms  to  illustrate  the  use  of  this  modern 
and  powerful  technique  on  small  computers. 

All  of  the  routines  are,  therefore,  written  in  a close  accordance 
with  the  ASHRAE  Task  Group  algorithms  and  made  into  many  subroutines, 
each  of  which  could  be  used  independently  for  other  programs. 

Attached  are  the  Fortran  listings  of  NBSLD.  The  program  in  the 
form  of  punched  cards  or  on  magnetic  tape  is  available  from  the  Environ- 
mental Engineering  Section  of  NBS  including  assistance  for  its  use,  if 
desired.  Figure  1 shows  the  logic  network  for  NBSLD. 

1.  ABCD2,  ABCDP2,  DERVT,  GPF,  MULT,  RESF,  RESFX,  RESPTK: 

These  routines  are  parts  of  response  factor  calculation 
package  and  are  needed  for  the  accurate  evaluation  of 
thermal  time  lag,  damping,  heat  storage  in  exterior 
facing  surfaces  as  well  as  the  internal  furnishings. 

2.  DPF : Calculates  dew  point  temperature  when  the  partial 

vapor  pressure  is  known. 

3.  GLASS:  Calculates  solar  heat  gain  through  glass  when 

given  the  shading  coefficient,  orientation  type  of 
glass,  type  of  fenestration. 
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4.  OUTSID:  This  routine  calculates  the  outside  surface  tem- 

perature and  wall  heat  gain  by  taking  into  account  solar 
heating,  back  radiation  to  the  sky,  convective  heat  loss 
to  the  ambient  air  and  transient  heat  conduction. 

5.  PSYl : This  is  a simplified  psychrometric  routine  that  de- 

termines the  thermodynamic  properties  of  moist  air  when 
given  dry-bulb  temperature,  wet-bulb  temperature  anci 
barometric  pressure. 

6.  PSY2 : This  is  the  same  as  PSYl  except  that  the  dew  point 

temperature  is  used  instead  of  the  wet-bulb  temperature. 

7o  PVSF : This  routine  determines  the  saturated  vapor  pressure 

as  a function  of  temperature. 

8.  SHG:  This  is  the  ASHRAE  routine  for  calculating  solar  heat 

gain  through  glass. 

9.  SUN:  Calculates  basic  sun  data  such  as  angles,  cloud  cover, 

direct  and  diffuse  radiation  needed  for  solar  heat  gain  and 
solar  heating  of  the  building  exterior  surfaces. 

10.  TAR:  Calculates  transmission  and  absorption  characteristics 

of  glass. 

11.  WBF : Approximates  the  wet-bulb  temperature  when  provided 

with  the  enthalpy  of  moist  air  and  the  barometric  pressure. 

12.  WF : Determines  the  cooling  load  by  multiplying  the  heat 

gain  by  the  ASHRAE  weighting  factors.  (This  routine  was  not 
used  in  the  version  listed  in  this  report  because  it  incor- 
porates the  basic  calculation  used  for  deriving  the  weight- 
ing factors . ) 
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13.  RMTMP : Determines  the  room  temperature  as  a balance  of 

heat  gains  and  cooling  capacity  of  an  air  conditioning 
unit.  Since  this  routine  is  not  available  in  ASHRAE 
Task  Group  Algorithms,  detail  is  given  in  the  following 
pages  o 

14.  SOLVP : Solves  simultaneous  linear  algebraic  equations 

needed  in  RMTMP. 

15.  WEATHE,  WD,  DEC0DE:  This  package  is  a weather  decoding 

program  and  was  not  included  in  this  version  because  the 
weather  input  to  this  version  is  implicitly  defined  in 
the  following  section  on  input  data. 

16.  CCM:  This  routine  modifies  the  solar  radiation  for  a 

cloudless  sky  by  the  instantaneous  cloud  cover.  (This 
routine  is  not  included  in  this  version.) 

17.  FO:  This  routine  calculates  the  outside  surface  heat 

transfer  coefficients  from  the  weather  data.  (This  routine 
is  not  included  in  this  version  where  the  coefficients 

are  considered  to  be  input  data.) 
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' gure  1 Logic  network  for  NBSL^ 


2 . RMTMP 


Room  Temperature  Calculation  Routine 


Input:  NS  = number  of  heat  transfer  surfaces  in  the  room 

2 

(S(I),  1=1,  NS)  = area  of  the  heat  transfer  surfaces,  ft 
(M(I),  I = 1,  NS)  = number  of  response  factor  terms  for  each 

heat  transfer  surface 

(IX(I),  I = 1,  NS)  = index  for  the  thermal  storage  effect  for 

each  heat  transfer  surface 

IX (J)  = 1 for  thermal  storage  surface 
: IX(J)  = 0 for  non-thermal  storage  sur- 

faces such  as  windows  and  door 
(CR(I),  1=1,  NS)  = common  ratio  for  the  thermal  response 

factor  of  each  heat  transfer  surface 


M(I)  = 1,  CR(I)  = 0 if  IX(J)  = 0 
((X(I,  J),  Y(I,  J)  for  1=1,  NS),  J = 1,  M(D)  thermal  re- 
sponse factors  for  each 
surface 

X(I,  1),  Y(I,  1)  = overall  thermal  conductance  of  the  non-thermal 

storage  surface  and  all  the  other  response 
factor  terms  should  be  treated  as  zero  if 


IX(I)  = 0. 


Note:  For  the  calculation  of  X(I,J),  Y(I,J),  the  surface  heat  transfer 

coefficients  (both  inside  and  outside)  are  not  included. 
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((T0(I,t-J)  for  I = 1,NS),  J = = outside  surface  tem- 

perature history,  °F 

((TI(I,t-J)  for  I = 1,NS),  J = = inside  surface  tem- 

perature history,  °F 

TA  = air  temperature  of  the  room 

(H(I),  I = 1,NS)  = convection  coefficient  of  the  interior 

surface,  °F 

(F(I,K),  I = 1,NS),  K = 1,NS)  = radiant  heat  exchange  factors 

between  surfaces  I and  K, 
where  F(I,K)  = 0 if  I = 

K 

(R(I,t),  I = 1,NS)  = heat  input  per  unit  indoor  surface  at 

time  t to  the  surface,  such  as  solar 
heat  or  radiation  heat  from  the  lighting, 
equipment  and  occupants  to  the  surface 
(E(I),  I = 1,NS)  = emissivity  of  the  surface 

Q(I,t-l)  = heat  flow  at  the  Ith  surface  at  the  previous  time 

2 

period  or  time  = (t-l)A,  Btu/hr,  ft 
A = time  increment 

t = time  index  for  the  elasped  time  tA  hours 
CFML  = outdoor  air  leakage , CFM 

CFMV  = ventilation  air  rate,  CFM/°F  (at  time  tA) 

DB(t)  = outdoor  air  temperature,  °F 

TV(t)  = ventilation  air  temperature,  °F  (at  time  tA) 
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QBQUP: 


convective  component  of  internal  heat  from 
equipment,  Btu/hr 
Q0CPS:  convective  component  of  internal  sensible 

heat  from  occupants,  Btu/hr 
QLITE:  convective  component  of  heat  from  lights 

suspended  in  air,  Btu/hr 

lo  Basic  heat  balance  equation  at  the  I surface  (at  time  tA) 

M(I) 

Q(I,t)  = S {X(I,J)  * TI(I,t-J+l)  - Y(I,J)  * T0(I,t-J+l)} 

J=1 

+ CR(I)  * Q(I,t-l) 

' Ns 

= H(I)  * (TA(t)  - TI(l,t))  + S G(I,K)  * (TI(K,t)- 

K=1 

- TI(I,t))  + R(I,t) 

where  G(I,K)  = 4 * E(I)  * F(I,K)  * (TA  + 460)^  * 0.1714E-8 

2.  Total  heat  balance  for  the  room  air 
Ns 

Z S(I)  * (TI(I,t)  - TA(t))  + 1.08  * CFM  * (DB(t 
1=1 

* (DB(t)  - TA(t)) 

+ 1.08  * CFMV  * (TV(t)  - TA(t)) 

+ QEQUP  + QOCPS  + QLITE  = 0 
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3.  Letting  matrix  elements 


Ns 

A(I,I)  = X(I,1)  + H(I)  + Z G(I,K) 

K=1 

A(I,K)  = -G(I,K),  A(K,I)  = -G(K,I),  for  I = 1,  NS 

A(I,  Ns+1)  = -H(I) 

M(I)  M(I) 

B(I)  = -Z  X(I,J)  * TI(I,t-J)  + Z Y(I,J)  * T0(I,t-J) 

J=2  J=1 

- CR(I)  * Q(I,t-l)  + R(I,t) 

A(Ns+l,K)  = S(K)  * H(K)  for  K = 1,  Ns 

Ns 

A(Ns+l,  Ns+1)  = -1.08  * (CFML  + CFMV)  - Z H(K)  * S(K) 

K=1 

B(Ns+l)  = -QEQUP  - Q0CPS  - QLITE  - 1.08  * (CFML  * DB(t) 

+ CFMV  * TV(t)) 

TI(I,t)  and  TA  can  be  obtained  by  solving  the  following  Ns+1 
simultaneous  equations 


A(l,l),  A(l,2)  ... 

A(l,Ns+l) 

TI(l,t) 

B(l) 

A(2,l),  A(2,2)  ... 

• 

A (2, Ns+1) 

"k 

TI(2,t) 

• 

B(2) 

• 

• 

• 

A(Ns,l),  A(Ns,2)  . 

. A(Ns,Ns) 

• 

• 

TI(Ns,t) 

• 

• 

B(Ns) 

A (Ns+1),  A (Ns+1 ,2) 

...  A (Ns+1, Ns+1) 

TA 

B(Ns+l) 

— 
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3.  Input  Data  Needed  for  the  Fortran  Listing  of  NBSLD 


Input  data  needed  for  the  heating/cooling  load  calculation  are 
listed  on  the  following  pages  but  not  necessarily  in  the  card  reading 
sequence  of  the  Fortran  version  listed  in  this  report. 


Building  Number  (BLDGNO) 

Ceiling  Height  (HT ) 

Floor  Area  (AG) 

Number  of  Floors  (N0FLR) 

Number  of  Occupants  (QCU) 

Winter  Window  Overall  Heat  Transfer  Coefficient  (UGW) 
Ground  Floor  Heat  Transfer  Coefficient  (UG) 

Air  Change  Per  Hour  (AIRCHG) 

Latitude  (LAT) 

Longitude  (L0NG) 

Time  Zone  Number  (TZN) 

Month  (M0NTH) 

Day  (DAY) 

Elapsed  Hour  Since  Midnight  of  January  1st  (ELAPS) 

Electric  Power  to  the  Light  Watts  Per  Square 
Foot  of  Floor  (QLITY) 

Electrical  Power  to  Equipment,  Watts  Per  Square 
Foot  of  Floor  Area  (QEQPX) 

Ventilation  Air  Rate  (CFMV) 

Air  Leakage  Rate  (CFML) 

Maximum  Temperature  of  the  Design  Day  (DBMAX) 

Daily  Temperature  Range  of  the  Design  Day  (RANGE) 
Design  Indoor  Temperature  Condition  (DBIN) 

Design  Outdoor  Wet-Bulb  Temperature  (WBMAX) 

Design  Indoor  Wet-Bulb  Temperature  (WBID) 

Design  Winter  Outdoor  Temperature  (DBMWT) 

Design  Summer  Ground  Temperature  (TG) 

Design  Winter  Ground  Temperature  (TGW) 

Total  Number  of  Exterior  Surfaces  to  be  Considered 
for  the  Heat  Gain  Calculation  (NEXP) 

Index  for  the  Room  Temperature  Calculation 
Index  for  the  Standard  ASHRAE  Task  Group  Calcula- 
tion in  the  Special  and  Detailed  NBS  Calculation 
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Repeat  the  following  cards  for  NEXP  times 
Type  of  Heat  Transfer  Exposures  (ITYPE) 

1.  Roofs 

2.  Walls 

3.  Windows 

4.  Doors 

5.  Floors 

Type  of  Response  Factors  to  be  Used  (IRF) 

lo  Heavy  roof  construction 

2 . Light  weight  roof 

3.  Heavy  weight  exterior  walls 

4.  Light  weight  exterior  walls 
5o  Heavy  ceiling/floor 

6 . Light  ceil ing /floor 

7.  Heavy  partition  wall 

8.  Light  partition  wall 

U Value  of  the  Exposures  (U) 

Area  of  the  Exposures  (A) 

Orientation  of  the  Exposures  (AZW) 

0.  South  facing 
90.  West  facing 
180.  North  facing 
-90.  East  facing 

Radiant  Heat  Exchange  Factors  Among  Exposure  Surfaces 


If  the  construction  of  roof,  wall  and  floor  is  non-standard,  the 
following  information  is  needed  in  addition  to  the  standard  data  indi- 
cated above. 
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1 

2 

3 

4 

5 

6 

7 

8 

9 

10 

11 

12 

13 

14 

15 

16 

17 


Roof,  Wall.  Floor  Data 

Time  increment  of  the  temperature  data 
Number  of  roof  layers  (NR) 

Thermal  resistance  of  the  roof  inside  surface 

Jl,  fi,  p,  c , and  resistance  of  the  1st  layer  counted 

from  inside  surface  ...  (NR-2 ) Cards 

Thermal  resistance  outside  surface  of  the  roof 

Description  of  the  1st  layer  of  the  roof 

Description  of  the  2nd  layer  of  the  roof 

Description  of  the  NRth  layer  of  the  roof 
Number  of  wall  layers  (NW) 

Thermal  resistance  of  the  inside  surface 

I,  fi,  p,  c and  resistance  of  the  1st  layer  counted 

from  inside  (NW-2)  Cards 

Thermal  resistance  of  the  outside  surface  layer 
Description  of  the  1st  layer  of  the  wall 

Description  of  the  2nd  layer  of  the  wall 

Description  of  the  NWth  layer  of  the  wall 
Number  of  layer  of  the  floor  and  the  semi- infinite 
layer  (NF)  index  (if  basement  floor) 

Thermal  resistance  of  the  inside  surface  f?,  p,  c 
and  Res  of  the  1st  layer  of  the  floor  counted  from 
the  inside  surface 
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18  ly  pi  ^ 3.nd  Res  of  the  2nd  layer  of  the  floor 
(NF-1)  Cards 

19  ^.y  pi  ^rid  c of  the  earth  ...  if  basement  floor 

20  Description  of  the  1st  layer 

21  Description  of  the  2nd  layer 

22  Description  of  the  NFth  layer 


RASE  HtATiNG/COOLiNG  LOAr  CALCULATION  PROGRAM 

INCLUDING  The  room  temperature  change  predictions 
1= EXPOSURE  NUMBER,  1 = 1, 2, -NEXP 
I r YPE  ( 1 ) , E xposuKE  Type  ni.imber 

1 ROOF 

2 EXPOSED  'JVALLS 

3 vM  N D 0 A*  S 

4 DOORS 

5 GROUND  heat  transfer  SURFACES 

6 furnishings, partitio t i wALL^j, Party  a- alls  and  floor/ceilings 
7 open  surface 

p exposed  floors 

PR  barometric  pressure  in  of-  HG 

ITEMP=  temperature  RISE  INDEX 

IHT(I)  HEaT  transfer  index 

AvEHTG - -average  liFAT  GAIi.  FQR  ^ITE 

TSITHT--TOTAL  SITE  hEAT  C.AIN  EOr  2H  HOURS 

IHT  = -1  glass  S»JREACE(TRaNsPAREnT) 

iMT  = n opaouE 

I HT= 1 others  I SE 

Q1--HEAT  FLO*  through  EACH  EXPOSURE 
QSU^i--5E^lSIBLE  HEAT  GAIN 
QTLAT--LATEnT  HEaT  GAIN 
TOTHT--TOTAL  HEAT  GaIN 
QC--SENSlBLECOOLINGLOAn 
siteos--enTire  site  sensible  hEat  gain 

SITEQL--ENTIRE  SITE  LATE^.T  hEAT  GAIN 
51TETH--ENTIRE  SITE  TOTA(  HEAT  GAIN 
tU.OMAX--aUlLDlNG  MAX  HEAT  GAIN 
QSUMT--AVEKAGF  HEAT  GAIN 
SITELO--ENTIRE  SITE  COOLirJC,  LOAp 
SITMAX--5ITE  MAX  HEat  GAIN 
AyE5IT--SlTE  AVERAGE  HEAT  GAIN 
sq*int--site  heat  loss 

IREU)  RESPONSE  FACTOR  NIjMRER  APPLICABLE  TO  THE  SURFACE 
ABSP(I)  SURFACE  solar  HEaT  COEf-plCENT 
SHADE(I)  SHADING  COEFFICIENTS 
U(i)  exposure  u value 

uT(i)--u  value  without  External  surface  resistance 
H(I)  exposure  exterior  surface  thermal  conductance 
All)  exposure  area 

v'.^AZ(I)  A/ all  azimuth  angle  measured  clockwise  from  south 
TG --ground  temperature  for  cooling  load  calculation 
TV  = ventilation  air  temperature 
UG- -GROUND  HEa.T  TRANSFER  COEFFICIENT 

a G - -GROUND  heat  TRANSFER  SlJRFACE(aO  WHEN  NO  GROUND  FLOOR) 

T G A - - W I N T E R GROUND  TEMP 
OHMW7-  - winter  OUTDOOR  TE^^P 

laTslatitijoe  degree 
1.  ONGsLONGITuDE  degree 

TZN--TIME  Zone  number 

MnNTH--MONTH  OE  YEAR 
0 A Y - - D A Y 

QLITX--MAXImUM  lighting  load  in  WATT/FT2 


^'KQPX--MAX^fOUlPLO^DIN^ATT/f^T2  bbU 

rjEXP--NUMBE.R  OF  EXTEPIOK  HfAT  TRANSFER  SURFACES  b60 

fU.  DGNO--BUILDINGNUMBER  b70 

MT--HUlLDINGORnAELLlNGUNlTHElGHT  bBU 

(>PSX--MAx  OCCUPANT  SENSn^LF  LOAq  BTU /HR, PERSON  b90 

'JPLX--MAA  occupant  latent  lqAO  BTU/HR,PERSO^J  6 0 0 

OP--OE^POlNTTEMp,F  610 

lJCU--^AX^)UMREROFOCCUPANTS  ^20 

FLAPS=UAYSELAPSEOSINCEJAN*l  630 

UGL  AS--A  I NTER  glass  heat  T R A N S E R COEFFICIENT  6M0 

Ml inner  surface  CONvECTIVE  heat  transfer  coefficient  6b 0 

hr  inner  surface  radiative  HEaT  transfer  coefficient  660 

g,gg  raoation  heat  exchange  Surfaces  shape  factors  6/o 

X,Y,Z  RESPONSEFACORS  680 

these  response  facors  Should  nut  inclde  outside  surface  690 

ThERMaLRE5ISTANCE  \HEN1TEMP,E(j,0  7 on 

they  SHOUO  not  Ir^CLuDE  BcTh  the  outside  and  inside  thermal  7 10 

resistances  vvhENITFMP.En.1  7?0 

CFML  A I R LE  AK  age  7 30 

CFMVVEfJllLATION  78  0 

R A fraction  OF  LIGHTirG  POWER  THAT  GOES  I TO  FLOOR  7S0 

DPMAX  DESlGNOuTDnORDRY-BULci  TEMPERATURE  760 

RANGE  lOAlLYRAfjGE  OF  THE  OUTDOOR  TEMPERATURE  770 

ARMAX  DESIGN  OUTDOOR  AFT-BULD  temperature  7b0 

A bid  desion  INDOOR  wet-Bulb  Temperature  ■ 790 

ORIN  DESIGNINDOOPDRY-BuLRTfMPRATURE  800 

ITK  Index  to  calculate  room  temperature  rise  mhEN  not  AIR  CnNr.  IT  810 

ITK=1'WHENN0TAIRc0NDITI0NED  620 

ITEMP  INOEXTOuSE  ASHRAE^EIgHTINGFACTUR  830 

IFITEMP=o  ASHARE  weighting  factor  880 

common  /CC/  X(lU,100),Y(i0,10n),Z(ia,100),lTYpE(lU),lHT(lO),lRF(l0  6bO 

l),ABSP(in),U(l0),H(i0),Ml(lQ),A(l0),UT(10),T0S(l0,‘T8),TlS(10,88),b  660 

2(lQ,10),TOY(8e),DB(2B),QLlTE{2M),QEQUP(2R),OOCPS(2‘i),QIIlO),CR{lO)  870 

3,NR(1Q),QgLaS(10,28),IThsT  88Q 

dimension  XX(I0O,10),YY{lO0,10),ZZ(ina,10),TNEW(2‘T),Tlx(2M),Tl(88)  890 

1 ,Q0CUP(28)  ,QTL(28)  ,XDUM(  ino)  ,YDUM(  100)  ,ZDUM(  100)  ,7DUM(  lOQ)  ,Q0(  10)  9(..n 

REAL  LG(ft),LX(8),LlS(8),0G(8),0x(8),0IS(8),QG7(8),UXZ{8),«ISZ(8),S  910 

1ITEQS(28).SiTEDL(28),SITeTh(28),SITELD(2R)iTIF(10)  920 

DIMFNSIO.N  0LlTX(28),QEQuv(2M),DDESlN(10,28),QPE0PL{28),DDES(in)  93  0 

DIMENSION  y5UN(l0»28),0SKY(l0,2H),SHADE(l0),AZW(lu)  980 

dimension  NAMEBD(6),VT(in),DR(l0),MR(l0)iQGX(l0)  950 

dimension  DBNBS(28)/.26,,20,.lS,.10,.0b,.0i,03i.li«l9,.30,.83,.s7»  960 

1.69,.eU,,9U,.96,.99,1.0,.97,.90,.7b,.b7».M3,.33/  970 

REALLAT,L0nG,M0nTH,N0Flp  980 

DIME  NSlOrjLTYPE(iO),GG(l  0,10) 

COMMON  /SOL/  LAT,L0?gG,TZN,WAZ,ATiCN,D5T,LPYR,S(3b)  lOOO 

REA0(b,90U)QLITX  10  10 

READ  (b,90U)  QEOUX  1020 

READ(b,900)QOCUP  lC30 

SIGMA=a,i7l8E-8  1080 

HR  = 8,*(S3b***3)*5lC)MA  lOSQ 

i3  079GlJKLMNsl,20  lOfoO 

REAn(b,9lU,END=flOO)NAMFBD  1070 

READ(b,BBU)lROT,ISKlP  1080 

IF(NAMEbD(i),EQ,»»)G0T0800  1090 

I F ( I SK I P , ne  • 0 ) GO  TO  30  1100 

DO  1 0 I a I , 1 0 1110 

00100=1,100  112 C 
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X ( I , J ) *G , 1130 

YII,J)«U,  IIMO 

in7(l,J)=n.  11  so 

n020J=l,2M  116U 

SITF.Q5(J)aU,  1170 

S I T L (0  L ( J ) =:  0 , 1 1 B 0 

S I T F T H ( J ) = U . 119  0 

S I T E L D { J ) = U . 1 2 n f j 

20  CONTINUE  1210 

S Q A 1 N T = 0 , 12  2 0 

CaLLKESFX  (X,Y,Z.XX,YY,7Z,MR»DK,VT,in)  1230 

' R I T F ( 6 . 8 2 0 ) . 12  8 0 

PH=29.92  12S0 

RF.  A[5  (S.vnO)  LAT, LONG,  TZ^, MONTH, DAY, ELAPS|UG,UGLAS  1260 

.V  R I T F ( 6 , fl  S 0 ) 12  7 0 

A rite  (6,flMn)  LAT,LnNG,T/N,MONTHiDAY,ELAPS,UG,UGLAS  12BU 

WFAD(b,90U)0LlTY,0EQPX,CFMV,CFML  1290 

A R I T F:  ( 6 , R 6 0 ) 13  0 0 

vVRITF(6,R80)QLITY,QEQPy,CFMV»CFML  l3lO 

READ  (b,900)  OBMaX,RANGE,DrIN,AbMAX,ABID,DBMAT,TG,TGA,tv  1320 

R I T F ( 6 , R 7 0 ) 1 3 3 0 

A RITE  (6,RM0)  DOMAX,RANGr,DBIN,wBMAX,ABID,DBMAT,TG,TGA,TV  1380 

CALL  PSyI  (DBMAX,'ABMAX,P!^,OP,PV,AOUT,HOUT,VnUT,KHOUT)  13S0 

CALL  PbYl  (UBlN,vVBlO,PB,rPlD,pV,MD,HIND,VlN,RHlN)  1360 

iA  V = A n U T . 13  7 0 

I N = V'/  I 0 1 3 P 0 

A = A0UT  1390 

T I 0 = P B I N 18  0 0 

3 0 READ  (b,900)  ROOMNO,HT,Af.  ,NOFLR,t>)CU,AlRCHG  1810 

v\RITF(6,830)  1820 

'RITE  (6,r8q)  ROOMNO,HT,aG,nOFUk,QCU,AIRCHG  1830 

RFAD(b,R90)MEXP,lTK,ITFMp,ITHsT  1880 

00110I=1»NEXP  18S0 

RFAD  (5,920)  ITYPE(1),IRF(I)|U(I),A(I),AZ'A(I),DUm,5HADE(I),AHsP(I)  I860 
REA0(b,9(j0)(G(I,J),Jsi,NFXP)  1870 

L T YPE  ( I ) = I T YPE ( I ) 18  80 

I F ( I T Y E ( I ) . E Q • 7 ) G 0 T 0 1 1 0 18  9 0 

K = I RE  ( I ) 1 boo 

IF{Y(K,i).GT.l.)lRF(I)=rIO,  1510 

NR  ( I ) =MR  ( K ) 1520 

UT ( I ) -V  T ( K ) 1530 

CR(I)=DR(<)  1580 

IF(NR(I).Gt.8r)NR(I)=8r  1550 

IF(1TYPE(I).EQ.3)AB5P(I)=0,  1560 

IF  ( 1 type ( I ) *EQt 5 ) A0SP(I)=O,  1570 

IF  ( I T YPE { 1 ) . EQ  . 6 ) ABSP(I)rO.  ISfcU 

I HT ( I ) = I 1590 

IF(1TYPe(I),EQ.3)IHT(I)=-1  1600 

H(I)=8.08  1610 

M 1 ( I ) =0 . 582  1620 

IE  ( I T Ype  ( I ) . EQ . 6 ) H(I)=n,  1630 

lE(lTYPE(l).EQ.5)Hl(I)st  0 ,162  168  0 

IF(U(I))80,80,50  1650 

MP  PUsl./OT(T)+l,/HI(l)  1660 

IF  ( I T YPE  ( I ) • NE  . 6 ) R U « R U + I , / H ( I ) 1670 

0(I)*1*/R0  1680 

SnCONTINOE  1690 

IF(X(K,2))110,60,110  1700 
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60  IF  70,80,70 

70  Ksl./U(I)-l,/H(I) 

GO  TO  VO 
BO  Rs 1 . /U ( I ) 

90  UT(I)*1./R 

IF  (ITFMP,Ne.O)  UT(I)«1*/(R-1*/(HI(I)+0*V)) 

I F ( UT ( I ) ) 1 00  » I 00  I I 1 0 

inn  UT(i)=ioo, 

iin  coruiNUE 

,V  R I T E (6,1170) 
no  120  I = 1 ,NEXP 
AZW(  I )=AZW(  I ) + IR0T 

IF  ( AZM  I ) *GT,  1 80.  ) AZW ( I ) =AZW ( I ) -360  . 

WRITE  (6.9J0)  I.ITYPEm.lHTm.IRFdl.ABSPlIl.Udl.Hdl.Am.AZWl 
1I),SHADE(1>.0T(1) 
ipn  CONTINUE 

^4EXPl=NExP-l 
NEXP2=NEXP-2 
DO  160  IsljNEXP 
G 5 U M a O . 

NE  X = N£  Xp  1 

IF  (I.EQ.NIEXP)  NfXsMEXPZ 
M E X s N E X ♦ I 
DO  130  J= 1 » NE  X 

IF  (I.CiT.J)  G(I,J)=G(JiI)*A(J)/A(I) 

130  G5UM=G5UM+G(IiJ) 

IF  (GSUM-l#)  1M0,1S0i150 

190  G(I»NEX)=i.-GSUM 

GO  TO  1 6U 
150  G(I,MEX)=0* 

160  continue 

ARITE(6,1180) 

^^.RITE  (6,ll90) 

DO  170  I=1»NEXP 

^MRITE  (6,1200)  I,(G(I,J),J=l,NExP) 

170  CONTINUE 

DO  IPO  1 s I , NE  XP 

00  IPQ  JsliNEXP 
180  G ( I , J ) *HR  *(3  ( I , J ) 

1 I sO 

DO  200  1 = 1 , NEXP 

IF  ( 1 T YPE  ( i ) .EQ  . 7 ) GO  TO  200 
11=11+1 

I TYPE ( I I ) a 1 T YPE ( I ) 

I RF ( I 1 ) = I PE ( I ) 

U ( I 1 ) =U ( I ) 

A ( I I ) »A  ( I ) 

AZW ( I I ) = AZW ( I ) 

SH ADE ( I I ) aSHADE ( I ) 

ABSP (11)  a A^SP ( I ) 

NR  { I I ) =NR  ( I ) 

UT ( I I ) *UT  ( I ) 

CR ( I I ) *CR  ( I ) 

IHT ( I I ) = I HT ( I ) 

H ( I I ) =H  ( I ) ■ 

H 1 ( I I ) «H  I ( I ) 

J J*0 

DO  190  Ja I ,NEXP 


17  10 
1720 
1 730 
1 790 
1 7bO 
1760 
1770 
1 7B0 
1 790 
1 800 

18  10 
1 820 
1 830 
I 8PQ 
1 8S0 
1860 
1870 
I 880 
1890 
1 900 

19  10 
1920 
1930 
1 9H0 
1950 
I960 
1970 

1 980 
1990 
2000 

20  1 0 
2020 
2030 
20P0 
2050 
2060 
2070 
2080 
2090 

2 100 
2 110 
2 I 20 
2 I 30 
2 IPO 
2 150 
2 I 60 
2 1 70 
2 180 
2 190 
2200 
22  10 
2220 
2230 
22P0 
2250 
2260 
2270 
2280 


i 
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IF  (LTYPE(J),EQ,7)  go  TO  l90 

2290 

J vJ  = J J + 1 

2300 

GG(II,JJ)=&{I,J) 

23  10 

1 9'^ 

C 0 N T I rj  U E 

2 320 

2nn 

C C ) (M  T I N U E 

2330 

NF  XP=  I I 

2 3 80 

A R I T E ( /,  , 1 2 1 0 ) 

23bu 

no  210  I=liMEXp 

2 3 60 

00  210  J=IiNEXP 

23  70 

2 1 r- 

G(IfJ)=GG{I,J) 

2 3 80 

R I T E ( 6 , 1 1 8 0 ) 

2390 

.'v  R I T E ( 6 , I 1 9 0 ) 

2 8 00 

0 (J  2 2 0 I = 1 , N E X P 

28  1 0 

2 ?n 

A K I T E (6,1200)  I , ( G ( I , J ) , J s 1 , N E X P ) 

2820 

,V  R I T E ( 6 , 1 1 7 n ) 

2 8 30 

f)0  230  1 = 1,  N EXP 

2880 

2 30 

A R I T E (6,930)  I , I T Y P E ( I ) , I H T ( I>  , I R F ( I ) , A B S P ( 1 ) , U ( I ) , H { I ) , A ( I 

) , A Z A ( 

2850 

1 I ) , SHADE ( I ) , UT  ( I ) 

2860 

II 

o 

• 

2 8 7 0 

«v  R I T F (6,980) 

28  80 

0 0 2 80  1 = 1,28 

2890 

2M0 

0B(  I ) = (DBMAx-RANGE)4-(RANGE*DBNbs(  I ) ) 

2 500 

5 U M = 0 * 

25  1 0 

DO  250  1=1,28 

2520 

250 

SUM  = 5UMi-dQ(  1 ) 

2530 

D f3  M = S U M / 2 8 • 

2580 

■A  R 1 T E ( 6 , 9 b 0 ) ( D B ( I ) , 1 = 1 , 2 8 ) , 0 B M 

2550 

DO  260  1=1,28 

2560 

?<sr 

T I X ( I ) =T  I 0 

25  70 

suM  = n . 

2580 

00  270  1=1,28 

2590 

2?n 

SUM  = 5UM-fT  I X ( I ) 

2 6 0 0 

TIM=SUM/28. 

26  1 0 

A RITE  (6,960)  (TIX(I),  1 = 1,28), TIM 

2620 

vRITE  (6,970)  QLITX 

2630 

A R I T E ( 6 , 9 B 0 ) Q E Q U X 

2680 

'^RITE  (6,990)  QOCUP 

2650 

CFM'A'T  = AG*HT/6n.*AlRcHG 

2660 

C F M L = 0 • S ♦ C F M W T 

26  70 

cem=cfml+cFmv 

26tj0 

0 L 1 T 0 = 0 L I T Y * A G • 3 , 8 1 3 * N 0 E L R 

2690 

QEQPO  = (JEijPX*AG*3.8l3*NOF|R 

2 700 

00  280  J=l,28 

27  1 0 

0LITE(J)=QLITX(J)*QLIT0 

2 7 20 

QEQUP(J)=QEqUX(J)*QEQPO 

2730 

QTL(J)*8880.*CFM*W0uT 

27HO 

2 80 

CONTINUE 

2 750 

DO  290  1=1,9 

2760 

U 0 ( I ) = 0 . 

2 7 70 

Q I ( I ) = 0 . 

2780 

2 90 

T N E //  ( I ) = 0 . 

2790 

C 

DRH  = TIM=  reference  TEMPEfx'ATURE 

2 800 

D B M = T I M 

28  10 

S ( 1 ) =L at 

2820 

S ( 2 ) = L 0 N G 

2830 

5 ( 3 ) =TZN 

2880 

S ( 8 ) =ELAPS 

2850 

S ( 6 ) * 1 . 

2860 

18c 


S ( 7 ) =0 . 2 

28  70 

5 ( 8 ) * 1 . 0 

2880 

S ( 3 3 ) = 1 . 

2890 

A' RITE.  (6,1220) 

2 9 0 0 

DO  3S0  I=l,NtXP 

29  1 U 

IF  (ITYPt(I),LT.s)  GO  TO  3lQ 

2920 

on  300  J=li29 

29  30 

(^5UN{  I ,J)=U, 

2 9 8 0 

0GLA5(  I ,J)=0. 

29b0 

3 0^^ 

QSKY( I iJ)=0. 

2960 

GO  TO  380 

2970 

3 1 C' 

A A Z = A Z W ( I ) 

29c0 

S ( 9 ) = A A Z 

29  90 

S ( 1 0 ) =9Q  , 

3 0 G 0 

IF  (ITYPF(I),Fq.i)  S(10):=0. 

30  1 0 

DO  330  J*l»29 

3020 

Q 5 K Y ( I 1 J ) = U , 

3030 

IF  (ITYPt{I),r0.1)  Q5KY(I,J)=2U. 

3080 

I I MF  = J 

3050 

S ( S ) = T I MF  • 

3060 

CALL  SUN 

3070 

IF  (5(25), GT.n.)  GO  TO  3?0 

3080 

0 S U N ( I » J ) = U , 

3090 

QGLAS(  1 ,J)=0. 

3 1 t;o 

GO  TO  330 

3 110 

3 7C 

QSUN ( I 1 J ) -S ( ?5  ) * aBSP ( I ) 

3 1 20 

Q G L A S { I , J ) = 0 

3 1 30 

IF  (IHT(i).GT.O)  go  to  330 

3 180 

CALL  glass  ( SH ADE  ( I ) , 1 . , 1 . , QGLAS ( I , j ) ) 

3 1 50 

3 3C 

CONTI NUE 

3 1 60 

3H0 

•".RITE  ( 6 , 1 0 0 0 ) I 

3 170 

A RITE  (6,1U10)  (QSUN(IiJ),Jx1,28) 

3 1 bO 

35H 

^RITE  (6,l0l0)  (QGLa5(I,J),J=1i2B) 

3 1 90 

DO  360  J=  1 *29 

3 2 00 

ri  (J)=TIX(28-J+1  )-TIM 

32  10 

Of)  360  I = 1|NEXP 

3220 

3A- 

TOSI  I ,J)=Db(28-J+I  ) -DRM 

3 2 30 

DO  370  J=?b,88 

3280 

TI  (J)=TI  (J-28) 

3250 

DO  370  IsliNEXP 

3260 

370 

T05 ( I , J ) =TOS ( I , J-28  ) 

3270 

IF  (ITEmP.NE.O)  go  to  39 

32  80 

DO  3 8 0 L = 1 » 8 

3290 

LG ( L ) =0 . 

3 3 0 0 

LX ( L ) -0 . 

33  10 

L I 5 ( L ) =0  , 

3320 

Q G { L > = 0 , 

3330 

Q X { L ) « 0 . 

3380 

380 

Q 1 S { L ) = 0 . 

3 350 

390 

CONTINUE 

3360 

D 08001  = 1, NEXP 

33  70 

f)0  800  J=l,88 

3 380 

800 

r I s ( I , J ) =0 . 

3390 

T A-T  I M 

3800 

DO  720  N=  1 , 7 

38  10 

IF(N.NE,7)G0T0Bi0 

3820 

Q S U M T = 0 . 

3830 

RLDM A X*0 , 

3880 
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M 1 n 

CONTINUE 

3450  1 

DO  77D  NK= 1 , 24 

3460  . 

DO  440  I = 1 , NEXP 

34  70 

i^O  420  NTT*2»4R 

34b0 

^ 7 n 

TOYirJTT  )=TUS(  I , NTT-1  ) 

3490 

DO  430  NTT=2,48 

3 5 0 0 , 

M 3H 

T 0 5 { I , N T T ) = T 0 Y ( N T T ) 

35  10 

M M Cj 

C 0 N T I fJ  U E 

3520 

IP  (ITLMp.NE.n)  GO  TO  49  f' 

3530 

r A*T 1 X ( NK  ) 

3540 

DO  450  L=2i8 

3550 

Q G Z ( L ) = Q G ( L - 1 ) 

3560 

Q X Z { I-  ) = Q X ( L - I ) 

3570 

GISZ(L)=QIS(L-1  ) 

35B0 

*4  b 0 

CONTI NUE 

3590 

on  460  Ls2»8 

3600  ; 

Q G { L ) = 0 G Z ( L ) 

3610. 

N X ( L ) = 0 X Z ( C ) 

3620  . 

QlS(L)=yiSZ(L) 

3630  : 

4 

c 0 N T UJ  u E 

3640  . 

on  470  NTT=2i4B 

3650  : 

4 7 0 

T0Y(NTT)=Ti  (NTT-1  ) 

3660 

DO  480  NT  T = 2 , 48 

3 6 7 0 ‘ 

4RP 

Tl(NTT)=roY(NTT) 

3680  . 

SUMQGsU , 

3690 

4 90 

CONTINUE 

3700  . 

DO  540  I=1,NEXP 

37  10 

< = I RF ( I ) 

3720 

0 0 5 0 0 J = 1 » 4 8 

3 7 30 

XDU^MJ)=X(^,J) 

3740  : 

YDUM{J)=Y{K,J) 

3750 

ZDUM(J)=Z(^,J) 

3760 

TDUM ( J ) = TOb { I , J ) 

37  70 

IF  (lTYPE(I).fcQ*6)  TDUM(j)=tIS(I,J) 

3 7 80 

IF  (ITYPE{I),EQ.5)  TDUM(j)sTG-TiM 

3790 

IF  (ITEMP.Nt.n)  T I ( J ) =T 1 s ( I , J ) 

3800 

IF  ( TDUM ( J ) , GT • 1 00  • , OK  , T I ( J ) . GT  , 1 no . ) GO  TO  760 

38  10 

IF  ( TOUM { J ) • LT . - 1 00 . .OR , T I ( J ) ,Lt • - 1 00 . ) GO  TO  76q 

3820 

b n i : 

CONTI NUE 

3830 

U X = U ( I ) 

3840 

IF  ( H ( I ) ) 520 , 520 ,510 

3850 

b 1 0 

RX= 1 . /U ( I ) - 1 , /H ( I ) 

3860 

LJ  X = 1 . / R X 

38  70  J 

b 70 

CONTI NUE 

3680 

CALL  0UT5IU  (XDUM,YDUM,ZrUM,CR(l),UX,H(I),DB(NK),TlM,QO(I),Ql,] 

) , 'V 

3 8 90 

iSUNI  I ,NK)  ,QSKY{  I ,NK)  ,T0U^-,TI  ,TNe^,(0,TA,ITLMP) 

3 9 0 0 . 

DO  530  Js 1 , 48 

39  1 0 : 

b3r^ 

T05 ( I , J ) =TDUM ( J ) 

3920 

TNEW( I IsTNLwO+TIm 

3930 

b40 

CONTINUE 

3940 

QOCPS(NK)=QOCUP(N<)*IG.*(100«-Ta)*QCU 

3950 

Q0CPLs10,*(TA-60,)*Q0CUP(NK)*QCu 

3960 

IF  (TA-lOO*)  560,550,550 

3 9 70 

bbO 

QOCPS{NK)sO, 

3 9 80  ] 

Q0CPL=40n.*Q0CUP(NK)*QCU 

3990 

GO  TO  580 

4000 

560 

I F ( T A-60  . ) 570,580 ,580 

4010. 

570 

Q0CP5(NK)a4Q0.*Q0CUP(NK).QcU 

4020 

2 0c 


I 


QOCPL*0,  *4030 

bPn  (jpeopl(nk)=qocpl  momu 

SUML  = QTL(N^)-48MQ,*CFM*;wIN  + Q0CPL  Musn 

■0TLAT*-5UML  M060 

c ,'SOM  INSTAmTANFOiiSHElATgAIN  ^07  [) 

C SUMQG  IMSlANTAMtOljSSOLARHEATGAlM  MO  80 

c 01  CONDUCTION  heat  transfer  MOVU 

wSUM=l.08*EFML*(TA-0B(NK))^.l,08*CFMV*(TA-TV)-QLITE(NK)-QEQUP(^'K)-''/  MlfiO 
I 0 C P S ( N K ) M 1 10 

0 0S9□I=;l,^JEXP  M12  0 

8 9 n 0 S U M = Q S U M + A ( I ) • { g I ( I ) - g G I A 5 ( 1 I N K ) ) M 1 3 0 

lF(^>.NE,8)GnT06S^  MIMU 

IF{NK.NE,1)G0T0800  MlbO 

\RlTF(6,lU2n)NAMER0  H160 

f ONT  I NUE  H 1 70 

.VRITF  (6,1030)  Nk,(TNEA(1),  1 = 1, 9), DB(NK)  M180 

IF  (ITEMP)  720,720,650  . M 1 90 

610  IF(MK.NE.1)G0T0620  M2  0 0 

VRITE(6,10M0)NAMERD  M2lD 

6? n continue  M220 

TOTHT=gL+QTLAT  M230 

rRITE  (6,1US0)  NK,(Ql(I),I=l,9),QSUM,QTLAT,gL,T0THT  M2M0 

D0630I=1,NEXP  M2S0 

QnE5TN(I,NK)sQl(i)*A(I)  M260 

6 3 n C 0 N T I N U E M 2 7 0 

S I TEDS ( Nk  ) =S I TEQ5 ( NK  ) ♦QS"M  M2B0 

SlTEQL(NK)=5lTEgL(N<)-*-gTL.  AT  M290 

S I TE T H ( NK  ) *S I TE T H ( NK  ) ♦ TOT Ht  M300 

S I TFLO ( NK ) =S I TELD ( Nk ) +QL  M31G 

lF(QL*GT.8LDMAX)G0T0^Mn  M320 

BLDMAX=QL  M33U 

T 0 THT  X = T0THT  M3mO 

I maX=NK  M3S0 

6M0  IF(N,EQ,7)g5UfMT  = QSUMT  + T0THT  M360 

G 0 T 0 7 2 0 M 3 7 0 

68nn06R0I=l,NEXP  M3fiO 

Df1  660NTT*2,M8  M390 

660  TOY(MTr)=TlS(I,NTT-l)  MMOO 

DO670NTT=2»9fi  MMlO 

6 70  T I S ( I , NT T ) =TOy ( NT T ) MM20 

6 H 0 C 0 N T I N U E M M 3 0 

TV  = D8(fMK)  MMMO 

CALL  RMTMP  (NEXP,NK,TV,CfNl,CFMViR,TIM,TA|TIF,QL,ITK)  MMSO 

IF  (TA.GT.TIN)  GO  TO  890  MM60 

DP  I rOP  ID  MM  70 

H 1 N = W B I D M M 8 0 

HINrHIND  MM90 

I N = V;  I D M 5 0 0 

G0T0700  M510 

690  CONTINUE  mS20 

Q0CPL=gOCPL/in60.  Mb30 

.«I  = (M.5*CFNL*^«/A  + M.b*CFMV*Wv-»-Q0CpL)/(M.5*CFML'^M.5.CFMV)  MbMO 

PvI=PB*AI/(Q.622+WI)  MbbO 

DPI=DPF(PvI)  Mb60 

CALL  PSy2  (TA,DPI,PB,WBI^!,pvI,WlN,HIN,VIN,RHIN)  *+570 

700  continue  ‘ MB80 

IF  (N.LT.6)  GO  TO  720  *4590 

IF(N,Eg.7)G0T06l0  *+6U0 
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I F ( NK  • nE  . 1 ) Go  TO  7 I 0 *4610 

rtKlTE(6,lU60)NAMEBD  ^620 

7 \ n A'RITF  (6,1U70)  NK,(TlF(I),lsl,V),TA,\ft8lN  "4630 

7 ?n  COMT  I MUE  M6M0 

(0SUMT=«SU'1T/2S,  '46bO 

.j/'lNT=1.08*CFMWT*(TlM-DB-v'vT)+UG*(TlM-TGV\/)*AG  B660 

Q v';  I N T = I . 0 8 * C F M M ♦ ( T I M - D b T ) *4  6 7 T) 

[)0  7‘4QI  = i,nEXP  B63U 

IF{IHT(l).LT,0)U(I)=UGlAS  H69Q 

iF(ITYPE(1).NE.S)GOTO730  *4700 

« 'T  HJ  T = « v-n  P4 1 + U G * A ( I)  # ( T I M ^ T G A/  ) M 7 1 U 

GO  TO  740  4720 

73 IF(ITYPE(1).EQ.6)G0T0740  H730 

IF(ITYPE(I),EQ.7)G0T0740  *4740 

0aINT  = QA/INT-»-A(I)*U(I)*(7iM-DBM1A/T)  4750 

74nC0NTINUE:  4760 

bQ^INTsSQMNT+QMNT  4770 

A RITE  (6.1140)  WWINT, TOT mTx.OSUhT  47fl0 

IF(lTK,fjE.0)GOT07V0  4790 

0 07501=1, N EXP  4800 

0GX(1)*QGLAS(I,ImAX)*A(I)  4810 

7 5 0 0 D E 5 ( I ) = (J  0 E S I M ( I , I M A X ) 4 8 2 0 

CFM  = CFML4-cEMV  4830 

CALL  output  (DBMAX,AjBMAx,DRlN,WBlD,WoUT,\MN,QGX,CFM,WLlTE(IPAx).QO  4840 
lCPS(IMAX),UPEnPL(IMAX),(vr'ES,AZW,ITYPE,NEXP,NAMEBD)  4850 

GOTO  790  4860 

76nu9ITF(6,l080)N  4870 

7RITE(6,l090)  4880 

007700*1,48  4890 

770  RlTE(6,lllO)J.(TOS(I,j),I  = l,lO)  49  (JO 

.vRITE(6,l0fl0)N  4910 

\RITE(6,llOO)  4920 

D07B0J*lt48  4930 

780  ARITE(6,1110)J,(TI5(I,,J),I  = 1,10)  4940 

79nC0NTINUE  4950 

800  CONT I MUE  4960 

A R I T E ( 6 , I 1 2 0 ) 4 9 7 0 

.V  R I T E ( 6 , 1 1 3 0 ) 4 9 8 0 

S I TMA  X = 0 . 4990 

T5AVE=0.  5000 

T5  I ThT  =0  , 50  10 

008101=1,29  5020 

SQLLD  = S I TEUL  I I ) AS  I TELD ( I ) 5O3O 

T S I ThT  = TS I 1 HT  + 5 I TETh { I ) 5040 

IF  ( 5 I TE TH ( I ) . LT  . S I TM AX  ) S I T M A X = S I T E T H ( I ) 5050 

IF(5QLlD,ET.TSAVE)TSAVF=SQLL0  5060 

vVRITF  (6,1150)  I,SITEQS(I),SITEqL(I),SITETH(I),SiTELD(I),SQlLd  5070 

810  CONTINUE  5080 

A V E H T G* T 5 I T H T / 2 4 • 5090 

A RITE  (6,ll60)  5qMnT,SItMaX,AVeHTG,TSAVE  5100 

.V  R I T E ( 6 , 1 1 2 0 ) 5 110 

STOP  5120 

C 5 130 

C 5 140 

82OF0RMAT(1h1)  5150 

830  FORMAT  (0Hl  0LDGNO,8X«HT»,8X»A6»,5X*NOFLR»,7X»QCu»,4X*aIRCH6*)  5160 

04OFORMAT(lOf'lO,l)  5170 

850  FORMAT  (7X*LAT*,6X»LONG*,7X»TZN»,5X*MOMTH*,7X*DAY»,5XtELAPS»,0X*UG  5180 
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]»,SX»UGLAS*)  bl90 


8^0 

FOKM  A T 

(SX’QLITY*  ,6X*QEypXf  ,6X»CFMV*  ,6X*CFML*  ) 

S2U0 

P 7 

r 0 W M A T 

(SX*DBmAX»5X»RANGf*  ,6X*UbIN*  ,5X*ABMAX*  ,6X»y\bID»  ,5X*DBMaT*  ,« 

S2  1 0 

1 X » T G • , 7 X » T G • , 8 ^ V ♦ ) 

S220 

9 « ■ 1 

format 

(1017) 

52  30 

fl  9 M 

format 

1 n I 7 ) 

5290 

V n 0 

FORMAT 

( 1 on . n ) 

525U 

9 1 n 

format 

( 6 A 6 ) 

5260 

920 

format 

( 2 I / , 6F  7 • D ) 

5270 

930 

r 0 R M A T 

{13»3llO,BFin,2) 

5280 

9 M Li 

F 0 R M A T 

(*l  ENVIRON MLNTAi  data*) 

5290 

gsO 

FORMAT 

(*OOb  temp •/ 1 2F 1 0 . 2/ 1 2F 1 0 . 2/ * QmeAN  V A L U E = » F 1 0 • 3 ) 

5300 

9 60 

forma  T 

( *0T1  */12Fn.2/12F  ln.2/»0MEAN  VALUE**F10.3) 

53  10 

970 

FORMAT 

( »OWLlTE*/(  12F1U.0)  ) 

5320 

9 p n 

format 

(•0DEUUP*/(1?F1D,3)) 

5330 

9 90 

format 

( *OWOCUP*/(  12F10.3)  ) 

5390 

1 n 0 n 

FORMA  T 

( I 1 0 , F 1 0 . 0 ) 

5350 

inn 

F 0 R A T 

( 2MF5 . n ) 

5360 

r^2n 

format 

{1H1?0X*ExP0SURE  surface  temperature,  degrees  F for  *6A6// 

53  70 

] 7 X ♦ T 1 Mt 

»SX*(1)*7X*{2)*7X»(3)»7XM9)*7X*(5)*7X*(6)*7X»(7)»7X»(8)*7X 

5380 

2 • ( 9 ) • 7X 

• C8  * ) 

5390 

n 30 

format 

( I 1 U , 1 OF  1 n . 2 ) 

5900 

1090 

FORMAT 

{1H130X6A6///7X*TiMe»1MX*EXPOSUPE  HEaT  FLuX  ’BOX*  HEAT  GApNl 

59  1 0 

IS  » 3 X • 

sensible  LOAD’IX*  Total  LOAD*/30X*BTU/hRiFT2*9SX*BTU/Hr*12X 

5920 

2 • P T U / HR 

♦SX»BTU/HR*/*rj»aix»sEN5lBLE*3x*LATENT*/19x*  ( 1 ) »9Xt(2)  ’MX*  (3 

5930 

3 ) ♦ M X * ( 9 

)*9Xt(S)*MX*(6)*9x*(7)*9x*(8)*MX*(9)»7X*HEAT*7X»HEAT»/) 

5990 

1 nbn 

FORMAT 

( I lU,9F7.2,9X,9Glr.9) 

5950 

1060 

FORMAT 

(1H120X*  INSIDE  SLRFACE  TEMPERATURE,  DEGREE  F FOR  *6A6//7x» 

5 9 60 

1 T I ME  ♦ S X 

•(l)*7X»{2)'7X*(3)*7X*{9)»7X*(5)*7X*{6)*7x*(7)*7X*(8)*7X*(9 

5970 

2 ) • 7 X » T A 

• 8 X • A R ♦ / ) 

5980 

10  7 0 

FORMAT 

( I 1 n , 1 1 F 1 0 . 2 ) 

59  90 

1080 

F'  0 K M A T 

(♦1CINVER5IONEKpOrATn=*IIO) 

5500 

1 G9Q 

format 

( • 0 TOS  • ) 

55  1 0 

1 inn 

format 

(•OTIS’) 

5520 

1 1 n 

FORMAT 

( 1 1 U , 1 OF  1 0 . 2 ) 

5530 

1 1 20 

F 0 R M A T 

( 1 H 1 ) 

5590 

1130 

format 

(30X*SITE  SUMMARY  » //7X  • T I ME  * 9X * he  AT  G A I N 5 ♦ 1 5 X • T 0 T A L hEAT'MX 

5550 

1 ’COOLING  1.  0AD»7X»T0tAL*/22x*BTU/HR*  19X»6TU/HR»9X»bTU/HR*7X»c00LING 

5560 

7 LOAD • // 1 4X  • SENS  I BLE  H E A T • 3 X * L ^ T £ N T HEAT*//) 

5570 

inn 

F 0 R M A T 

(////•  heat  LOSS*iOx»COOlING  L0AD*7X*AvERaGE  heat  GAIN*3X// 

5580 

1 1 X 3 ( G 1 0 

. 9 , 1 OX  ) ) 

5590 

1 1 so 

format 

(6(GiO,9,SX)  ) 

5600 

l 1 60 

format 

(////•  HEAT  LOSS’iQX’MAX  HEAT  GAIN*7X*AVEraGE  HEAT  GaIN*3x* 

56  1 0 

I MAX  total  cool  L0AD»//lXt4(Gl0.9,lDX)) 

5620 

1170 

format 

( • n SURFACE  NO  1 T Y PE * 9 X , • 1 HT  • 7 X , • I RF  • 7 X , * ABSP * 6X  , * U * 9X  , • H t 9 

5630 

1 X , * A • 9 X 

,*AAZ*SX,  ’SHADE'S X,tuT*  ) 

5690 

1180 

FORM  A T 

(’0  RADIATION  I'TE'RCHAniGE  FACTORS’) 

5 6 5 0 

1190 

format 

(•nSURFACE  1 2 3 H 

5660 

1 S 

6 7 8 9 10’) 

5670 

1200 

format 

( I 1 0 , 1 0 F 1 G • 3 ) 

5680 

1 2 1 0 

format 

(•0  modified  surface  data*) 

5690 

1 220 

format 

(»i  SOLAR  Data  ( ^Sun/q^3L  ass  ) » ) 

5700 

END 

57  1 0 
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subroutine  ABCD2  (Z,K,L,G,A,B,C,D,NL) 
DIMENSION  AX( 10) ,BX{ 10) tCX( 10) tDX( 10) ,6( lO) 
RF.  ALK(10)»L(10) 

P I = H . • A T A N ( 1 • ) 


PP  = 

P I *0  , b 

00 

bO  1=1, 

NL 

I F 

( G ( I ) ) 

MO 

, M 0 

1 0 

I F 

( Z ) 30  , 

30 

,20 

2 0 

ZQ  = 

S Q R T ( Z / 0 ( 

I ) ) 

ZQL 

= ZQ*L ( I 

) 

CO=SIN(ZQL) 

C 1 =C05 ( 2QL  ) 

S1=C0/ZQL 

S2= ( S 1 -C  1 ) /ZQL/ZQL 
A X ( I ) =C  1 

BX(I)=L(I)/K(I)*Sl 

CX(l)=-ZQL*Kn)/L(I)*CO 

!:)  X ( I ) = C 1 
s 0 TO  b 0 
3 0 A X ( I ) = 1 * 

C X ( 1 ) = 0 , 

D.  X { I ) = 1 . 

•3X(I)=L(I)/K(I) 

G 0 T 0 b 0 
MO  A X ( I ) = I • 

!}  X ( I ) = 1 / K ( 1 ) 
c X ( I ) = 0 . 

0 X ( I ) = 1 • 
bO  CONTINUE 

A = A X ( 1 ) 

3 s B X ( 1 ) 

C = C X ( 1 ) 

OsDX  ( 1 ) 

IF  (NL.LT.2)  CiO  TO  60 

CALL  mult  (AX,BX,CX,DX,A,B,C|D,NL) 
^0  RETURN 

E NO 


I 0 
20 
30 
MO 
bO 
6 0 
70 
BO 
90 
100 
I 10 
I 20 
I 30 
1 MO 
1 SO 
160 
1 70 
1 80 

1 90 
2U0 
210 
220 
230 

2 MO 
2bO 
260 

2 70 
2B0 
290 
300 

3 1 0 
320 
3 30 
3M0 
3bO 
360 
370- 
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SUBROUTINE  ABCDP2  (Z»K,L,G,AP,BP,CP,DP) 

lU 

REAL  K ,L 

20 

PIs4.*ATAn(  !•) 

3 0 

IE  (G)  3C,-iO»lO 

^Ci 

1 p 

PP  = P I /M  , /(a 

bO 

IF  (7)  M0,Hn,20 

60 

?n 

ZOsSQKT  ( Z /G ) 

70 

ZQLs/Q*L 

80 

X=L*L*U,S/G 

90 

RE5=L /K 

1 no 

C 0 = 5 I N ( Z G L ) 

1 i 0 

C 1 = C 0 5 { Z G L ) 

120 

S 1 =C0/ZQL 

1 30 

52=(51-C1  J/ZQL/ZqL 

1 MO 

A P = X • S 1 

1 bO 

^P  = X •RLS  *S2 

160 

CP=X*(bl+Cl)/RES 

170 

OP  = X ♦ 5 1 

1 80 

(,  0 TO  bo 

1 90 

3P 

• 

o 

II 

CL 

< 

2 00 

R P = 0 . 

210 

• 

c 

II 

a. 

220 

0 P = 0 . 

2 30 

GO  TO  b 0 

2M0 

MO 

CONTINUE 

2S0 

X=L*L*0,^/G 

260 

AP  = X 

2 70 

tiP  = X*L/K/3 

280 

CP  = K /L*  X • 2 • 

290 

OP  = X 

300 

GO  TO  b 0 

3 1 0 

bC 

RETURN 

320 

END 

3 30 
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subroutine  DERVT  ( a ,B ,C ,n , AP f bp ,cp ,0P , APP ,BpP ,cpp  ,DPP  ,N ) 10 

DIMENSION  A(N),B(N),C(N),D(N),Ap(N),BP(N),CP(N)iDP<N),AT(lO),BT(lO  20 

1),CT(10),DT(10)|ATT(10),RTT{IO»,CTT(10),DTT(10)  30 

D 0301  = 1, N 

DO  20  J= I , N 50 

I F ( I • to  . J ) 60  To  1 n 60 

AT(J)=A(J)  70 

BT(J)=B(j)  80 

C T ( J ) =t { j ) 90 

DT ( J ) =0 ( J ) 100 

GO  TO  20  I I 0 

in  AT  ( J ) =Ap ( J ) 120 

0T(J)=BP(j)  130 

CT ( J ) *CP ( J ) 180 

DT ( J ) =DP { J ) 150 

20  CONTINUE  160 

30  CALL  MUlT  ( A T , R T , C T , D T , A T T ( I ) , B T T ( I ) , C T T ( I ) , D T T ( I ) , N ) 170 

APP  = A T T ( 1 ) 180 

BPP=RTT(1)  190 

CPP=CTT ( 1 ) 200 

DPP=DTT(1)  210 

DO  HO  1=2, N 220 

APP  = APP-^aTT  ( I ) 230 

0PP=BPP+BTT(I)  280 

CPP  = CPP-»-CTT  ( I ) 250 

qn  OPP=npP+DTT ( I ) 260 

return  270 

END  280- 
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function  dPF  (PV)  id 

c This  subroutine  calculates  dlw-point  temperature  For  given  vapor  pre  20 

YsLOG { RV ) 30 

IF(PV.GT.0.1fl36)GOT0)0  MO 

DPF  = 7 1 t98  + 2M.  873*Y  + o,8927*Y*Y  E.0 

GO  TO  2U  60 

in  0PFs79.0M7  + 30.S79*Y+  1 .8893«Y*Y  70 

20  return  80 

END  90 
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subroutine  glass  { SHDCF  ,gLTYP  f «3LAZE  .SH6F  ) 

10 

DIMENSION  TR ( 9 ) I SH ( 25  ) 

20 

common  /SOL/  LAT,LONG,TZ^,  vvaZ,WT»CN,DST,LPYR,S(35) 

30 

TR ( 7 ) *S ( 19) 

90 

TR ( 8 ) xGLTYP 

50 

TR ( 9 ) xGLAZE 

60 

CALLTAR(Tr) 

70 

SH ( 1 ) xS  ( 29 ) 

80 

SH  ( 2 ) xS ( 22 ) 

90 

SH ( 3 ) xS ( 23  ) 

100 

SH ( 9 ) x5 ( 1 9 ) 

1 10 

SH ( S ) xO . S 

120 

SH ( 6 ) =0 , S 

1 30 

5H ( 7 ) aO  , 25 

1 90 

SH ( 8 ) *0 , 

1 50 

SH ( 9 ) xO  , 7 

160 

SH ( 1 0 ) * 1 t Q 

1 70 

SH  ( 1 1 ) *SHDCF 

1 80 

SH(12)»TR(1) 

1 90 

SH(  13)=TR(2) 

200 

5H(  19)=TR(3) 

210 

SH(  15)=TR(5) 

220 

SHI  1 6 ) =TR ( 9 ) 

230 

SH(17)=TR(6) 

290 

CALL  SHG  (SH) 

250 

SHGFxSH (18) 

260 

RETURN 

270 

END 

260 
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subroutine  GPF  (UiZLiZ)  10 

D I MENS  I ON  z ( I ) 20 

PI=M,*ATAN(1*)  30 

SQTPIsSURTlPl)  MO 

P I 2 = 2 . /P I SO 

ER=n,00|  60 

DB=0 , 1 70 

A'RITE(6,3U)  80 

ARITE(6,^0)  90 

Z ( 1 ) =2*2L *^ORT ( U ) /SQTP  I 100 

ZZ=Z ( I ) 110 

Z ( 2 ) =Z ( 1 ) * ( SORT ( 2 . ) -2 , ) 120 

00  1 0 K = 3 , 5^0  130 

ZK  = K 1 MO 

in  Z ( K ) sZ  ( l ) • ( SORT  ( ZK  ) -2  . *SqRt  ( ZK- 1 ) ^-SQRT  ( ZK-2  • ) ) 150 

DO  20  K = 1 , bo  160 

20  ARITE(6,SU)K,Z{K)  170 

RETURN  180 

C . 190 

C 2 00 

30  format  (SOHO  response  FaCTORS  FOR  S E M I - I N F I N I T E BEO  , 2l0 

MO  FORMAT  (bOMQ  K Z (< ) ) 220 

50  FORMAT  ( I I 1 0 , 3F  1 0 . b ) 230 

END  2M0- 
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subroutine  MUUT  ( a , B iC  ,D  , at  , bt  ,ct  ,DT  , N ) ID 

dimension  a ( N ) ,B ( N ) , C ( N ) , D ( N ) 20 

ATT=A ( 1 ) 30 

BTT=B ( 1 ) MO 

C TT  = C ( 1 ) BO 

DTTaO ( I ) 60 

I F { N ,LT  , 2 ) GO  TO  20  70 

DO  10  J*2,N  80 

AT®ATT*A  ( J ) ■►BTT*C  ( J ) 90 

RT= ATT*B ( J ) +BTT*D ( J ) 100 

CT*CTT*A ( J ) >DTT*C ( J ) 1 1 0 

DT=CTT*B ( J ) +0TT*D ( J ) 120 

ATTsAT  130 

BTT=BT  IMO 

CTT=CT  IbO 

10  DTT=DT  160 

GO  TO  3q  170 

20  aT=ATT  180 

HT=BTT  190 

CT=CTT  200 

DT=OTT  210 

30  return  220 

END  230 
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subroutine  output  ( DB  , WB  , Db  I t I I WA  , W I , QGX  , CFML  » QL  I TE  , QoCS  . QOCL  , 0 » 10 

1 WAZ , I type  ,NEXP .name ) 20 

dimension  W(l0)iWAZ|i0),lTYPE(l0),NAME(6),QGX(l0)  30 


Q v;  S = n . 

MO 

y 'A  W = 0 • 

SO 

Q A N = 0 . 

60 

Q A E = 0 . 

70 

0 G S = 0 . 

80 

Q G E * n • 

90 

y G w = 0 . 

100 

y G N a 0 • 

1 10 

Q D S = 0 . 

1 20 

Q D A = 0 • 

1 30 

Q D E a 0 • 

1 MO 

QDNsO  . 

1 SO 

A R I T E 

(6,20)  (NAME(I),Ia|,^) 

160 

AR  I TE 

(6,30)  DB.VWB.AA 

1 70 

AR  I TE 

(6, MO)  DBI.WBI.WI 

1 80 

DBOaDB 

-DB  I 

1 90 

/'DeWA- 

A I 

200 

WRITE 

( 6 , SO  ) 

2 1 0 

write 

(6,60)  DBD  , AD 

220 

DO  1 n 

la)  , NE  X P 

230 

y ( 1 ) a - 

U ( 1 ) 

2M0 

l I a I TYPE  ( 1 ) 

250 

IF  (II 

.EQ,3)  Q(I)aQ(I)'*-QGX(l) 

260 

I 'A  A Z a A A z ( I ) 

270 

IF  (II 

• EQ  . 1 ) QROOFaQ ( I ) 

280 

IF  (II 

•Ey,5)  QFLOORaQ(l) 

290 

IF  (II 

• Ey,6)  QFLOORaQFLOnR  + Q ( I ) 

300 

IF  (II 

•EQ,2*AND,IAAZ.EQ.n)  QWS»Q(I) 

3 1 0 

IF  (II 

.EQ.2.AnD#IAAZ.EQ,90)  QWWaQ(l) 

320 

IF  (II 

.Ey.2.AND.lWA7,Ey,-90)  QWE»Q(I) 

330 

IF  (II 

.EQ.2,AnD.  IWAZ.FQ,  i8n)  Q«N  = 0(I) 

3M0 

IF  (II 

.Ey,3.AND,rAAZ.EQ.n)  06S=Q(I) 

3S0 

IF  (II 

. EQ  • 3 , and  . I A A 2 . EO • 90  ) QGWsQ(I) 

360 

IF  (II 

• EQ  , 3 . AND  • I AA  Z . EQ  • -90  ) QC^E'QII) 

370 

IF  (II 

• EQ*3»AND*  IAAZ*FQ«  ] 0Q)  Q&N®Q(  I ) 

380 

IF  (II 

•EQ.M.AND.IWAZ.EQ.O)  QnS*Q(I) 

390 

IF  (II 

• EQ.H,  And.  I AAZtEQ.90  ) QDAaQ(I) 

MOO 

IF  (II 

.EQ,M,AND.IAAZ.EQ,-90)  QD£aQ(I) 

M 1 0 

IF  (II 

• EQ.M.AnD,  IAAZ.EW,  1 80)  QONeQ(I) 

M2Q 

AP  I TE 

( 6 , 3S0 ) 

M30 

A RITE 

(6,70)  QROoF 

MMO 

write 

(6,80)  QWIS 

MSO 

AR  I TE 

(6,90)  QWE 

M60 

WRITE 

(6,100)  QWW 

M 70 

WRITE 

(6,110)  QWN 

M80 

WRITE 

( 6 , 1 BO ) 

M90 

WRITE 

(6,190)  QDS 

500 

WRITE 

( 6 ,200  ) QDE 

5 1 0 

WRITE 

(6,220)  Q D A 

520 

AR  I TE 

( 6 , 2 1 0 ) QDN 

530 

AR  I TE 

(6,120)  QFLOOR 

5M0 
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WRITE 

(6,130) 

SbO 

AR  I TE 

(6,180)  QGS 

B60 

WRITE 

( 6 , 1 bO ) Q6E 

570 

WRITE 

(6,160)  QGN 

580 

WRITE 

(6,170)  Q G W 

590 

Q 1 NF  I L 

« 1 • 0B*CFML  *DBD 

600 

vVATTxQLITE/3.81B 

6 1 0 

WR  I TE 

(6,380) 

620 

SuM  = QROOF-^WAS4-QAE-*-OAA4QAm  + qploOR'*-QDS  + QdE*^QDw-^QDN4-QGS'»-QGA-^(>|GE*pGN 

630 

WRITE 

(6,270)  SUM 

680 

WRITE 

(6,230) 

650 

write 

( 6 , 250  ) WATT  , QL I TE 

660 

WRITE 

(6,260)  QOcS 

670 

write 

( 6,280)  CFmL , DBO  , 0 I Np  I U 

680 

WRITE 

( 6 , 380 ) 

690 

sum=sum+qlite+qinfil+qocs 

700 

WRITE 

(6,280)  SUM 

710 

Q I NF  I L 

= 8880 • *WD  *CFML 

720 

write 

(6,320)  CFmL,WD,QU:Fil 

730 

WR  I TE 

(6,330)  QOCL 

780 

WR  I TE 

( 6 ,380  ) 

750 

SUMLxQlNFIL+QoCL 

760 

WRITE 

(6,2*70)  SUML 

770 

SUMTxSUM+SUML 

780 

WRITE 

( 6 , 300  ) 

790 

WRITE 

(6,310)  S U M T 

800 

return 

8 1 0 

c 

820 

c 

830 

?0 

format 

( 1 H 1 1 OX • summary  Of  CALCULATIONS  F0R»6A6) 

880 

30 

format 

(*o  outdoor  COND I T I ons • • , 

• . »F5, 1 , • 

DB*F10,l,*WB«Fl0.8i*HUMl 

850 

) DTY  ratio*) 

860 

HQ 

FORMAT 

(*0  SPACE  COND I T I oNs  • t . • • 

. . *F5,  1 , • 

DB'FIO,  1 , *WB*F10,8  , *HUmI 

870 

1 D T Y RATIO*) 

860 

BO 

FORMAT 

( 28X  * • , 20X  • • ) 

890 

60 

FORMAT 

( * 0 D I FFERENCE  • . • 

•••*F5,l, 

F25 . 8 ) 

900 

70 

format 

( • roof  = ♦ 36X I F 1 0 . 0 ) 

9 1 0 

80 

FORMAT 

(•  SOUTH  WALL=*36X,FlO,U) 

920 

90 

FORMAT 

(*  East  wall  =*36XiFio,0) 

930 

ion 

format 

(•  NORTH  WALL=*36X,F10,U) 

980 

1 lO 

FORMAT 

(•  WEST  wall  =*36x,F10.0) 

950 

1 20 

FORMAT 

( • Floor  = • 36x ,f i Q«o  ) 

960 

1 30 

FORMAT 

(*0  SOLAR  heat  GAtN  aND  TRANSMISSION  THROUGH  GLASS*) 

970 

1 80 

format 

( * South  ■ • 36x , f i o , o ) 

980  i 

1 BO 

FORMAT 

(•  east  =*36X,F10.0) 

990 

1 60 

format 

( * NORTH  « • 36X  , F 1 0 • 0 ) 

1 000 

170 

FORMAT 

( * WEST  =*36X  »F10,0) 

1010 

1 80 

format 

( * qooors  • ) 

1020  1 

1 90 

format 

(*  South  s»36XiF10,U) 

10  3 0 

200 

FORMAT 

( * East  =*36x ifio.O) 

10  8 0 : 

2 1 0 

FORMAT 

( • north  **36X ,F10.0) 

10  5 0 

220 

FORMAT 

(*  WEST  =*36X,F10#0) 

1060  1 

230 

format 

( • G internal  load  * ) 

1 0 7 0 ! 

280 

FORMAT 

(*  I NF I LTRAT I ON  * F 1 0 , 0 » • CFMX  1,08  X 

*F10. 1 ,F13*Q) 

1080  ! 

250 

FORMAT 

( • LIGHTS'FIO.  1 , »x  3.81-* 

23X  ,F 1 0*0 ) 

1090 

260 

FORMAT 

(*OPEOPLE  « »37X,F10#0) 

1 100 

270 

FORMAT 

( 8 8X , F 1 0 • 0 ) 

1110 

280 

FORMAT 

(*o  total  sensible  Space 

LOAD ♦ 20X 

,F10,0) 

1 120 

290  FORMAT  (*0  TOATL  LATENT  SPACE  U 0 A D • 2 0 X , F I 0 . 0 ) I i 30 

300  FORMAT  -•  llHO 

] ) I 1 SO 

310  FORMAT  (•  GRAND  TOTAL  lOaO  •3qX,F10.0)  1160 

320  FORMAT  (♦0INFILTRATI0N»F'^.1,»CEm  X 48*40  X»F6,*4,**»lOX,FlO.O)  1170 

330  format  (•  PEOPLE  = •27X,F20.0)  1180 

3*40  FORMAT  (50X,» *)  1190 

END  1200- 
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subroutine  OUTSID  (X,Y,Z,Cr,UX»FO,DB,TIM,QO,QI  ,QSUN,QSKY,TO,Ti  ,TOn  10 

1 EW  , T A , I TEMP  ) 20 

DIMENSION  TO ( 1 ) , T I ( 1 ) , X ( 1 ) , Y ( I ) , Z ( I ) 30 

XNUMsQSun-USK  Y-fFO*  ( DB-T  I M ) MO 

I F ( X ( 2 ) ) bO  , 1 0 , BO  BO 

1 0 I F ( FO ) 20  » 20 , 30  60 

20  TONEWsTO ( I ) 70 

GO  TO  Hq  80 

3nTAMrTA-TlM  90 

T0NEWs(XNUM4.UX*TAM)/(UX  + r0)  100 

MO  CONT  I NUE  1 1 U 

QO*U X * ( T A M- TONEW  ) 120 

IF  ( I temp  . EQ  , 0 ) 0 1 «Q0  1 30 

TO ( I ) sTOnEA  I mo 

return  ISO 

50SUMZ=0t  160 

SUMYs Y { 1 ) *T  I ( 1 ) 170 

SUMXsX  ( I ) • T I ( 1 ) 180 

SUMXYsO.  190 

DO  60  Js2,Hfl  200 

SUMY=SUmy+Y(J)*TI(J)  210 

SUMXaSUMX>*.X(J)*Tl(J)  220 

SUMX Y = SUM X Y + Y ( J ) * TO ( J- 1 ) 230 

60  SUMZ=SUmZ+Z(J)*T0|J-1)  2M0 

XNUM  = SUMY-SUMZ-»-CR*QO  + XNUm  2B0 

TONEWaXNUM/ ( Z ( 1 ) ♦FO  ) 260 

IF(F0)70»70,80  270 

70TONEWaT0(l)  280 

00  TO ( 1 ) xTONEW  290 

SUMZxbUMZ + Z ( 1 ) *T0 ( 1 ) 300 

5UMXY=SuMXY>Y(1)*T0(1)  310 

Q0=SUMY-SUNZ+CR*Q0  320 

IF(ITEMP,EQ*0)QI*SUMX-sUmXY^CR*QI  330 

RETURN  3M0i 

END  3B0i 


I 
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subroutine  PSYl  ( DB  , WB  I Pp  , DP  I p V , W ,H  , V , RH  ) lO 

This  subroutine  calculates  vapor  pressure ( pv  ),  hum i d i ty  ratio  (W)  20 

enthalpy ( H ),  VOLUME { V ) jRELAT  I VE  HUMIDITY(RH)  AND  DDEW-POINT  30 

temperature  when  the  DRY-BuLB  TEMpERATURE(DB), wet-bulb  TEMpERATuR  MO 

(WB)  AND  barometric  PrESSuRE(PB)  ARE  GIVEN  SO 

PVPsPVSFi'Mb)  60 

I F ( DB-AB ) 30 , 30  , 1 0 70 

Q A'5TAR  = 0.6?2*PVP/(PB-PVP)  80 

I F ( WB-32  . ) 20 , 20  I MO  90 

>0  P V = PVP-S , 70mE-M*PB* ( DB-AP ) / 1 , 8 100 

G 0 T 0 5 0 1 1 0 

JO  P V«P VP  120 

G 0 T 0 5 0 1 3 0 

JO  CDBa ( DB-32 • ) / 1 • 8 1 MO 

C W B 3 ( A B - 3 2 • ) / 1 . 8 ISO 

HL=597.31+0,MM09*CDR-CWB  160 

CH=0.2M02+0,MM09*ASTAR  170 

EX=(WSTAR-CH*(CUB-CWB)/Hi_)/n.622  180 

PV=PB*EX/ ( I . +EX ) 190 

bO  AsO  • 622*Pv/ ( PB-PV ) 200 

Vs0.7SM*(0B  + MS9.7)*(l+70n0*A-/M360)/pR  210 

H = 0.2M*DR^(l06l+0.M‘4M*DB)*'rt  220 

DPsDPF ( P V ) 230 

rh  = PV/PVSF ( OB  ) 2M0 

RETURN  250 

END  260 


subroutine  PSY2  ( D B , D P , P R , B » P V , KX  , H , V , R H ) 10 

c This  subroutine  calculates  the  FOlLOXXINGS  when  DRT-BULB  temperature  20 

c (OB), DEW -POINT  TEMPERaTURE(Dp),ANd  BAROMETRIC  PRESSURE(PB)  ARE  GIVEN  30 

C AB  WET-BULBTEMPERAtURE  RO 

c A HUM  IDITYRaTIO  SO 

C H enthalpy  60 

C V VOLUME  70 

c pv  vapor  pressure  80 

C RH  relative  humidity  90 

I F ( DP-DB  ) 20  , 1 0 , 1 □ 100 

1 D DpsDB  110 

20  P V=PVSF ( DP ) 120 

PV=PVSF(DP)  130 

P VS*PVSF ( DB ) 1 MO 

RH=PV/PVS  ISO 

A=0.622*Pv/(PB-PV)  160 

V»0.75M*(DB+M59.7)*(l+70r0*w/Mi60)/PB  170 

H = 0,2M*DB4(1061+O.MMM*DB)*'*(  180 

A B = W B F ( H , P B ) 19  0 

return  200 

END  210- 
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FUNCTION  PVSF  (X)  10 

dimension  A(6)/-7.90298,c;.02808,-I.3aUE-7»lI,3HH|8,l328E-3,-3,M9l  20 

1*49/,B(M)/-9, 09718, -3tB665M,0. 876793, 0,0060273/|P(M)  30 

T= ( X+M59 , 688 ) / 1 . 8 MO 

IF  (T.LT.273. 16)  go  to  in  50 

7=373 • I 6/T  60 

P ( I ) =A ( 1 ) * ( Z-  I ) 70 

P ( 2 ) =A ( 2 ) *LOG 1 0 ( Z ) 80 

Z 1 =A ( M ) * { 1 - 1 /Z  ) 90 

P(3)sA(3)*(lO**Zl-l)  100 

ZI=A{6)*(Z-1)  110 

P(M)=A(5)*(10**Zl-l)  120 

GO  TO  20  130 

10  2*273. 16/T  IMO 

P ( 1 ) =B ( 1 ) * ( Z-  I ) ISO 

P ( 2 ) *B ( 2 ) *UOG I 0 ( Z ) 160 

P(3)*B(3)*(1-1/Z)  170 

P ( M ) *LOg 1 0 ( B ( M ) ) 180 

20  SUM=0  190 

DO  30  I = I , H 200 

30  SUM  = 5UM-^P  ( 1 ) 2 1 0 

PVSF»29 , 92 1 • I 0**SUM  220 

RETURN  230 

END  2m0- 
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subroutine  rest  ( XX , YY . Z7 I I run ) 10 

c this  program  is  developed  by  T*KUSUDA  of  the  national  bureau  of  2C 

c standards  For  CAlCUlTING  ThE  thermal  response  factors  for  30 

C composite  walls , floors , ROOpS , basement  walls  basement  floors  40 

C AND  internal  FURNIShINGS  Op  SiMpLE  SHAPES  SO 

c RESPONSE  FACTORS  ARE  USED  IN  THE  FOLLOWING  MANNER  60 

C X,Y,Z  are  RESPONSE  FACTORS  70 

c Q I sX*T I -y*to*gma  inside  Where  r is  minimum  so 

c qo=y*ti-z*to  outside  Where  r is  maximum  90 

C TI  inside  TEMPERARURE  where  R is  minimum  100 

C TO  OUTSIOETEMPERaTUREWhERERISMAXIMUM  IIQ 

C K thermal  conductivity  120 

c G thermal  DIFFUSIVITY  130 

c L THICKNESS  IMO 

C lN=0  finite  THICK  WALL  150 

c IN*1  SEMI-FINITE  wall  160 

CINs2S0LID0BJECT  170 

c IF  response  Factors  of  the  solid  cylinder  or  sphere  of  homogeneous  ibo 

C PROPETY  are  desired,  TREAT  ThE  PRqBLEM  OF  MULTILAYER  BUT  WITH  THE  190 

c identical  properties  for  All  the  layers  except  the  radius  200 

REAL  K ( 1 0 ) ,G ( I 0 ) ,L ( 1 □ ) ,KG  2l0 

dimension  X(lD0),Y(l00),Z(i00),Cn0),D(10),RES(la),RMK(l0,H)  220 

DIMENSION  RMKG(M),F(100),XX(10D,1),YY{100,1),ZZ(I00,1),FF(100,20)  230 

10  READ  (S,240)  DELTAT  240 

IRUNsO  2S0 

20  READ  (5,230)  NLAYR,IN  260 

IF  (NLAYR.EQ.o)  go  TO  20o  270 

IRUN=IRUN+1  280 

IF  ( NLA YR .GT  • 1 0 ) GO  TO  200  290 

NNL  A YR*NL  A YR-f  1 300 

IF  (NLAYR.EQ.O)  GO  TO  40  3l0 

DO  30  I s 1 , Nl A YR  320 

30  READ  (5,2H0)  L ( I ) , K ( I ) , D ( I ) , C ( I ) , R E 5 ( I ) 330 

IF  ( I N • EQ  . 2 . AND  . I M , E(3 . 0 ) Gq  TO  50  340 

c Readk,rho,andcofgrounhifin*i  3 so 

C followings  are  ground  thermal  conductivity,  density  and  SP.HT  if  360 

C IN=2,  otherwise  the  SAME  PROPERTIES  OF  THE  INTERNAL  SLAB  370 

40  IF  (IN.NE.O)  READ  (5,240)  kG,DG,CG  380 

C AG  THERMALDIFFUSIVITY  GFEARTh  390 

IF  (IN*NE*0)  AG*KG/C6/DG  4D0 

IF  (NLAYR.EQ.O)  GO  TO  lOD  410 

IF  (IN.EQ.2)  READ  (5,330)  ( R M K G ( J ) , J * i , 4 ) <(20 

50  DO  60  I = 1 , nl A YR  430 

60  READ  (5,330)  ( R M K ( I , J ) , J e I , 4 ) 440, 

IF  (IN.EQ.I)  read  (5,330)  ( R M K G ( J ) , J « I , 4 ) *(50 

D09QIs1,NlAYR  46  0; 

I F ( L ( I ) ) 80 , 70 , 80  M7C 

70  G ( I ) »0*  H80 

K ( I ) » 1 . /RES ( I ) M90 

GO  TO  90  500 

80  6 ( I ) *K ( I ) /C ( I ) /D  ( I ) 51C 

90  CONTINUE  520 

100WRITE(6,3b0)  530 

CALL  RESPTK  ( K , L , 6 , A 6 , K G , X , Y , Z , NL A Y R , DE L T A T , NR T , C R , U T , I N , p ) 54C 
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WRITE  { 6 t 220  ) IRON 

S50 

WRITE  (6,360) 

S60 

WRITE  (6 , 2b0  ) 

570 

WRITE  (6,260) 

5bO 

WRITE  (6,210) 

590 

IF  (NLAYR.EQ.q)  go  to  I3n 

600 

IF  ( I N .EG • 2 , AND  , I M , me  . 0 ) WrjTE 

(6,370)  KG  , DG  , CG , ( RMKG ( J ) 

, J * 1 , M ) 

6 1 0 

no  120  I = 1 , NL A YR 

620 

IF  ( L ( I ) ) 1 20 , I 1 0 , 1 20 

630 

1 n 

K ( I ) »0  t 

690 

20 

WRITE  (6,270)  I,L(I),K(I), 

D ( I ) , 

C ( I ) , RES ( I ) , ( RMK ( I , J ) , Jx 1 

.9) 

650 

IF  (IN*£Q.l)  WRITE  (6,370) 

KG,Og,CG, (RMKG(J) ,jxl ,H) 

660 

30 

WRITE  (6,2^0)  DELTAT 

670 

write  (6,280)  UT 

680 

WRITE  (6,300) 

690 

WRITE  (6,210) 

700 

IF  ( I N • NE  , 0 ) GO  TO  1 5Q 

710 

WRITE  (6,310) 

720 

X X ( 1 , IRON  ) =FLOA  T ( NRT  ) 

730 

YY(1,IRUN)=FL0AT(NRT) 

790 

ZZ ( 1 , IRON ) =FLOAT ( NRT  ) 

750 

X X ( 2 , I run  ) =CR 

760 

YY(2,IRUN)=CR 

770 

ZZ ( 2 , I RUN  ) =CR 

780 

XX  ( NRT-»-3  , I RUN  ) sUT 

790 

no  190  n=i,nrt 

800 

XX ( N + 2 , I RUN ) xX ( N ) 

8 1 0 

Y Y ( N4-2  , IRON  ) xY  ( N ) 

820 

ZZ(N-f2,lRtjN)«Z(N) 

830 

JNxN-  1 

890 

MO 

WRITE  (6,320)  JN,X(N),Y(n) 

,Z(N) 

850 

GO  TO  190 

860 

50 

WRITE  (6,380) 

870 

IF  (IN.EQ.l)  GO  TO  l70 

880 

IF  (IN#EQ,2)  GO  TO  170 

890 

XX ( I , I run ) =FLOAT (NRT  ) 

900 

X X ( 2 , I Run  ) *CR 

9 1 0 

X X ( NRT-^3  , I RUN  ) XUT 

920 

DO  160  Nxl,NRT 

930 

JN=N- 1 

990 

X ( N ) x-X  ( N ) 

950 

X X ( N-»-2  , IRON  ) xX  ( N ) 

960 

AO 

■WRITE  (6, 3*^0)  JN,X(N) 

970 

GO  TO  190 

980 

70 

DO  180  Nx 1 , NRT 

990 

JNxN- 1 

1 000 

FF(N+2,IRUN)xF(N) 

10  10 

80 

WRITE  (6,3^0)  JN,F(N) 

1020 

FF ( 1 , I RUN ) =FLOAT ( NRT  ) 

1030 

FF ( 2 , I RUN  ) =CR 

1090 

FF ( NRT  + 3 , 1 RUN  ) xUT 

1050 

90 

WRITE  (6,2l0) 

1060 

WRITE  (6,210) 

1070 

write  (6,390)  CR 

1080 

GO  TO  20 

1090 

:nn 

RETURN 

1 100 

1110 

1 1 20 

■ 
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2 1 0 

format 

( 2H0  ) 

220 

FORMAT 

( 1 OH  I I run.  110) 

230 

FORMAT 

(1017) 

280 

FORMAT 

{ 1 OF7.O  ) 

250 

FORMAT 

(77MO  layer  t(I) 

1 I ) DESCRIPTION  ) 

260 

format 

( 7 7H  NO 

1 OF 

layers  ) 

270 

FORMAT 

(1I6,1F11.3,1F10.3i1F10,2 

280 

FORMAT 

( 58HO 

1 U=1F7. 

3 ) 

290 

FORMAT 

( 89H0 

300 

format 

( 5QH0 

310 

FORMAT 

( 1 20h0  J 

1 

Z 

2 ) 

3 20 

F ORMAT 

(1117,  1F23.8 ,2f 15.8  ) 

330 

FORMAT 

( 8 A6  ) 

380 

FORMAT 

( 8 8H0 

3 50 

FORMAT 

( 2H  1 ) 

360 

format 

(50H0  wall  composition 

3 70 

format 

( 1F27.3,  1F10.2,  lFlO,3f  lOX 

380 

format 

(50H0  J 

390 

FORMAT 

END 

( 1 I 28  , 1 F2  1 .5  ) 

1 1 3U 
1 IMG 
1 ISO 
1 1 60 

K(l)  (I)  C(I)  Rr.S(  1170 

I I 80 
1 1 90 
1200 

1 F 1 0 • 3 I 1 F8 , 2 I 2X  , ^A6  ) 12  10 

thermal  conductance  1220 

1230 

TIME  Increment  dt=ip3,d)  i2mo 
response  FACTORS)  1250 

X Y 1260 

1270 
1 280 
1290 
1 300 

COMMON  RATIO  CR*1F7.5)  1310 

1320 

) 1330 

8A6  ) 1380 

F ) 1350 

1360 

1370 
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subroutine  RESFX  ( X , Y , Z , y X , y Y , ZZ  iNR  ,cr  , ut  , ne XP  ) 10 

DIMENSION  Xx(lOO,lO),YY(lOG,l0),ZZ(inO,lO),X(lO,lOD),Y(lO,l0D),Z(l  20 

10,100)»NR{10),CR(10),UT(10)  30 

DO  I 0 K=  1 , I 0 40 

D01GJ=1,100  SO 

X X ( J , K ) =0  60 

YY ( J ,K  ) sO  70 

1 0 ZZ ( J , K ) =0  80 

CALL  RtSF  ( XX  , Y Y , ZZ  , I RUN  ) 90 

D030<=1,NEXP  100 

I =K  110 

IF  (K.GT.TRUN)  GO  TO  30  120 

X(I,1)*XX(3,K)  130 

Y ( I , 1 ) =Y Y ( 3 , K ) 1 HO 

Z ( I , 1 ) =ZZ ( 3 ,K  ) 150 

NR ( I ) =X  X ( 1 » K ) 160 

CR(I)=XX(2iK)  170 

J J J=NR ( I ) +3  180 

UT(I)=XX(JJJ,K)  190 

NMAXsNR(I)  200 

DO  20  J=2,NmAX  210 

J3  = J-*-2  220 

J2  = J-*-l  230 

X(I,J)=XX(J3,K)-XX(J2,K)»CR(I)  2H0 

Y(I,J)=YY(J3iK)-YY(J2,K)»Cr(I)  250 

20  Z { I , J ) =ZZ ( U3  ♦ K ) -ZZ ( J2  ,K  ) *CR ( I ) 260 

30  CONTINUE  270 

return  280 

END  290 
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subroutine  RESPTK  ( K ,L  » G , Ag  , kg  , X » Y , Z , NL I dt  , nr  , CR  , U , I S ,F  ) lO 

DIMENSION  K(10),L(I0),G(10),X(I00),Y(10O),Z(10O),AP(I0)iBP(IO),CP(  20 

110), 0P(l0),A(10),B(i0),C(ln)|D(l0),ZRl(3),ZR2(3),RB(3),RAP(3), ROOT  3D 

2(100), KA(3»I00),ZRK(3, 100), RX(IqO),RY(IOO),AZ(IOO)»F(100)  HO 

REALK,L,kG  bO 

PI=H,*ATAN(1.)  ao 

M 3 * 3 7 0 

IF(IS.NE.1)G0T010  faU 

ZL»KG/  1 0 , 90 

UY*100./AG/DT  100 

CALL  GPF  (Uy,ZL,AZ)  llO 

IF(I5«EQ.1.AnD.NL.EQ,0)GoT0330  120 

10  CALLABCD2(0.,K,L,G,AX,nX,cX,Dx,NL)  130 

R B ( 1 ) = D X I H 0 

RB ( 2 ) = I . ISO 

RB  ( 3 ) a A X 160 

U* 1 . /BX  170 

D02Glsj,NL  180 

P X aO  190 

CALL  ABCDP2  ( P X , K ( 1 ) , L ( I ) , G ( I ) , A P ( I ) , B P ( I ) , C P ( I ) , D P ( 1 ) ) 2oO 

?0  CALL  ABC02  ( P X , K ( I ) , L { I ) , G ( I ) , A ( 1 ) , B ( I ) , C ( I ) , D ( I ) , I ) 210 

IF(^'L.LT.2)G0T03D  220 

CALL  OtRVT  (A,B,C,D,AP,Bp,CP,DP,APP,BPP,CPP,DPP,NL)  230 

GOTOHO  2R0 

30APP=AP(1)  250 

BPP«BP(1)  260 

CPP=CP ( 1 ) 270 

DPP*OP ( 1 ) 280 

M0RAP(1)=DPP  290 

RAP ( 2 ) *0.  300 

RAP(3)*APP  310 

DO  SO  I * 1 I 3 320 

C 1 *R AP ( I ) /tix/DT  330 

C2=RB ( I ) *BPP/BX/BX/0T  3M0 

ZR2  ( I ) =-C  1 '•’C2  3S0 

SO  ZR 1 ( I ) *-ZR2 ( 1 ) ♦RB ( I ) /B X 360 

l^'RITe  (6,H80)  370 

*\RITF(6,h90)(ZR1(I),I  = i,m3)  380 

ARITE(6,H90)(ZR2(I),I*1,m3)  390 

CRoOTSOFB(P)xO,  MOO 

NMAXxlO  MlO 

TESTMX*‘40  , M2D 

PX=0,00l  M30 

dpo=o.1/dt  mho 

DL  XxO . 000  1 MSO 

N s 0 M 6 0 

R I T E ( 6 , S 0 0 ) 7 0 

60  DL=DP0  M80 

CALL  ABCD2  ( P X , K , L , G , A X , p X , c X , D x , N L)  M90 

70  PXPsPX^Dl  500 

CALL  ABCD2  ( P X P » K ♦ L , G , A X P , B X P , L x P , D X p , N L ) 5l0 

IF  ( ex*BXP  ) 90  , 1 1 0 , flO  520 

80  PXxPXP  530 

BX«BXP  5M0 
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• 90 
' 100 

i 1 10 

I 120 
i 

j 130 
I 

j 1 HO 

I 

i I'^o 


1 60 
170 

I 80 
1 90 


200 


? 1 0 


220 


i 230 


2M0 

2S0 


260 


i 

V 

t 

i 


testx=px*dt 

IF  ( TEST X-1 E5TMX ) 70fl70,l70 

I F ( dl-dlx  ) mo , I MO  , 1 00 

DL=DL/2 , 

GO  TO  70 

IF  ( BX ) 1 30  , 1 20  , 1 30 

RXX=PX 

GO  TO  IBG 

K X X = P X P 

GO  TO  ISO 

AR=ABS(8X/RXP) 

RXX=(PX  + AB*PXP)/(  l,-*-AB) 

N s N ♦ 1 

ROOT  ( N ) sRXX 

IF  (M.GT.l)  OPO=ROOT(N)-rOOT(N-1) 

NR  T = N 

JVRITE  (6,S10)  N,R00T(N) 

PX  = RXX-*-DLX 

tfstx=pxx*dt 

IF  { TEStX-TESTMX  ) I60,l6n,i70 

I F ( N .LT  .NMaX  ) go  to  60 
vV  R 1 T E ( 6 , S 2 0 ) 

IF  ( root ( NRT  ) - 1 00  . ) 1 90  ,i8q,180 

NrTsNRT-i 

DO  250  J J= 1 , NRT 

P XsROOT ( J J ) 

DO  200  J=  1 »NL 

CALL  ABCD2  {PX,K(J),L(J),G(J)|A(J),B(J)|C(J)iD(0)il) 
CALL  ABCDP2  ( p X i K ( J ) i L ( J ) i G ( J ) t A P ( J > » B P ( J ) » C P ( J > i 0 P ( J ) ) 
CALL  ABC02  (PX»K,L,G,AX,t^X,CX,OxiNL) 
IF(NL.LT,2)GOT0210 

CALL  DERVT  ( a ,B ,C ,D , AP ,BP ,CP iDP , APP  ,PPP  ,CPP  , DPP  iNL  ) 

GO  TO  220 
APP= AP ( I ) 

OPPsBP ( 1 ) 

CPP=CP ( I ) 

DPP=DP 1 1 ) 
py*bpp*px*px*dt 
R A ( I , J J ) =0X /PY 
RA{2,JJ)=1«/PY 
R A ( 3 , J J ) = AX/PY 
PZ=PX*0T 

IF  ( PZ-20  . ) 2Mn  » 2*^0 , 230 
R X ( J J ) -0  . 

RY ( J J ) =2S  . E 1 6 

GO  TO  250 

RX ( J J ) =EXP ( -PZ  ) 

RY ( JJ)  * ( 1 .-EXP(PZ  ) ) **2 

A rite  (6,S30)  R00T(JJ)»(‘=A(M,JJ)iM=1,M3) 

DO  260  J Js  I ,NRT 
DO  260  Ms  1 , M3 

ZRl  (M)»RA(M,JJ)*RXUJ)+Zr1  (M) 

11  = 1 
I I I =2 

WRITE  (6.SM0) 

WRITE  (6,5S0) 

IF  ( ZR  I ( 2 ) *LT  ,0  ) ZR1(2)*P« 


5S0 
560 
570 
580 
590 
600 
6 1 0 
620 
630 
6M0 
650 
660 
670 
680 
690 
700 

7 10 
720 
730 
7M0 
750 
760 
770 
780 
790 
800 

8 1 0 
820 
830 
8M0 
850 
860 
8 70 
880 
890 
900 
9 1 0 
920 
930 
9M0 
950 
960 
970 
980 
990 

I 000 

lOlO 
1020 
1030 
I OMO 
1050 
1060 
1070 
1080 
1090 
1 100 
1110 
1 120 
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WRITE  (6,560)  I I , ( 2R 1 ( M ) , Ms  1 ,M3  ) 

1 1 30 

WRITE  (6,560)  I I I , (ZR2(M)  ,M«l  ,M3) 

1 190 

DO  270  Mal,M3 

1 1 50 

ZRK  ( M , 1 ) aZf<  1 ( M ) 

1 1 60 

270 

ZRK ( M , 2 ) =7R2 ( M ) 

1 1 70 

N T = 1 0 0 

1 1 80 

DO  300  N=3iNT 

1 190 

NR  = N 

1200 

DO  2fl0  Ms  1 ,M3 

12  10 

2fl0 

ZRK ( M , N ) =0 . 

1220 

DO  290  MsI,M3 

1230 

DO  290  J Js  1 , NRT 

1 290 

PZ=(RX(JJ)  )**N 

1250 

2 90 

ZRK(M,N)=ZRK(M,N)^PZ*RY(jJ)*RA(M,JJ) 

1260 

WRITE  (6,560)  N,(ZRK(M,N),MsI,M3) 

1270 

IF  (N,LT.5)  Go  To  3oo 

1280 

TE5T1=ZRK(  1 ,N)/ZRK(  1 ,N-l  ) 

1290 

TEST2  = ZRK(  1 ,N-l  )/ZRK(  1 ,N-2) 

1 300 

TEST3=ABS(TEST1-TEST2) 

13  10 

IF  ( TEST3-0. 0000  1 ) 3 1 0,310,300 

1320 

300 

CONTI NUE 

1330 

3 1 0 

DO  320  Ns  1 , nr 

1390 

X ( N ) sZRK ( 1 , N ) 

1 350 

Y(N)sZRK(2»N) 

1 360 

320 

Z ( N ) sZRK ( 3 » N ) 

1370 

CR=TEST2 

1380 

WRITE  (6,570)  CR 

1390 

IF  (IS*EQ.2)  GO  TO  950 

1 900 

IF  ( I 5 • NE  , 1 ) GO  TO  970 

19  10 

330 

IF  (NL#EQ,0)  go  TO  390 

1920 

GF  = 2*KG/SQRT(DT*AG*PI  ) 

1930 

IF  (NR.lt,  50)  60  TO  350 

1990 

DO  390  Js5U,NR 

1 950 

Z J = J 

1960 

390 

AZ(J)=GF*(SqRT(ZJ)-2.*SQpT(ZJ-1.)-^SQRT(ZJ-2,)) 

1970 

NRRsNR 

1980 

GO  TO  370 

1990 

350 

DO  360  JsNR,50 

1 500 

Z ( J-^  1 ) =Z  { J ) *CR 

1510 

X ( J-»-  1 ) ax  ( J ) *CR 

1520 

3 AO 

Y ( J-*-  1 ) a Y ( J ) ♦CR 

1530 

NRR=50 

1 590 

370 

DO  380  Js  1 , NRR 

1550 

3R0 

F(J)aX(j)-Y(J)*Y(J)/(Z(J)’^AZ(JM 

1 560 

NRaNRR 

1570 

GO  TO  HIO 

1580 

390 

DO  900  Js 1 , nr 

1590 

900 

F ( J ) aAZ  { J ) 

1600 

9 1 0 

WRITE  (6,580) 

16  10 

CR  1 a 1 , 

1620 

DO  930  Js 1 , 50 

1630 

CR  = F ( J+  1 ) /F { J ) 

1 690 

TESTCRaABS ( CR-cR  1 ) 

1650 

IF  ( TESTCR-0. 0000  1 ) 990  ,^90,920 

1660 

920 

CR  1 aCR 

1670 

JJ* J- 1 

1680 

9 30 

write  (6 ,590  ) JJ  ,F ( J ) 

1690 

990 

NRs  J 

1 700 

CR«CR  1 

17  10 

GO  TO 

*♦70 

1720 

i/x’  R I T E 

(6,580) 

1 730 

DO  M60 

J = 1 » N R 

1 7H0 

F(J)=X(J)^Z(J)-2.*Y(J) 

i 750 

J J = J-  1 

1760 

A'  R I T E 

( 6 ,590  ) J J , F ( J ) 

1 7 7U 

M70 

return 

1 7b0 

C 

1790 

C 

1 800 

M80 

FORMAT 

(bOHo  residues 

at  Pan 

) 

18  10 

M90 

FORMAT 

( 3F20  • 6 ) 

1820 

soo 

FORMAT 

(SOHO  ROOTS  OF 

B ( P ) aO 

) 

1830 

5 1 0 

FORMAT 

( I 1 U , 1 F20 . 6 ) 

1 8H0 

5?n 

format 

(SOHO  REsODuES 

AT  paROoT  ( N ) 

) 

1850 

530 

FORMAT 

( MF  20 . 6 ) 

1 860 

580 

FORMAT 

(bOHO  RESPONSE 

Factors  oe  finite  slab 

) 

1870 

bsn 

FORMAT 

( 1 ?uho  j 

X ( J ) 

Y ( J ) 

1880 

1 

Z ( J ) 

1890 

2 ) 

1900 

560 

FORMAT 

(1)0, 3F20.6  ) 

19  10 

570 

FORMAT 

{ 1 OHO  CR= 1 F 1 n . 6 ) 

1 920 

580 

format 

(5QH0  J 

F 

) 

1930 

590 

format 

( 1 I 1 0 , 1 F20.5 ) 

19H0 

END 

1950 
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subroutine  RMTMP  (NpXP,Ny  ,TV,CFML,CFMV,R,TIM,TAiTIF,QL,  ITk)  10 

COMMON  /CC/  X ( 1 0 , I DO ) I Y ( 1 0 , I 00 ) , Z ( 1 0 , 1 00 ) , I T YPE ( 1 0 ) , I HT ( I 0 ) , I Rp ( 1 0 20 

1),ABSP(10).UU0),H(10),HI(10),A(I0),UT(10)|T0S(ID,M8),TIS(I0,‘48),G  30 

2(10,10),T0Y(M8),DB(2M),QL  lTE(2‘+),QEQUP(2H)iQ0CPS(2*^),QI(l0),CRn0)  HO 

3,NR(  10)  ,QG«-AS(  10.2H)  , ITHsT  SO 

DIMENSION  AA(2Q.20),BB(2n),TT{20),TlF(20),A2(20»20),B2(20),B3(20),  60 

16SUM(20)  70 

DBNXsDB(nx)-TIM  80 

TU  = T\/-T  I M 90 

NPXP2=NEXP+1  lOQ 

DO  1 0 I = I , NEXP2  I I 0 

RB { I ) =0 , 120 

B 2 ( I ) = 0 . 13  0 

DO  1 0 J= 1 , NpXP2  1 HO 

A2(I,J)=0.  ISO 

10AA(I,J)=D.  160 

ShG=0 . 170 

HSUMsO,  180 

ASUMsO.  190 

ASUMT=0,  200 

D 0701=1, NEXP  210 

SHG  = SHG  + QGLa5 ( I , NX  ) *A ( I ) 220 

ASUMT  = ASUMT-»-A  ( I ) 230 

GSUM { I ) =0  t 2H0 

DO  20  J= 1 ,NEXP  250 

20  GSUM  ( I ) =GSUM  ( I ) ■►G  { I , U ) 260 

IF(ITYPE(I).NE.3)ASUM  = aSijM-^A(I)  270 

HSUMxhSUM+H I ( I ) *A ( I ) 280 

I R=  I RF  ( I ) 290 

CRX=CR ( I ) 300 

I F { X ( I R , 2 ) ) HO , 30 , HO  3 10 

30  X { I R , 1 ) =UT ( I ) 320 

Y ( I R , I ) =UT ( I ) 330 

CRX=D,  3H0 

Z { I R , 1 ) =UT ( I ) 3S0 

HO  AA(I,I)=X(lR,l)-t'Hl(I)4.GSuM(T)  360 

DO  SO  J= I , NEXP  370 

IF  ( I .EG)  • J ) GO  TO  50  380 

A A ( I , J ) =-G ( 1 , J ) 390 

SOCONTINUE  HOO 

AA ( I ,NExP2)=-HI  ( I ) H 10 

SUMYxY ( I R , I ) *TOS { I . 1 ) H20 

SUMX=0.  H30 

00  60  J=2,H8  HHO 

SUMY=SUMY+Y(IR,J)*TOS(I,j)  hso 

60  SUMXxSUMX  + X ( I R , J ) *T I S ( I , J ) H60 

B3(I)=SUMY-CRX*GI(I)-SUMx  H70 

70  AA ( NEXP2  , I ) xA ( I ) *H I ( 1 ) H80 

QLT=QL I TE ( NX ) /ASUMT^R  H90 

D080I=1,NEXP  500 

SHF=SHG/aSUM  510 

IF  ( I type ( 1 ) .EQ.  3 ) 5HF»0,  520 

BO  BB  ( I ) =B3  ( I ) ■►SHF^QLT  530 

AA(NEXP2,NExP2)=-1.08*(CfMl^CFMV)-HSUM  5H0 
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subroutine  SHG  (SH)  10 

D I MENS  I ON  SH ( 20 ) 20 

C SH ( 1 ) * I NTENS I TY  OF  DIRECT  NORMAL  SOLAR  RADIATION  30 

C SH ( 2 ) = I NTENS I TY  OF  DIFFUSE  5KY  RADIATION  HO 

C SH ( 3 ) = I NTENS I TY  OF  GROUND  REFLECTED  DIFFUSE  RADIATION  BO 

C SH(H)=C0SINE  of  INCIDENCE  OF  DIrECT  SOLAR  RADIATION  60 

C sh(S)=formfactorbetweenthe»indowanothesky  70 

c sh(6)=form  factor  betweem  the  Window  and  the  ground  bo 

C SH ( 7 ) =THERMAL  resistance  at  outside  surface  90 

c SH ( 8 ) ^Thermal  resistance  at  the  air  space  (double  glazinG)  loo 

C SH ( 9 ) sThERMAL  resistance  at  the  inner  surface  110 

C SH(10)=SUnLITAREAFACToR  120 

C SH  ( 1 1 ) =SHAD I NG  coefficient  ,NOn-ZERO  VALUE  WILL  BE  GIVEN  qNlY  130 

c when  THE  WINDOW  Is  ShADEd  BY  DRAPES  OR  BLINDS  OR  IF  IT  H^s  IHO 

C AN  InTERPANESEPA(?AtiONOFMORETHAN  1-INCH  150 

C SH { 1 2 ) =TRaNSM I SS I ON  FACTOR  FOR  DIRECT  RADIATION  160 

C SH ( I 3 ) -TRaNSM I SS I ON  FACTOR  FOR  DIFFUSE  RADIATION  170 

C SH ( 1 H ) =ABS0RPT I ON  FACTOR  FqR  DIrECT  RADIATION  (OUTER  PANE)  I 80 

C SH ( I 5 ) =AB50RPT I ON  FACTOR  FqR  DIrECT  RADIATION  (INNER  PANE)  ' 190 

C SH ( 1 A ) =ABS0RPT I ON  FACTOR  FqR  DIFFUSE  R A D I A T I 0 N ( 0 U T £ R PANE)  200 

C SH ( 1 7 ) =A8S0RPT I ON  FACTOR  FqR  DIFFUSE  R A D I A T I 0 N ( I N N £ R PANE)  210 

C SH(18)=50LAR  heat  GaIN  220 

COMMON  /SOL/  LAT,L0NG|TZ^i,WAZ.'AT*CN,DST,LPYR,S(3B)  ^ 230 

REALNI.N0  ^ 2H0 

N I = ( SH  { 7 ) +SH  ( 8 ) ) / ( SH  ( 7 ) ♦SH  ( 8 ) •►Sh  ( 9 ) ) 250 

N0=(SH(7))/(SH(7)"»-SH(8)+SH(9))  260 

D = 5H  ( 1 0 ) *SH  { I ) *SH  ( M ) • ( SH  ( 1 2 ) ■►NO^SH  ( I M ) +N  I *SH  ( 1 5 ) ) 270 

D0= ( SH ( 2 ) *SH ( B ) ♦SH ( 3 ) ♦SH ( 6 ) ) * ( Sh ( I 3 ) +N0*SH ( 1 6 ) *N I *5H ( 1 7 ) ) 280 

I F ( SH ( I 1 ) ) 20  I 1 0 , 20  290 

in  SH ( 1 8 ) =D^DD  300 

GO  TO  30  310 

20  SH ( 1 8 ) » ( D + DD  ) *SH ( 1 1 ) 320 

30  RETURN  ’ 330 

END  3 M 0 - 
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subroutine  50LVP  ( M , N , C , n , X , I ) 

10 

c 

This  is  a routine  for  solving  simultaneous  linear 

EQUAT IONS 

20 

c 

THE  ROUTInL  was  DEVELOPEo  BY  B«A*  PEAVY  OF  NBS 

30 

c 

ROUTINE  fails  WHEN  ANY  OF  THE  DIAGONAL  ELEMENTS 

IS  ZERO 

MO 

dimension  A(lOO,l01),C(I,l),D(l),X(l) 

50 

no  in  I X = 1 , M 

60 

DO  in  I Y=  1 1 M 

7 0 

1 0 

A(IX,IY)rC(IX,lY) 

80 

DO  20  I Z=  1 ♦ M 

90 

20 

A { I Z , N ) =D ( I Z ) 

1 DO 

L=  1 

1 10 

30 

A A = A { L ,L  ) 

1 20 

DO  MO  K.  = L,N 

1 30 

MO 

A(L,K)*A{L»K)/AA 

1 MO 

DO  6 0 K = 1 , M 

150 

IF  (K.Eq.l)  go  To  60 

1 60 

A A A ( K , L ) 

1 70 

DO  50  I A = L 1 N 

1 80 

SO 

A(K,IA)sA{K,IA)+AA*A(LiIa» 

1 90 

60 

CONTI NUE 

200 

L = L+  1 

2 1 0 

IF  (L.LE.M)  GO  TO  30 

220 

DO  7 0 I P X 1 . M 

230 

70 

X(IP)=A(1P,N) 

2M0 

return 

250 

END 

260 

49c 


c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 


subroutine  sun 

DIMENSION  A0(5)/,302,-.0p02, 368. 1717,0. 09  05 /,AI(5)/-22.93.,m9 
17, 2M. 52, -,03^*4,-, OHIO/, A2(5)/-*229, -3, 2265, -I.IH,, 0032,. 0073/, A3(5 
2)/-. 2H3, -.0903,-1. 09, .00?M,.00i5/tBl  (5)/3.851  *-7.351  ,.58,-.00H3,-. 
3003H/,B2{5)/.D02,-9,3912,-,18,0.,O.OOOH/,B3(5)/-.055,-.3361,.28,-. 
H008,-.0006/  ‘ 

COMMON  /SOL/  LAT,L0NG,TZ  r:  ,A/aZ,AT.CN,DST,LPYR,S(35) 

REAL  LAID, long, MERIO,LONr) 

S(  I )s  latitude, 0E6RE£S(+ north, -SOUTH) 

S(2)=  longitude, DEGrEES{+’/*eST, -EAST) 

5{3)s  time  zone  number 

standard  time  Daylight  saving  time 


atlant  I c 

EASTERN 

central 

MOUNT  A I N 

PACIFIC 

5 ( H ) = OA YS ( FROM 
5(5)= 

S ( 6 ) = 

5(7)  = 

S ( 8 ) = 

S ( 9 ) = 


H 

5 

6 
7 
B 

start 


3 

H 

5 

6 


TIMEiHOUR  after  ■... 
da YL  I GHT  SAVING  T I M 

ground  reflectivity 

ri  PAKMrCC:  MMMQPk 


OF  YEAR) 

M I r.  N I G H T ) 

T I ME  INDICATOR 


FROM 

FROM 


S(11)=SUN  RIS 
5(12) =SUN  SET 
5 ( 1 3 ) =C0S7 
5 ( 1 H ) =C05N 
5(15) =COS ( S ) 
5(16)  =ALPHA 
5(17)  =BET A 


SOUTH 
HOR I ZON 


ground  reflectivity 
Clearness  number 

WALL  AZIMUTH  ANGLE,  DEGREeS 

s(1o)=wall  tilt  angle,  degrees  

S(11)=SUN  RISE  TIME  (HOURS  AFT£R  MIpNiGHT) 
5(12)  =SUN  SET  TIME 

D I RECT  I ON 
D I RECT I on 

D I RECT  I ON  _ 

DIRECTION  cosines 


COS  I NES 
COSINES 
COSINES) 

COSINES  Normal  to  surface 


surface 


5(17)  =BET A 
5(18) =G AMMA 

S ( I 9 ) =C0S ( tTA  ) COS  I NE  of  INCIDENCE  ANGLE 
5(20)»S0LAR  altitude  ANGLE 
S(21)=S0LAR  AZIMUTH  ANGLE 

S ( 22  ) =D I FFUsE  SKY  RADIATION  ON  HORIZONTAL 
S(23)=DifFUsE  ground  REFLECTED  RADIATION 
s(2H)=dir£Ct  normal  RADIaTIQN 
S(25)=I0TAL  solar  RaDIATION  INTENSITY 
5(26)=DIFFUsE  sky  RADIATION  INTENSITY 
S(27)sGR0UNd  reflected  D I F F U S E R A D I A T I ON  INTENSITY 
S(28)=SuN  declination  ANGLE  , DEGrEES 
S ( 29  ) =EQu A T I ON  OF  TIME  , HOURS 
S(30)=A  solar  factor 

5(31)=  SOLAR  factor 

5(32)=  solar  factor 

CLOUD  COVER  mOd I F I ER 

intensity  of  direct  solar  radiation  on  surface 
hour  angle, degrfe 

c o 7 7 


5(32) 

S ( 33  ) = 

S ( 3H  ) 

5(35)  n-u 

P I =3.  1 H 15927 
X=2*PI/366**S(H) 


C 1 «COS ( X ) 
C2«C0S ( 2*X  ) 
C3*C0S ( 3*  X ) 


10 
20 
30 
HO 
50 
60 
70 
80 
90 
100 
I 10 
120 
I 30 

1 HO 
150 
1 60 
170 
1 80 
1 90 
200 

2 I 0 
220 
230 
2H0 
250 
260 
270 
280 
290 
300 

3 1 0 
320 
330 
3H0 
350 
360 
370 
380 
390 
HOO 
H 1 0 
H20 
H30 
HHO 
H50 
H60 
H70 
H80 
H90 
500 
510 
520 
530 
5H0 


50c 


5 I «S I N ( X ) 550 

52*SIN(2*X)  560 

S3»SIN(3*X)  570 

DO  10^=1, b 580 

KS=(K-l)+2«  b90 

in  S(KS)=AO(<)+Al(K)*Cl-*'A2(t^)*c2^A3(K)*C3  + 01(K)*5l^B2(K)*S2  + B3(K,*s3  600 

S(29)*S(29)/60.  610 

L ATO  = S ( 1 ) 620 

LONGsS ( 2 ) 630 

MEK  I Ds 1 5*S ( 3 ) 6S0 

LONDsLOnG-MEK I D 650 

Y = 5 ( 28 ) *P I / I 80  . 660 

YY=LATO*P 1/180.  670 

HPs-TAN ( Y ) *TAn ( Y Y ) 680 

TR=  1 2/^' I * aCoS  ( HP  ) 690 

S ( 1 1 ) = ( I 2- T R ) -S  ( 29  ) -fLOND  / 1 5 • 7oO 

5 ( 1 2 ) =2M , -5  I 1 1 ) 7 10 

H=15*(S(5)-l2-fS(3)*S(29)-S(6))-S(2)  720 

S ( 35 ) *H  730 

Sl3  = SIN(YY)*SlN(Y)-*-C05(Yv)*c0S(Y)*C05(H*Pl/l8n.)  780 

5 ( 1 3 ) =5 1 3 750 

hP1s10O,*aC-OS(HP)/PI  760 

X1=AB5(HP1)  770 

X 2sABS ( H ) 780 

I F ( X 1 -X2 ) I 30 , 20 , 20  790 

20  S ( 1 8 ) =COS ( Y ) *S I N ( H*P I / I 8p . ) 800 

S(15)sSQRT(1.-S(i3)*S(13)-S(18)*S(18))  810 

STESTsS ( I 5 ) 820 

STE5TlaC05(H*PI/180,)-TAf(Y)/TAN(YY)  830 

I F ( STEST I ) 80 , 30  I 30  880 

30  5(15) sSTEST  850 

GO  TO  So  860 

80  S ( I 5 ) s-STEST  870 

50  5(20)=ASINJ(S(13))  880 

I F ( 5 { 1 5 ) ) 70 , 60 , 60  890 

60  S ( 2 1 ) *AS I M ( S ( I 8 ) /COS ( 5 ( 2n ) ) ) 900 

GO  TO  80  910 

70  5 ( 2 1 ) =P I - AS  I N ( S ( 1 8 ) /COS ( S ( 20 ) ) ) 920 

80  S ( 20 ) * 1 80 , * S ( 20 ) /P I 930 

S ( 2 1 ) = 1 80 . *5 ( 2 1 ) /P I 980 

S{28)=S(30)*5(8)*S(33)*EyP{-S|3i)/S(l3))  950 

5 ( 22 ) *S ( 32 ) *5 ( 28  ) /S ( 8 ) /S ( 8 ) 960 

S{23)=S(7)*(5(22)+S(28)*S(13))  970 

ATsS ( 1 0 ) *P 1 / 1 80 . 980 

S(16)sC0S(WT)  990 

AA*S(9)*Pl/l  80.  1 000 

S ( 1 6 ) =Cos  ( lAT  ) I 0 1 0 

S(17)sSlN(WA)*SlN(AT)  1020 

S ( 1 8 ) sCOS ( W A ) *S  I N ( ^T  ) 1030 

S(19)=S(16)*5(13)^S(17)*S(18)+S(10)*S(15)  1080 

S ( 38  ) =5 ( 28  ) *S  ( 1 9 ) 1050 

YaO . 85  1 060 

IF  (S(l9)+0,2)  I00,l0n,9n  1070 

90  Y*0 . 55^0 , 837 *S ( 1 9 ) ♦O . 3 I 3 *S ( 1 9 ) •♦2  1080 

100  IF  (S(19))  Il0tll0,l20  1090 

1 1 0 S ( 1 9 ) «0.  1100 

5 ( 38 ) «0,  1110 

1 20  CONTI NUE  1 120 


51c 


5 ( 26  ) »5 ( 22 ) ♦Y 

1 1 30 

S ( 27  ) «S ( 23 ) • ( 1 -5 ( 1 6 ) ) /2 • 

1 1 HO 

S ( 2B  ) «S  ( 3M  ) ♦S  ( 26  ) >*.5  ( 27  ) 

1 ISO 

GO  TO  IBO 

1 1 60 

1 3 0 

DO  1 HO  Js  1 H , 26 

1 1 70 

1 HO 

S ( J ) = 0 • . . 

1 1 80 

5 ( 3H  ) =0 

1 1 90 

1 BO 

RETURN 

1200 

end 

1210 

bUbKUUI iNL 
REAL  ‘ ‘ ■ 


A 1 ( 6 ) /a  . a 1 I BS 
A2  ( 6 ) /O.O  1 636 
A3(6)/0 •01837 
AM ( 6 ) /O  . 09902 
A5(6)/0«01712 
A6 ( 6 ) /O  . 0 I M06 
A7(6)/0-^‘ 

A8 ( 6 ) /O 
T 1 ( 6 ) / 

T2 ( 6 ) / 

T3 ( 6 ) / 

TM(6)/-0, 01218,  1 
T5 ( 6 ) /-0,0 1 056  , 1 .297  1 1 , 


0.7767M,-3.9i4657, 8,5788  1 ,-8 
1.90783,-6,79030,19.37378,- 
1.92997,-8.89139,18.90197,- 
2,359i7,-10.9715,21.29322,- 
3,50839,-13.8639,26,39330,- 
9.15958,-15.0628,27.18992,- 
.55946,-15.9329,26.70568,- 


/O  . 0 1 056 , 9 
/O  , 008  1 9,5 
/O  . 00670 , 5 
/O , 00228 , 0 
/O  , 00  1 23 , 0 
/O  , 0006 1 , 0 
/O, 00035 ,0. 
/-O . 00009 , 0 
/-O  , 000 1 6 , 0 
/-O  . 000 15,0 
/-O.OOO 12,0 
/-G , 0090 1 , 0 
/-O . 00938 , 0 
/-O . 0Q928 , 0 
/-O  , 0090  1 , 0 
/-0.00279 ,0 
/-O, 00  1 92 ,0 
/-n  . 00  1 36  » 0 
’-0,00098 , 0 


15049,-0.27590,0.25618, 
10579,-0.15035,0,06987, 
07717,-0.09059,0,00050, 
057q6, -0.05878, -0.0 1855 
7^050,7.20350,-20.11763 
57818,7.92065,-20.26898 
95797,7.91367,-19,92009 
36698,7.27329,-19.29369 
1 6968,6.  177  15  ,-  1 5.898  1 1 


•23.89896 ,8.17372/ 
•23.8851 8 , 8 .03650/ 
•22.87993, 7 .57795/ 
•21.69106,7.08719/ 
,9.75995 ,-3.89922/ 

, 1 .5  1 798  ,-9.52069/ 

1 2.99  1 6 1 ,-9.83285 / 
12.83698,-9.95199/ 

1 1 .95889  ,-9.59880/ 
».87 1 52 ,-3.73328/ 

' . 80697  ,-2.99959/ 
.00596  , -2 • 28  1 62/ 

- 1 1 , 78092 , 9 ,20070/ 

- 1 8 ,56099 ,6.53990/ 
-22 . 69299 , 7 .89959/ 
-25.  1 I 877  i8. 68895/ 

1,-27.81 955 ,9.36959/ 
,-26.63898,8.79995/ 
,-29.76915,8 .05090/ 
,-22.996  1 9 , 7 .38  1 90/ 
2.05287 ,0.72376/ 

1 . 90040 , 0.483  1 6/ 
0.97138,0.32705/ 

0 . 67666 ,0.22  1 02/ 
-0.12919,0 .02859/ 
0.02759  ,-0.023  1 7/ 
0.067  1 1 ,-0.03399/ 
,0.06837  ,-0.03  191  / 

. lo.6ftflP4.-6, 79585/ 


fUfLJOoj  / f 

, l 9 ,68829 
,19.79706 
, 1 9 ,90969 
,18,75908 
, I 5 ,28302 

. I 1 * 9 74QC: 


— O.  / “DOS/ 

-6,7961 9/ 
-6,6660^/ 
-6,93968/ 
-5.23666/ 
-9.07787/ 


REAL  TD8(6)/-0, 00098, 0.02576, 3.0090 
dimension  TR(9)  ,A(8,6)  ,Tf8,6)  ,Aq(8, 
TR(1)=  TRANSMISSION  FACTOR  , DIRECT 
TR(2)=  transmission  FACTOR  .DI^FOSE 
TR(3)a  ABSORPTION  FACTOR  , DIRECT, 


, 1 1 .92995  ,-9.0778 
9,  16022  ,-3.  1 2776/ 
6.98747 ,-2.38328/ 
TD ( 8 , 6 ) 


10 
20 
30 
40 
50 
60 
70 
80 
90 
1 00 
1 1 0 
120 
1 30 
1 90 
1 50 
1 60 
1 70 
1 80 

1 90 
200 

2 I 0 
220 
230 
240 
250 
260 
270 
280 
290 
300 

3 1 0 
320 
330 
390 
350 
360 
3 70 
380 
390 
900 
9 1 0 
920 
930 
940 
950 


c 

TR ( 9 ) a 

,0  IFF  USE 

, OUTER 

960 

c 

TR  ( 5 ) a 

, DIRECT 

, inner 

9 70 

c 

TR ( 6 ) a 

, 0 I FFuSE 

, INNER 

980 

c 

TR ( 7 ) a 

COSINE  OF 

incident  angle 

990 

c 

TR  ( 8 ) 

aTYPE  OF  glass 

500 

c 

TR  ( 9 ) 

*ID  CODE 

FOR  THE  glazing 

510 

c 

I D 

1 SINGLE 

glazing 

520 

c 

I D a 

2 DOUBLE 

glazing 

530 

DO  1 0 

1 , 6 

590 

53c 


A ( 1 , J ) *A  1 ( J ) 

550 

A ( 2 , J ) *A2 ( J ) 

560 

A { 3 , J ) =A3 ( J ) 

570 

A ( H , J ) =AM ( J ) 

580 

A(5,J)=A5(J) 

590 

A(6,J)=A6(J) 

600 

A(7,J)=A7(^) 

6 1 D 

A(8,J)=Aa(J) 

620 

T { 1 , J ) =T  1 ( J ) 

630 

T ( 2 , J ) ®T2 ( J ) 

680 

T ( 3 , J ) =T3 ( J ) 

650 

T(M,J)=TM(J) 

660 

I { 5 , J ) ( J ) 

670 

T (6, J)=T6( J) 

680 

T ( 7 . J ) =T  7 ( J ) 

690 

T ( 8 , J ) =T8 ( J ) 

700 

A0(1,J)=a01(J) 

7 1 0 

A0(2,J)=A02(J) 

720 

A0(3|J)=A03(J) 

730 

AO ( H , J ) X a04 ( J ) 

780 

A0(5,J)=A0b(J) 

750 

A0(8,J)=A06(J) 

760 

AO ( 7 , J ) xAO  7 ( J ) 

770 

A0(8,J)xA08(J) 

780 

A I ( 1 , J ) xA  I 1 ( J ) 

790 

AI (2,J)xAI2(J) 

800 

AI (3,J)xAI3(J) 

8 1 0 

AI  { H fJ)  ( J) 

820 

AI(B,J)xaI5(J) 

830 

AI (6,J)=AI6(J) 

880 

AM7,J)xaI7(J) 

850 

AI(8,J)=AI8(J) 

860 

TD(1,J)sT01(J) 

870 

TD ( 2 , J ) xjD2 ( J ) 

880 

TO ( 3 , J ) =TD3 ( J ) 

890 

TD(M,J)xTD'^(J) 

900 

TD(B,J)=TD5(J) 

9 1 0 

TD ( 6 , J ) xTD6 ( J ) 

920 

TD ( 7 , J ) =TO  7 ( J ) 

930 

TD ( 8 , J ) =TD8 ( J ) 

980 

PT A = TR ( 7 ) 

950 

LxTR { 8 ) 

960 

I DxTR ( V ) 

970 

IF  (ID.EQ.2)  Go  TO  30 

980 

TR ( 1 ) xT  ( L , 1 ) 

990 

TR ( 2 ) =T ( L , 1 ) /2, 

1 000 

TR ( 3 ) =A ( L , 1 ) 

lOlO 

TR ( 4 ) xA ( L , 1 ) /2, 

1020 

DO  20  Jx2,6 

1030 

TR(  1 )xTR(  1 )-».T(L,J)*(ETA*^(j-.i  ) ) 

1080 

TR(2)xTR(2)+T(L,J)/(J^l  ) 

1050 

TR(3)xTR(3)+A(L,J)#(ETA#*(j-1  ) ) 

1060 

TR  ( 8 ) xTR  { ♦ A ( L , J ) / ( J+  1 ) 

1070 

TR  ( S ) xO 

1080 

TR  ( 6 ) xO 

1090 

GO  TO  50 

1 1 00 

TR ( 1 ) xTD ( L f 1 ) 

1110 

TR ( 2 ) xTD ( L 1 1 ) /2. 

1 120 

54c 


TR(3)»A0(Li1)  1130 

TR ( H ) -AO ( L I 1 ) /2 • j j 

TR  ( 5 ) *A  I ( L » 1 ) i 1 ,^0 

TR(6)«AI(l»1)/2.  11 

DO  MO  J«2  » 6 j j 

XsETA**(j-l)  IIBO 

TK { 1 ) =T  R { 1 ) +10 ( L , J ) *X  1 1 9Q 

TR ( 2 ) =TR ( 2 ) +10 ( L , J ) / ( 1 ) 12U0 

TR  ( 3 ) =TR  ( 3 ) •►AO  ( L , J ) *X  12iq 

TR  ( M ) =Tr  ( R ) ♦ AO  ( L , J ) / ( 1 ) 1220 

TR  ( 5 ) =TR  ( c,  ) ♦A  I ( L , J ) *X  1230 

MO  TR(6)sTR(51+AI(L|J)/(J^1)  12MU 

E>0  TR  ( 2 ) s2*TR  ( 2 ) 1 250 

TR(M)=2*TR(M)  1 2^0 

TR { B ) r2*TR ( 6 ) 1270 

return  12ao 

END  1290 


55c 


function  WBF  (H,PB) 

10 

c 

THIS  PROGKAm  approximates  The  wET-bULB  temperature  when 

20 

c 

enthalpy  is  Given 

30 

IF  (PB-29,92)  10,30,10 

HO 

1 0 

YsLOG ( H ) 

50 

IF  (H.GT,  11,758)  GO  TO  2 n 

60 

ARF=0,60HI'^3,H8Hi*Y-*.1.360Uy*Y  + 0*97307*Y*Y*Y 

70 

GO  TO  loo 

80 

20 

WBF=3D.918B-39,6820Q*Y+2n,588l*Y*Y-l,758*Y*Y*Y 

90 

GO  TO  loo 

100 

30 

'.V  R 1 * 1 5 0 . 

1 10 

PVlsPVSF(A/Bl  ) 

1 20 

W 1=0,622*PV1/(PB-PVI  ) 

1 30 

Xl*0,2H*\fVBl-^{  1061+0, )*W1 

1 HO 

X 

1 

X 

il 

>- 

150 

HO 

A B 2 a W B 1 - 1 

160 

PV2=PVSF(WB2) 

170 

A2=0.622*PV2/{PB-PV2) 

180 

X2=0,28*rfiB2+(106l+0,HMH»vB2)*W2 

1 90 

Y2=H-X2 

200 

IF  ( Y 1 *Y  2 ) 90,60,50 

2 1 0 

SO 

W B 1 = W B 2 

220 

CM 

>- 

II 

>- 

230 

GO  TO  HO 

2H0 

60 

IF  ( Y 1 ) 80,70,80 

250 

70 

0 F = W B 1 

260 

GO  TO  100 

270 

80 

A'  B F a W B 2 

280 

GO  TO  100 

290 

90 

Z = AB5 ( Y 1 /Y2  ) 

300 

ABFa { WB2*Z^WR 1 ) / ( 1 +Z ) 

3 1 0 

1 on 

return 

320 

END 

330 

56c 


subroutine  WF  ( QG » QX  , Q I S , LG  , lx  , L I S I ql ) 

1 0 

c 

THIS  routine  takes  heat  gains  To  heat 

LOSS  BY  weighting  FaCTOr 

20 

c 

QG--HISTORY  OF  SOLAR  HEAT  GAIN 

30 

c 

OX--HISTORT  OF  LONG  WAVE  LENGTH  HEAT  GAIN 

HO 

c 

OIS--HISTORy  OF  LIGHTING  POWER  INPUT 

50 

REAL  QC3(8)  iQX(8)  ,QlS(8)  ,lG(8)  ,Lx(8)  ,LIS(B) 

60 

REAL  AG(8)/0*2060»-0#398fi, 0*2287, -0,02*45,-0,0026,-0, 0006, -0*0Q  0 2,, 

70 

10,0001/ 

80 

REAL  BG(8)/1* 000, -2, 8586, 2, 0078,-0, 5887*0.  i0.,0.,0./ 

90 

REAL  AX(8)/0. 6258, -1.2892*0. 7932, -0.1573,-0. 0003, 0**U,,0,/ 

100 

REAL  BX(8)/1. 000, -2. 0676,  1,3651,-0, 2837, 0.,0,,0.,0*/ 

1 10 

REAL  AlS(8)/0, 2902,-0.  1866,0.  ,0,  ,0.  ,0. 

, 0 * , 0 . / 

1 20 

REAL  BlS(8)/l. 000, -0.8781  ,D.,0.,0,,0,, 

0.  ,0./ 

1 30 

DIMENSION  QZG(8),QZX(8),qZiS(8) 

1 HO 

DO  1 0 L = 2 , 8 

1 50 

UZG ( L ) »LG ( L- 1 ) 

1 60 

OZX ( L ) *LX ( L- 1 ) 

1 70 

Q Z I S ( L ) = L I ( L - 1 ) 

180 

10 

CONTINUE 

1 90 

DO  20  L=2,8 

200 

LG ( L ) sUZG ( L ) 

2 1 0 

LX  ( L ) =OZX  ( L ) 

220 

LIS(L)«QZIS(L) 

230 

20 

CONTI NUE 

2H0 

SIJMAG  = AG(  1 )*QG(  1 ) 

250 

SUMBG=0 , 

260 

SUMAXsAX(  1 )*QX(  1 ) 

270 

SUMBXsO. 

280 

SUMAI5*AIS(1)*QIS(1) 

290 

SUMB 1 S*0 . 

300 

DO  30  L=2,8 

3 1 0 

SUMAG  = SUMAC3^AG(L)*QG(L) 

320 

5UMBG  = SUMBG4-8G(L)*LG(L) 

330 

SUMAX=SuMAX+AX(L)*QX(L) 

3HQ 

SUMBX=5umbX+BX(L)*LX(L) 

350 

5UMAIS  = SUmAIS>*-AIS(L)*QIS(L) 

360 

SUMBIS  = SUMbIS-*-Bl5(L)*LIS{L) 

370 

30 

CONT  I NUE 

380 

LG(  1 )=5UMA6-SUMBG 

390 

LX  ( 1 ) sSUMAX-SUMBX 

HDO 

L I S ( 1 ) *SUMA I S-SUMB I 5 

H 1 0 

QL  = LG(  1 )+LA(  1 )^L.  IS(  1 ) 

820 

RETURN 

830 

END 

HHO 

57c 
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Appendix  D 


Computer  Program  Used  in  Evaluation  for  the  Prototype  Building 


T.  Kusuda 
D.  M,  Burch 

Environmental  Engineering  Section 
Sensory  Environment  Branch 
Building  Research  Division 
Institute  for  Applied  Technology 
National  Bureau  of  Standards 
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[H15  ^^OGRAV]  chFCKS  T^1E  MBS  INSIDEOIJT  HOUSE 
XX»YY,7Z  are  response  EACTOR5  CALCULATE^  rr^j  RPSPTK 
rJHiCH  HAb  BEEN  DEVELOPEO  BY  T.^USlJDA  OF  MBF 

x»Y»£  are  AJGjMENTED  response  factors  to  S^^ORTt^N  THE  CALCULATION 

A(l)  AREA  OF  THE  ROOF  IN  SO. FT 

A(^)  area  of  the  i^ALLS  IN  SO.  Fj 

A (3)  AREA  OF  THE  FLOOR  IN  SQ.  FT 

A(4)  area  OF  THE  INTERNAL  FURNiSHIMSS 

CF'^  AIR  Leakage  in  cu.ft.per  min. 

JO  Overall  he/'t  transfer  COEFFTCEINT  of  the  o^OR  in  BTH/HR»bQ.FT.  »f 
AO  area  of  TH^  OOOR  In  SQ.FT. 

Jw  overall  heat  transfer  COEF^^ICTENT  fo^  WInOOaS 
A/;  total  window  area  SQ.FT 
10  OROJNO  TEMPERATURE  IN  - 

TA  AbS'JMEO  initial  AIR  TEMPERATURE  AT  Tivir  7E00  IN  F 


'-parameter  R::4rS=lO(‘»T=4R.U  = T+l  » V-2^T 

01  PENSION  X(R»S)»Y(R»S)»Z(R»S).XX(S»R)»YY(5»R)»Z^:(b»R)rCR(R)»Q(R)» 
110  ( S)  » T 1 ( S ) » 1 PX  ( T ) » A ( R ) , 5UM(R  ) »TOY(S)»TTV’E(T)rTIX(T) 

3Uul  00  700  K=l,R 
00  70J  jri,s 


XX  ( J»  K. ) rO  . 
T Y ( J r K ) = 0 . 


/ 0 ( ■ Z / ( J » < ) r 0 . 

CALL  RESF  ( IRUNf  XX»  YYrZZOEL) 
00  1 IrlrlRUN 
X ( I » 1 ) rXX  ( ,3»  1 ) 

Y { 1 » 1 ) rYY ( 3»  I ) 

Z ( I » 1 ) rZZ ( 3»  I ) 

CR ( 1 ) =XX ( ?»  I ) 

•xlvAXrxx  ( 1 » I ) 


00  I J=P»NMAY 


X(I.J)r:XX(J3»I)-XX(J2»I)*CR(I) 

Y(1»J)=YT(J3»I)-YY(J2»I)^CR(I) 

1 Z ( 1 » J) =Z^ ( J3» I ) -ZZ ( J2» I ) ♦CR ( I ) 

•0  0 3 WRITE(6»b0b) 

100  r ORMAT ( 1 1F7.0 ) 

3000  READ(bflOO)  RUN 

ARlTE(6'»0n07)  ^UM 
btii)7  ^ ORMA  n IHl  ♦ RUN  NOr’FlO.O) 

IF(RUN)  boon f 6000 » 7000 
/UOO  0EN=0. 

00  102  I=1»T 

1 IME  ( I ) ='JEN 
1U2  0EN=DEN>0EL 

READ  (b»101)  ( TOX ( I ) » I-l » T) 

RE AO  (SrlOl)  (TIX(I)»I=1»T) 

READ  (S»100)  ( A ( I ) » 1=1 » R ) » CFM » UD» AO » UW r AW » TG» TA 

101  F0RMAT( 12Fb. 1 ) 
wRiTECb^'+nno) 

4UU0  fopnaI(30HO  OUTDOOR  TEMPERATURE  CYCLE 
WRITE  (6rfS0  9)(T0X(I)»in»T) 

IF ( T IX ( 1)  ) 8001  » 8002»  8001 
HO 01  WR1TE(6»3003) 

WRITE  (b»609)  ( Tlx (I ) » 1=1 » T) 

_NTI=1 


Id 


{ HOMNSrOE  AIR  TEMPERATURE  CYCL<^M 
o(i')  OPMA  T { 12F10  . ? ) 

UU  1 :)  FOlO 

HDlir; 

bub  ^ r ( pun  ) 

^ ( l^OUO  A(l)  M?)  AC^)  A(t+)  CFM 

1 j)  A1  UW  AW  TG  TA  ) 

M j i 0 C v)‘J  11  J )E 

!/•/  R 1 r F ( b » o 0 b ) 

WRl  TE(n»bO^^)  ( A(  I ) » I— 1 »R)  »CFM»UG»AD»UW»AW»TQ»TA 
JU  SlUi  j::i»T 

( UT I , UE. 0 ) r I ( J) rT I X ( T- J+1 ) 
but)  I 0 ( J ) - T bx  ( T - JM  ) 

JO  bOl  jr'JrV 

1-  ( JT I .U£,0 ) T I { J) rT  I { J-T) 
bjl  ( 'J  ( J ) rrO  ( J-T  ) 

i-(Jn.NE.(M  GO  TO  GOG  4 
JO  0 1 = 1 fV 
G 11(1 ) ztA 
h0ij4  JO  SO 2 J:=1»V 

lO(J)-TO(J)-TG 
sob  i I ( J)=TI  ( J)-r , 

JO  SOj  J=l»T 

IF ( J ( 1 . JE. 0)  T I X ( J) rTlX ( J) -TG 
suS  1 OX ( J) =TOX ( J) -TG 
JO  4 1=1 » 2 
0 { 1 ) ::iJ . 

JO  S J=1»V 

S 0(  I ) = j(  I ) + ( XX ( J + 2r I ) ♦TT ( J) -YY ( J + b» I) tTbl J)  ) 

4 G OUT  I J' IE 
JO  b I=3»4 
0(1)=J. 

JO  b Jr  1 » V 

b 0 { I ) =0(  I ) + VX  ( J^-2»  I ) tTI  ( J) 

SuOl  forma 1(?H1  ) 

aRI  TE ( b»  SOOl ) 

WR 1 TE ( b»  bono ) 

boon  F JRMAT { lOX » ?bu  TEMPERATURFSfF  »3nx»30H  HEAT  FlOaS 

1»HIJ/HR  ) 

WRITE(6»400) 

40  n forv|/\[  ( f T IME  • » GX  • TO*  » GX»  TI  * »bX*ROOF»  »bX*'VALL*  »bX*FLOOR*  »4X*  SOLID 
1 • »bX»oOOR*  »4X* WINDOW*  »SX* INFIL’ »7X*NET* ) 

JO  ?nn  N=i»io 
JO  ?no  N\=1»T 
Lr(N-i)*T4MK 

1F(N7I.'JE.0)  TINEw  = TIX  (NK) 

I NEW=TOX ( NK ) 

JO  10  NTT  = 2r\/ 

10  TOY (NrT)=TO(NTT-l) 

JO  ll  NTT  = 2»\/ 

11  lO(Nln=TOY(NTT) 

DO  1?  NTT  = 2»\/ 

12  I OY ( NTT ) =TI ( NTT-1 ) 

JO  13  NTTr2»V 

13  TI{NTT)=T0Y(NTT) 

10(1) =T  JEW 

IF(NTI.NE.O)  TT(1)=TINEW 
JO  7 l=l»2 

bJM (l)=o(I)*CR(I)-Y(Irl)*TO(l) 


2d 


\)yaxrxv(l,l) 

JO  B Jz?»N'^iAX 

^ SJ’-'K  I ) -S  JM(  I ) f { X ( I » J)  ♦II  ( J)  -r  ( I r J)  ♦TO  ( J)  ) 

f 'o  J'-  ( I ) ZSJ^^  ( I ) ♦ A ( I ) 

J ■)  ^ n I r 3 » 4 

ojv  ( 1 ) -z<{  T ) ♦')  ( 1 ) 

Nj  ;i  A X = X y ( 1 » 1 ) 

JO  ^ J-?»0'^AX 

'kJ’M  1 )rSJV(  I)fv{I,J)^TT(J) 

Jf)  { I ) rS  JV  ( I ) ♦ A ( I ) 

i - ( NT  1 . ‘.Ji  , n ) TO  BOOS 

Jc  l~  )J4  A J+  iw^A  /’+l  , OB^C  - VI 
J^vir  )•- ‘j 
JO  1 h T - 1 » 

1 S JiN  = JI'j  + A ( 1 ) ♦ X { 1 » 1 ) 

‘•♦TO  ( 1 ) 

JO  IB  I r 1 . ^ 

IB  XNE^^.=  X‘jE'^-SJV  ( I ) 

11(1)  -X*I£V/ Jr  M 

«IIUB  JO  17  IzlfR 

17  0 ( 1 ) = J'M  n + A ( I ) ♦ X ( I » 1 ) ♦ T I ( 1 ) 

00= JJ*A 0+ { r I ( 1 ) -TO (1 ) ) 

Ort  = + A A ♦ ( TI  ( 1 ) - TO  ( I ) ) 

0I  = l.nB  + Ct^'^^(TT(l)-T0(n) 

I 0‘l=T0(  1 ) + T G 
I 1 \J=T  i ( 1 ) + tg 
bJV10  = J'^  + Oa  +01 
JO  BUuB  l = l»fX 
Bbuf)  jv'0  = bjvio  + ^ ( I ) 

JO  ?i)i  1 = 1 

l-(A(I))2'UrJ0l»2  0? 

U 2 J ( I ) = 0 ( I ) / A ( 1 ) 

201  COJTlOjr 

I - ( N. . 1 0 ) GO  TO  ?no 

rt^IT£  (B»3nn)  JK  » TON»  TT  N»  { 0 ( I ) » 1 = 1 » 9 ) » 00»  Oi»' » 01 , 5JM0 
0 ION  OJTGISE  air  Tfv^ry/^:^  j:^p 

G TIJ  INSUE;  AI9  TE'^^ERATJPE 

G G(l)»Irlr4)  heat  f^LOX  IN  JTU  =>E9  HO'IRrEO.t^T 

G OJ  TOTOAL  heat  TRANS^E9  through  0009 

G Oa  total  heat  transfer  THROUGH  /JlNJJ'^WE 

G 01  TOTAL  HEAT  TRANjS-ER  OJE  TO  INPILTpaTION 

2J(i  GONTiNjE 
3U0  ‘-'OR'^MT  ( Ib»  1091  0.2) 

GO  10  3000 
bliUO  STOP 
EN  J 
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(.'I  LJlUUrLDlOO 


L 
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C 


ms  r>;nSRAVi  IS  DL'VELO^ED  ry  T.KJSUOA  of  tht  national  bureau  of 

•oT  ANOARn:^  -OR  CALCOLTlNO  HE  THERMAL  RESPO^.ISE  ^ACTORS  FOR 
COMPOSITE  WALLS»  f^LOORSf  ROOt^SHASEMENT  WALLS  RASEMENT  FLOORS 

AN)  Internal  - irnishinss  of  simple  sha^fs 

RES^O  4SE  FACTORS  ARE  USED  IN  THE  FOLLOWiMO  MAN^JER 

X»YrZ  ARE  RESR'ANSF  FACTORS 

'j:)rY*Tl-Z*TO  OIJTSIOE  WHERE  R IS  viaxHD’^^ 

Ti  Inside  temperarjre  where  r is  mintvi'jvi 
TO  OilTSIOF  temperature  WHERE  R IS  MAYlvn)'^ 

K.  thermal  conductivity 

0 THEPVAL  DIPPUSIVITY 


L I HICK  NESS 

l\i-n  OR  '^LABK  plane  WALL 

iM-i  cylindrical  wall 

IM^Z  SRHERICAL  .,'ALL 
iNro  FINITE  THICK  WALL 
IN=1  SEMl-FlNlTE  WALL 
INr?  SOLID  OBJECT 

response  factors  of  THE  SOLID  CYLINDER  or  SPHERE  OF  HOMOGENEOUS 
PRORETY  are  DESIRED»  treat  the  problem  op  multilayer  but  WITH  THE 
JDr.NTICAL  PROPERTIES  FOR  ALL  THE  LAYERS  pYCFdt  THE  RADIUS 
IF  IHEATrO  NO  TEMPERATURE  DATA  THUS  NO  HEAT  C ALC JLATlOiNi 
IF  IHEATrl  PERIODIC  30UNDATRY  CONDITIONS 


RoO  t-ORMATIPHn  ) 

parameter  S=lonrT=10»UrT+l»TVr2*T 
real  K ( T ) r G ( T ) » L ( T ) » kg 

dimension  X(S) »Y(S) »Z(S) »C(T) »D(T) »R(U) »RES(T) »RMK(T»T) »RMKG(T) r 
lF(S)»XX(S»rv)»YY(S»TV)»ZZ(5»TV) 

1 -ORMAT ( 1017) 

D format ( 1UF7.0 ) 

lUO  FORMAl  (14HO  EXPOSURE  NO=  110) 

iul  F0RmaT(77H0  layer  L(D  K(I)  (I)  C(I)  RE'Sd 

1)  DESCRIPTION  ) 

iiiP  porviaT(77H  no 


P OF  layers 
1U3  FjRvia  r ( 1 16»  1F1  ] .3 
iu4  ‘-ORMAKSBHO 
3J=1  >=  7,3) 
lUb  F0RMAT(49Hn 
iOb  FORMAKbUHO 
107  FORMAKiaOHO 

I z 


lFin.3»lP10.2»lF10.3»lFA.P»?Yr4A6) 

THERMAL  CONDUCTANCE 

TIME  INCREMEfJT  DT  = 1F3.1 
response  FACTORS 

J X Y 


IJB  P0RMAT(1I17»1pp3.4»?F1s.4) 

112  -ORMAT { 4A6) 

117  F0RviAT(44Hn  COMMON  RATIO  CR  = lF7,b) 

700  READ(br2)  DELTA! 

1 RUN  = l) 

300  READ(b»l)  NL A YR » NTEST » I M » I N 
IP ( NLAYR.EO. 0 ) GO  TO  SOO 
IRUNr 1RUN4 1 

IF ( NLAYR.GT. 10 ) GO  TO  BOO 
NNLAYRrNLAYR+l 
1 F { NLAYR . EQ.  n ) GO  TO  bOO 
JO  200  1=1 » NLAYR 


) 


) 

) 

) 


4d 


dui)  L ( I ) » K ( I ) » D(  I ) r C ( I)  » RES(  I ) 

iF(  IN.^'J.2.  AM!).  IVI.EQ.O)  GO  TO  ^01 
IF  ( IN.ME.O  ) REAO(^^»?.)KG»OG»CG 
1 F ( Pvl.  ME  . (3  ) AG  = KG/CG/:)5 
IF (ML AYR, E 0,0)  GO  TO  bOl 
1F(  Iv|.EQ«n  ) GO  TO  :^oi 
REAOiM»P)(^(l)r I=1»NMLAYR) 

00  TO  M)?. 

001  R(l)=10. 

JO  01)^  I=2»MMLAYR 
JOO  1 ) =R ( I-l ) +L ( T ) 

00  2 IM.E0*2.  AMO.  Iyi,mF.(3)  REA0(G»112)  (R'^KGf  J)  , J=1  »4) 

JO  115  I=1»MLAYR 
no  E A 0 ( 0 » 1 1 2 ) ( (I  » J ) » J=  1 » 4 ) 

IF(IM.EO.l)  READ(0»112)  ( R vik  G ( J ) » Jr  1 » 4 ) 

JO  10 J I=1»MLAYR 

1 F ( L ( I ) ) 1 1 0 » 1 11  » 1 1 n 

ill  G { 1 ) =0  , 

■<v  ( 1 ) =1  ,/RES(  I ) 

GO  TO  109 

110  G(  I ) r<  ( I ) /c  ( 1 ) /!)(  I ) 

109  CO  IT  1 IJE 

Di)]  GMAr  ( R ( MMLAYR  ) / R ( 1 ) ) ♦*TVI 

CALL  R^SPTK  (K»L»R»G»  AG»KG»X»  Y»?»MLAYR»OELTAT»MRTrCP»  J1  » I » I N r F ) 

XX  ( 1 » I'-^UM)  -FLOAT  ( MRT  ) 

Y Y ( 1 » IRJM) =FLOA  T ( MRT ) 

^2 ( 1 » IRU  J) =FLjAT ( MRT ) 

XX ( ?» IRUM) =CR 

YY  (2»  IR'JM) -CR 
.^2  ( RrlR.JM)  =CR 
XX ( MR  r + 0» IRJM) r JT 

Y Y ( MRT  + 0»  IRJM)  riJT 
^^2  ( MRT4-0»  IRUM)  rUT 
^RlTE(GrlOn)  IRUM 
iF(lM.FQ.O)  WRlTE(8r70l) 

701  f^ORYAKSOHO  PLANE  WALL  ) 

lF(l'^.rQ,])  WRlYE(6»7n2) 

ru?  -ORMAKGUMO  CYLINDRICAL  WALL  ) 

IF  ( I vi.iTQ,  2 ) iA/Rl  TE  ( 6»  70  X) 

/UO  FORVIAT  ( SOHO  SPHERICAL  WALL  ) 

WRITE(6» 101 ) 

WR1TE(6»102) 

WR I T£ { S»  LOO ) 

IF ( NLAYR.EQ, 0 ) GO  TO  S02 

IF  ( IM.E0,2. AMD. IM.NE.O)  WRITE  (SrlRO)  K G r DG » C G » ( R MK G ( J ) r jr 1 » 4 ) 

JO  202.  jri,NLAYR 
IF ( L ( I ) ) 202  » 203»  202 
203  X ( I ) ro , 

202  WRlTE(SrinO)  IrL{I)»<{T)»D(I)»C(I)rRES(T)r(RVK(I»j),Jrir4) 

1F(IM.EQ«1)  WRITE(6»120)  KG  » DG  » CG  » ( R^IKG  ( J ) » J=1  » 4 ) 

12  0 EORV)Ar(lF2  7.3nFl0.?»lF10.3rinx»4A6) 

SO 2 CONTINUE 

1F( IM.NE.O)  GO  TO  153S 
JO  114  N=1»MRT 
JM=M-1 

XX(NM-2.  IRUM)rX(N) 

Y Y ( N + 2 » IRUM) ry ( N) 


5d 


it  4 LCVJTI  JJE 
L-.:)  TO  >',0  4 
IOJ^'3  CO'JTiM.iE 

( SUHO  J F 

( I \J.EO.  1)  GO  TO  <0909 
GO  TO  qQ9B 
9990  blGN^l.n 
GO  TO  50b 

9'^;»  ^ ( I 'J.  ro.  2 • A'MO.  EO.  0)  30  TO  n997 

GO  TO  qQ96 
qqo/  blGN=-i.u 
vj»0  TO  bOb 
990;,  LOJTIvluF 

00  bOb  N-WN9T 

JOrNJ-1 

X ( 0 ) = - X { M ) 

XX  ( M + 2»  IHUO):iX  (0) 
bob  COOT  PJJE 
0 TO  S 0 4 

bOb  JO  bU9  vjri,NRT 

XX  ( M+2»  IRl)‘l)  zSTCiN^F  ( NJ) 

J\I  = N-1 
bU9  C OOT  1 Oi 

b')b  r O^viAT  (1  124»  1F?1  .b) 
b U 4 C 0 ''1 T I >1 IJ  F 
G TO  b n 0 
o 0 1 1 C 0 I T T ' i i,  I F 

bi)0  bFTUGb 

E 9 0 


6d 


oJ  1 

L 


L 


C 

L 

V 

L 


L 


“D  J-iP>')  JT  I NE  ^ES'^  ( I » XV  » Y Y » ZZ  r DELTA!  ) 

I 'US  r’POoRA'^  IS  DEVElD^'^ED  '3Y  T.K'JSUDA  OF  THF  NATIONAL  BUREAU  OF 
standards  -or  CALCULTING  THE  THER^'^AL  RESPONSE  J^ACTORS  FOR 
LO'^POSITii  sfrALLS»ELOORS»ROOFS»RASEMENT  ^^'ALLS  BASEMENT  FlOORS 
AND  I JTERNAL  FjRniiSHINGS  0^  SIMPLE  SHAPFS 
^^c^y')  jr^c  p-(\CTr)RS  ARE  USED  IN  THE  FOLLOaIING  MANNER 
x,Y»Z  ARE  RESPONSE  FACTORS 
Rizx*  fi-Y^TO  + Gv-iA  INSIDE  WHFRE  R IS  ^InTv»iJN 

DlrX  ♦!  l-Y^^TD^GMA  INSIDE  WHERE  R IS  MIv|TvfU‘'^ 

00-Y+II-Z  + TD  OJISIDE  WHFRF  R IS  MAXI'^U'^ 

II  Inside  T'Nperarure  where  r is  '^inthuv^ 

ID  3 U^IDE  temperature  where  R IS  MAvlv^’iJM 
< T U.RMAL  CONDUCTIVITY 
G THERMAL  DIF"USIVITY 
L THICKNESS 

OR  BLAB‘S  plane  WALL 
I‘1=l  cylindrical  WALL 
IMl?  SPH”RICAL  /.all 

iN=o  finite  Thick  wall 
INri  SpmI-FINITF  wall 
ri=2  SOLID  ohuect 

IF  response  factors  of  the  solid  cylinder  or  sphere  of  rIOMOGENEOuS 
pROPETY  are  desired^  treat  the  problem  of  vjltILAYER  t3UT  aITH  THE 
identical  properties  for  all  the  layers  except  the  radius 
1-  UUATpO  no  temperature  data  thus  no  hpat  calculation 

1^  IHEATPI  PERIODIC  BOjNDATRY  CONDITIONS 
4Ul)  ^ORMAT(?HO  ) 

real  K(iu)»G(in)»L(in)»R(in»<5 

J I PENSI  ON  X { 20n ) » Y ( POO ) » Z ( ?nO ) » TI { lOOn ) » TO  M onn ) , c ( 10 ) r D( 10 ) » RES ( 1 
10)  » R M^  ( I 0 U4  ) » R VK  G { 4 ) » F ( 2 0 0 ) 

I F'jPviat  ( lUI  7) 

R format ( lOF 7. 0 ) 
lUO  FORMAT! lOHl  ) 

iOl  F0RMAI(77Hn  layer  L(I)  K(T)  (T)  C(I)  RES(I 

1)  description  ) 


UP  F0RMAT{77h  no 

P OK  layers  ) 

10  3 FORVlAI(ll6»lFn.3rlFl0.3»lFlO.?»lF10.3r  IFS  . 2 » PX  , 4Ab  ) 

1U4  format (S3H0  thermal  CONDUCIANCE 

3U=1F  7, 3) 

iUS  formaT(4RHO  time  increment  DT=1F3.0  ) 

106  FORMAi(SOHi)  RESPONSE  FACTORS  ) 

107  format (IPOHO  U X Y 

1 Z ) 

lUB  FoRmai (1I17»1F23.4»PF15.4) 

il2  format (4A6) 

117  F0RMAT(44H0  common  RATIO  CR=lF7.6) 

RfA[)(6»1)  IHEAT 

iF( iheat.ne.o)  call  TDATA ( T0» TI » Np» THEAT) 

700  RFAD(b»2)  DELTAT 
300  READ(brl)  NLAYR»NTESTr TM, IN 
IF ( NLAYR . GT . 10 ) GO  TO  600 
NNLA YRzNLAYR+I 
IF ( NLAYR . EQ. 0 ) GO  TO  bOO 
■DO  POO  1 = 1 f NLAYR 


7d 


^UU  ^tAD(b»2)  L( 1 ) »K ( I ) 0( I ) »C ( T ) »RES( I ) 

1F(  INI.FO.?.  AN3.  I^.EQ.n)  GO  TO 
bOO  1F(IN,mE.O)  RFA0(5»2)KG»DG»CG 
1F(1N.NE*0)  AGrKG/CG/DG 
IF(NLAYR.EQ.O)  go  to  501 
IF(  Ivi.rQ.O)  GO  TO  301 
READ(b»?)(R(I)»l::i»  NNJLAYR) 

GO  TO  -^0  2 

301  R ( 1 ) =10 . 

DO  303  I=2»NNLAYR 
3U3  R ( I ) =R ( I-l ) +L ( T ) 

302  1F(  JN.FQ.p.AND.I^.NE.O)  RE  AO  ( 5 r 1 1 2 ) ( R^'^’ G ( J ) » J=  1 » 4 ) 

JO  113  I=lrNLAYR 

113  REAr>(b»  1 12)  (R^1K  ( I r J)  r Jrl  ,4) 

1F(1N.E0«1)  read  (5»  112)  (RV1KG(  J)  r J-1  »4) 

DO  100  I=1»NLAYR 

1F(L( I) ) 110» 111 »110 

111  G(I)=0. 

<(l)=l./RE5(I) 

GO  TO  100 

110  G(  I ) =•<  ( I ) /C  { I ) /D(  I ) 

100  COMTINJE 

bOl  GviAr  ( R ( NMLAYR  ) /R  { 1 ) ) ♦♦  TM 
*VR1  TE  ( 6»  20  7) 

2 0 7 FOR'^AnaHl  ) 

CALL  RE5PTK  (K»L^R»G»AG»KG»X»  Y»Z»^JLAYR»nELTAT»^I^T»CR»  JT»  I v) , I N » F ) 
rtRlTE(5»100) 

iF(lM.EQ.o)  WRTTE(6»701) 

701  FORVIAMbOHO  PLANE  WALL  ) 

IF ( n^.EO* I ) WRTTE(5»  702) 

702  PORYAfCbOHO  CYLINDRICAL  WALL  ) 

1F(  iv|.r^0,2)  WRITE(6»  703) 

703  PORMAT(bOrin  E^HERICAL  WALL  ) 

WR I TE ( 5» 101 ) 

v’^RlTE(&»  10?) 

A'  R I T E ( 5 » 4 0 0 ) 

IP(NLAYR.EQ.O)  go  to  '^02 

IFCIN.FO.?)  WRTTE(5»120)  K G » DG  » C G r ( R VIK  G ( J ) » J=  1 » 4 ) 

JO  202  I=lfNLAYR 
1 F ( L ( I ) ) 2 0 2 r 2 0 3 » 2 0 2 
2 0 3 ( I ) = 0 . 

202  WR1TE(5»103)  I»L(I)»K(T)»D(T)»C(I)»RE5(I)»(RMK(I»J),J=1»4) 

IFdN.ri.l)  WRITE(5»120)  KG  » DG  » CG  » ( RMKG  ( J ) » J=1 » 4 ) 

120  FOR''!AT(  1F27.3»  1F10.2»  1F10.3»10X»4A6) 
bO?  WRITE {5» 105)  DELTAT 
WRI TE ( 5» 104 ) UT 
wRlTE(o>106) 
wRITE(5r400) 

1F( IN.mE.O)  GO  TO  1535 
WR1TE(6» 107) 

JO  114  M=1»NRT 
JN=N-1 

114  WRITE(5»10B)JN»X(N)»Y(N)»Z(N) 

GO  TO  504 

1535  WRITE(5»55b) 

bbb  format (SOHO  J F ) 

IF(IN.EO.l)  GO  TO  505 
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1 F ( PJ . EO  • 2 • an  J . I EQ  . n ) 30  TO  305 

00  5()d  N::1  » NRT 
J\I  = N-1 
X { N ) ( N ) 

bUb  WRlTE(b»bO^)  JN»X(N) 
oO  TO  bn 4 

bUb  00  500  N=1»NRT 

b L)  9 w R 1 T E ( b » b 0 ^ ' ) J N » F ( N ) 
bJb  format  ( 1 1?4 » IF?1 .5) 
b()4  ARITE(b»4nO) 
w R I T E ( b » 4 0 0 ) 

WRITE (5^117)  CR 

IF  ( njTEM  . EO.  0 ) GO  TO  300 

CALL  9^AT(X»Y»Z»Tl»TOOELTATfNR»NRT»3V!A»C:R) 
GO  TO  Xf)U 
oun  broR 
- N.) 


o)i  i“L)  Rr^  ArA 

b.JORO  JT  I NE  RES'^T<  ( < » L » R » 3 » A 3 » K G » / » Y » Z » NL  O T » NR  , E R » J » I v, , I b » ^ ), 

C CALCULATES  RESPONSE  FACTORS  BY  MARING  USE  OF  THICKNESS,  THERMAL 

C CONDUCTIVITY,  DENSITY,  AND  SPECIFIC  HEAT  OF  EACH  LAYER  OF 
C COMPOSITE  WALL 

Jl'^E'iraTON  <(10)  »L(10)  »G(10)  »X(10O)  »Y(100)  »Z(10u)  »AP(10)  r 5=*{1 

I U ) » C ( 1 0 ) r OR  ( 1 0 ) » A ( 1 0 ) » 0 { 1 0 ) » C ( 10)  , n ( n ) » 7R  1 ( 3 ) » Z 2 ( 3 ) » R b ( 3 ) r R A ^ ( 3 

2)  fROOI  (luu)  rRA(3»100)  »7R<(3»100)  »RX(10O)  ,RY(100)  »AZ(1i)0) 

3 » F ( 1 .-)(j  ) 

real 

Ri:=3.  I4lbn?7 
Y 3 = 3 

i~  ( IS.E0»2.  ANO.  I^.NE.O)  ^/t3“l 
1 - ( IS. \E . 1 ) GO  TO  613 
oOR  Zl=KG/R ( NL ^ 1 ) 

JY-R(NL-»-1)^*2/AG/0T 
CALL  GPF(jYrZL»  AZ  ) 

1F( IS.EO* I • ANO.NL.EO.O ) GO  TO  901 
ul3  CALL  A9C02  ( 0 . » < » L»  R » G»  AX  » RX  » CX  » OX  r T'-S  NL) 

RR ( 1 ) ^OX 
R T ( 2 ) =1  . 

R R ( 3 ) - A X 
J-1 . /RX 
JO  1 Irl»NL 
RX  = n 

CALL  ARC JP2 { PX  » K ( 1 ) » L( I)  » R ( I ) r G( I ) » AD{ T ) » RD( I ) , CP( I ) 0^( I ) » IM) 

1 CALL  ARCO? (PX »<(!)» L(I)»R(I)»G(I)»A(I)»R(I)»C(I)»0(I)»IM»1) 
1F(NL.LT.2)  go  to  502 

CALL  JERVT(  A»0»C»D»AP»R'^»C^»0^»A^P»RF'^»C^b»OF^»'jL) 

GO  TO  503 
bU2  APP=AP(1) 
rppzRP ( 1 ) 
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CPP=CP( I ) 

□ PPriJP  ( 1 } 

bO.3  IFl  IS.ME.2)  go  to  bOl 
IF  ( IVI.EO.O)  GO  TO  501 
LALL  bOLlD(0.  »R(1)»KG»AG>  IM»hF»HC^P) 

( 1 ) = ( -CPP  + HFP*AX  ) /DX/OT 
( 1 ) r-ZRl ( 1 ) 
mut)  FORMAT  (4F2U. 5) 
bO  TO  1?12 
bUl  RAP(l)=OPP 
R A ( 2 ) r n . 

RAP( i) zAPP 
DO  2 i=lr3 
Ci::RAP(  I ) /RX/OT 
C2  = R5( I ) *BPP/BX/RX/DT 
^:R2(I  )z-Cl+C2 
2 ^R1  ( 1 ) Z-ZR2  ( I ) +RR  ( I ) /F^X 
1212  CO^JTl'lilE 
iOU  poRV'AT(3F2U.6) 

L ROOTS  OF  n(P)=0. 

212  MMAXZi+O 

IF ( IS.ro,2. and. I^.ME.O ) NMAX^lOO 

px=n .uni 

[)PQ=n.i/OT  _ ^ 

1 F ( IS.FO*  2 . AND.  IM.NE.O  ) 0P0=3. 1416^3.l4lG  + AC7/R(l)/RU)*0.2b 

OLXrU.OOOl 

( IS.EO-2-  and.  IM.NF.O)DLX  = OPO/10nO 
N=n 

11  JL  = [)PO 

CALL  AHCD2(PXrK»L»R»G»AXrBX»CX»DXfI'^»NL) 

IF ( IS.EQ*2. AND. IM, NE.O ) CALL  SOLDX (PXf  » ( 1 ) » K G » A G r I M » BX r DX » TEST  1 ) 
15  PX^=PX+DL 

call  ARCD?(PXD.K#L»R»Gr AXP»RXP»CXP»DXPf  Tv»NL) 

IF( IS.NE.2)  go  to  215 
IF ( IM.EQ-0 ) GO  TO  213 

CALL  SOLOX  (PXPrR(l)rKGrA0rIvifBXP»nXP»TF5T^) 

IF ( TEST1+TEST2)  112» 113» 1 14 
114  PX=PXP 

1 EST1ZTEST2 
bO  TO  15 

112  IF(DL-DLX)  130rl3n»117 
117  DL=0L/2. 

b 0 TO  15 

113  IF(TESTl)  llBrllPfllB 
119  RXX=PX 

bO  TO  31 
113  RXX=PXP 
bO  TO  31 

130  A3zA;3S(TEST1/TE5T2) 

PXX=(PX+AB*PXP)/(1+A3) 
bO  TO  31 

213  1F(RX*RXP)  12»13»14 
14  PXZPXP 

BX=BXP 

TE5TX=PX%DT 

IF(TE5TX-100. ) 15»43»43 

12  IF(DL-DLX)30»3n»17 


lOd 


17  DL=DL/2. 
bO  TO  15 

13  1F(0X)  18»19»1R 
19  RXXZPX 
00  TO  31 

18  RXXrPXP 
GO  TO  31 

30  A3=A05(BX/8XP) 

RXX=(PX'»'AB  + PXP)  / ( l.+AB) 

31  N = N|+1 
KOOT(N)=RXX 

IF(N.GT. 1 ) DP0=R00T(N)-R00T(N-1 ) 

NRT  = Ni 

<41  FOR'^ATC  I10»1F20.6) 

PX=RXX4-DLX 

TESTV1XZ40 

TE5TX=RXX*DT 

1F(  TESTX-TESTviX)  <42»42r43 

42  IF(N.LT.NMAX)  GO  TO  11 

43  CONTINUE 

DO  600  JJ=1»NRT 
PX=ROOT( JJ) 

DO  61  J=1»NL 

CALL  ABCD2(PXrK ( J) »L(J)»R(J)»G(J)»A(J)fR(J)»C(J)rU(J)»IM»l) 
bl  CALL  ABCDP2(PX»K( J) »L(J)»R(J)»G(J)»AP(J)»RP(J)rCP(J)pDP(J)»IM) 
CALL  ABCD2(PX»K»L»R»G»  AXrBXrCX»DX» IM»NL) 

1F(NL.LT.2)  go  TO  504 

CALL  DERVT( A»0»CrD» AP»BP»CP» DP »APPfBPP»CPP»DPP»NL) 

GO  TO  505 

504  APP=AP(l) 

3PP=BP( 1 ) 

CPPrCPd) 

DPP=DP( I ) 

505  1F(IS.NE.2)  go  TO  214 
1F(  I'^.EO.O)  GO  TO  214 

CALL  50LID(PX»R(1) r K G » AG » I M r HF r HFP ) 

IF(HF)  401»400»401 

401  PYS  =(hF*AX-CX)/PX/PX/(DPP-HFP^BX-HF*bPP)/DT 
GO  TO  402 

400  PYS=0. 

402  RA(1>JJ)=PY5 
GO  TO  601 

214  PY  = BPP*PX^PX4cDT 
RA(1» JJ)=DX/PY 
RA(2» JJ)=1 ./PY 
RA(3» JJ)=AX/PY 
bOl  PZ  = PX4c[)T 

IF  (P2. LT. 40. ) GO  TO  52 
RX(JJ)  = .0 
RY(JJ)  = 1.E30 
GO  TO  600 

52  RX( JJ)=EXP(-PZ) 

5 RY( JJ)=(1.-EXP(PZ) )^*2 
600  CONTINUE 
54  F0R'^AT(4F20.6) 

DO  154  JJ=1 »NRT 
DO  154  M=1»NI3 


lid 


( V!»  JJ)  *RX  ( JJ)  f ZRK  M) 

ib4  ^R2(^)  = RA  ( VI,  JJ)  ♦RX  ( Jj)  ♦ ( OX  ( JJ) -2.  ) ♦ZP?  (M) 

11  = 1 
1 1 I='Z 

BO  f-’ORViAT  ( bUHO  R^SPO'JSE  FACTORS  OF  FINITE  SLA3 

Bl  FOR'^AT  ( 120H0  J X(J)  Y(J) 

1 Z(J) 

701  f-'OR^Al  ( 120H1  response  FACTORS  FOR  SOLTO  CYLINORICAL  OBJECTS 
1 

/U2  ^'OR^AT  ( laOHl  RESPONSE  FACTORS  FOR  SOLTO  SPHERICAL  uBJECTS 
I 

1P(ZRI  (2)  ,LT,[))  ZRl  (2)=0, 

•JO  67 

ZRK  (VI,  1 )=ZR1  (\i) 
o7  ZRK  ( V|,  p ) ZZR2  ( VI) 
bS  FQRVIA  r ( I io  , 3F2n  . 6) 

NT  = lf)U 

JO  b.H  N = JrNT 
NR  = N 
JO  61 

ol  ZRK  ( vi,'j)zn. 

JO  6 7 v^zl,vi3 
JO  67  JJ=1,NRT 
->Z=  ( RX  { JJ)  ) ♦♦  J 

6 7 ZRK  ( VI,  ) Z7RK  ( vu  N)  4PZ  + RY  ( JJ)  ♦RA  ( V!,  JJ) 
l-(N.LT.b)  GO  TO  66 
I ESTI =7R< ( I ,N) /ZRK (I ♦ N-1 ) 

I EST?  = 7^^<  { 1 , M-1  ) / ZR<  (1  , N-2  ) 

TEST  5 = ARS( TESTl -TEST2) 

IF ( TEST3-0 " 00001 ) 69»bQ,5R 


6 9 

continue 

6 9 

JO  60  ^JZl,NR 

X ( N ) ZZRK ( 1 » N ) 

T(N)=ZRK(2»N) 

60 

Z ( N) =ZRK ( 3 » N) 

CR=TEST2 

o2 

“ORVIAT  ( lOHO 

CR=1Fi0.6) 

1F( IS,FQ.2. ANJ. 

I'^.EQ.O)  GO  TO  ROn 

IF  ( IS.^)E•  1 ) GO 

TO  900 

90  1 

1F(NL.EQ*0)  GO 

TO  905 

Gh=?’«c<3/SQRT  ( JT^AG^PI  ) 

1P(NR.LT.60)  go  TO  610 
JO  P04  J=bO»NR 
ZJ  = J 

EU4  AZ ( J) =GF* ( SORT ( Z J ) -2 . ♦ SOR T ( ZJ-1 . ) +SORT ( Z J-P . ) ) 
NRRZNR 
GO  TO  300 

biO  JO  301  J=NR»bO 
Z(  J + DzZC  J)*CR 
X{JfUzX(J)’»=CR 

3J1  Y ( J+1 )zY( J)+CR 
NRR=6U 

300  JO  206  J=1»NRR 

206  F(j)zx(J)-Y(J)*Y(J)/(Z(J)-^AZ(J)) 

NR=NRR 


DO  904  J=1»NR 
904  P(J)=A?(J) 

9Ub  CONTINUE 

^:^J7  FO-’{'^AT  ( 50HU  J F 

C91=l . 

DO  200  J=l»bO 
L9=F ( J+1 ) /F { J) 

1 ESTC9zAB5(CR-CR1 ) 
if^(TE3TCR-0.0  0n01)  611f6ll»6l2 
bl2  CR1=CR 

2U0  CONTINUE 

209  PORV]AT(llin»lFPn.E) 
ull  NR=J 

CR=CR1 
00  TO  900 
hOO  CONTINUE 

DO  210  J=1»NR 
F(  J)-2*Y( J)-(X( J)^Z( J)  ) 

210  CONTINUE 
900  RETURN 

END 


a)i  ►"OR^^  d»B 

bUBROUTiNE  DER\/T  ( A » R » C f D » AP  r BP  » C ^ » DP  » APf’ » r , DPP  mn|  ) 

: COMPUTES  DERIVATIVE  OF  MATRIX  ELEMENTS  FOR  PLANE  LAYER 

DlviENSION  A(N)»R(N)rC(N)rD(N)rAP(N)fB^(N)rCP(N)rDP(N)»AT(10)rBT(10 
1) ^CT(in) »DT(10)  »ATT(10) »BTT(10)  rCTT(lO) » DTK  10) 

DO  1 1 = 1 r N 
DO  2 J=KN 
IF (I .EQ. J)  GO  TO  3 
AT ( J)=A( J) 

BT( J)=R( J) 

CT ( J)=C( J) 

DT( J)=D( J) 

GO  TO  2 

3 AT(J)=AP(J) 

BT( J)=BP( J) 

CT( J)-CP( J) 

DT( J)=DP( J) 

2 CONTINUE 

1 CALL  viULT  ( AT  » BT»  CT»  DT  r ATT(  I ) » BTT  ( I ) f CTT  ( I)  r DTT  ( I ) » N) 

APPZATT ( 1 ) 

Bi:»P  = BrT(  1) 

CPP=CTT(1) 

DPP=DTT( 1) 

DO  4 I=2»N 
APP  = APP-»-ATT{  1 ) 

BPP=3PPfBTT( I ) 

CPP=CPPfCTT( I ) 

4 DPP  = OPP*»-DTT  ( I ) 

RE  T URN 


END 
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o)l  C»C 

INE  A»B»C»D»  IV)»ML) 

C COMPUTES  MATRIX  ELEMENT  FOR  MULTI-IAYER  PMNE  AS  SHOWN  IN  TABLE  I 
C OF  KUSUDA*S  PAPER 

br''E"'lbI  0^4  AX  ( 1 0 ) » RX  ( 10  ) » CX  ( 1 0 ) » OX  ( 1 n ) f o ( 1 0 ) f G ( 1 0 ) 

JOllBLE  PRECISION  D3E  J » ORE  Y r ZQ 1 » ?Q? 

REAL  K(10)fL(10)»J01»J0R»Jll»Jl? 

Pl=3. 141S927 
K>PZP1  ♦(!  , b 

IP ( NL.LT .?  ) R ( ? ) =R ( 1 ) +L ( 1 ) 

JO  U Ini  » NL 

IF  ( G(  I ) ) lf)3»  ln3»  1 0? 

I:)^  i-(7) 

il)l  ZO:=SQRT  ( Z/G  ( I ) ) 

ZQl -ZjtR ( 1 ) 

Z02  = Z0  + R ( ) 

ZOL=Zj*L( I ) 

IF  ( IVI.NJP.  1 ) GO  TO  3 
J01=r)B'^J(ZQl»l)) 

J1  1='J  3EJ  ( ZOl  » 1 ) 
jn2  = i)RFJ(Z02»  0) 

Jl?=JREJ(Z02r 1 ) 

YUl-JRFY(ZOl»0) 

YllrjRFY(Z01»l) 

Y02nDtiFY(ZO2»0) 

Y12=DRFY(Z02f 1) 

AX(l)--PP*ZQ2+(J01*Y12-YGltJ12) 
bX  ( I ) ::pp  + R (I*-1)/K(I)*(  -YO  1 ♦ JO  2 t JO  1 ♦ YO?  ) 
CX(1)=K(I)/R(I'»-1)*(-J11=^Y12  + Y11*J12)  ♦PP*Z02^Z02 
JX( 1 )=PP*Z02^( J11*Y02-Y1 1 *002) 

GO  TO  4 
3 COrE.lN(ZQL) 

C1=COG(ZOL) 

bl=:CO/ZOL 

b2“{bl-Cl)/ZQL/Z0L 
1F( IM.fO.2)  go  to  S 
AX ( I ) =C  1 

bX ( 1) =L(1 ) /K  ( I ) *51 
CX(I)--ZOL+K(I)/L(I)^CO 
JX ( I ) -Cl 

GO  TO  4 

b GM  — Rtlfl)/R(I) 

AX  ( I)  ::GM*  ( C1-L(  I ) /R  (ITi  ) *S1  ) 

RX(  I)rL(  I)/K(  I)*G\/]  + Sl 

CX( T)-L(I)*L(I)/R(I)/P(I)*<<I)/1-(I)+(-(Z01*Z02t1)*S1+C1) 
JX(I)=Gvit(ClTL(l)/Rn)^51) 

GO  TO  4 
1 AX ( I ) =1 • 

cx ( I) =n . 

JX  ( 1)  = ( R ( T 1)  /R  t I)  ) 1 

IF(IM.F0*0)  RX  ( I ) :pL  ( I ) /K  ( I ) 

ir(ivi,FQ*l)  RX(I)=RnTl)/K(T)+LOG(R(I^l  )/R(I)) 

IF  ( I '>1 . FQ . 2 ) yx  n ) =L  ( I ) /K  ( I ) ♦ ( R ( I T 1 ) /R  ( T ) ) 

GO  TO  U 


n n 


X I +S2/R 

^ 1 = Z 1 
2 

CPzx*(L/«)*^2/PE5^((?.^'^*Rl/L/L+l)^^l-^?P1*7f^2^l.)^b.>) 

■);j::X^^(p1/R*S1  + (L/R)^(^1/'^)*^?) 

bo  ro  u 

T3  AP^Xtbl 

3'orx  ♦ri^rs»S2 
bP-X*(Sl  + Cl  )/^'^5 
JP^X^bl 
bO  TO  4 
103  APrn. 
riP  = n • 

CP-n. 
jp  = o . 
bO  ] 0 4 

101  1F(  I'^.NE-n)  GO  TO  6 
X-L+L^n , b/G 
AP  = X 

3P  = X t L/K/3 
LP=< /l*X*2 • 

JP::X 
b 0 T 0 4 

b 1- ( I^.NE* 1 ) 00  TO  7 

R i::R-»-L 

AP-(n,S*(R^R“*Rl*^l  )4Rl4fRl*L0G(Rl/R)  )^0*0/G 
bP-Rl/4/3/K+(  (Rl^Rl-»-R^R)  ♦LOG  (Rl/R)-(Rl^^l-^*^>^ 
CP-</Rt0.b/G*(Pl^Rl-^*P) 

jPzo«b/G^(0.b*(Rl^Rl“R^R)^^l/^‘"^^^l-*E05(Pl/P)^ 

GO  TO  4 
7 Xi:L^L*n.b/G 
■a::R  + L 

A p = x/  3,^(2^RI/^ 
bPZL/<*Rl/R*X/3* 

CP::</;_*X^L/R*L/R*(2»^R’*^Rl/E/L^0*Eb6G67) 

•jorx/i.*Rl/R^(Rl/R  + ?) 

4 RFTJRM 
EMO 


o)l  t OR  » ♦ Of  0 

bJRROJTI  'JE  A0C'^''*2(  Z,<»Lf  Rf  Gf  APfRPfCPf  DPf 

COMPUTES  MATRIX  ELEMENT  FOR  SINGLE-LAYER  PLANE  AS  SHOWN  IN  TABLE  I 
OF  KUSUDA*S  PAPER 

OOJBLE  PRECISION!  Z 0 1 r 2 02  r DBE  J f DRE  Y 
REAL  X r L f JO  1 f J02 f Jl 1 f J1 2 
PI=3.l4lb927 
li^tG)  100fl03fl04 

104  PPZP1/4./G 

IF (2)  loif inif ]0b 
iOb  ZO  = SORT  ( ^i/G  ) 

ZOL=ZO*L 
Z 0 1 = Z J R 
Z02=Z01 +ZQL 
IF  ( I'^.^^lF.  1 ) GO  TO  3 
XzRt ( R f L ) 
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Y = ( 4-L  ) ♦ *2 

/.l=:(KfL  ) /K 
jni=D‘-;rj(zm  »o) 

jOPr  J -jr  j ( ^'0?  , 0 ) 

Jii“  )‘iEJ(z  n»  1 ) 

] ) 

YO  lrLY:iFY  ( Zivl  f f)  ) 

Yn;^=:Dr^E  Y ( ZQ2*  0 ) 

Y 1 iz  jr-sry  ( Z';l  » 1 ) 

Y1:^  = !-):3EY(Z.'0?»  I ) 

APz(-x+(J]l+YlP-Yll+JlP)+Y+(jni’^Yn?-Yni^jnP))*np 
-jPi:(x*(Jll  + Y()P-YlltJ0P)^Zl/7Q?  + Y*(J01*Yl?-Y01^J12)*Zl/ZQ2)+PP 
♦ ( y ♦ ( J(U  ♦ Y t ?-Yni  ♦ J1  ? )+Y*(Jn  + Ynp-Yll^jn2)) 
jPz(Xt{-J01+Yn:>  + Y01  + J()P)-Y*(-Jll+Yl?-*'Yll*JlP))*PP 

o 0 r i)  4 

3 xzL^L  + n.b/r. 

PES::L/K 
L0  = S1  j ( Z;:JL  ) 

Ll=COS( ZQL) 
bi::C0/70L 

4?=(Sl-Cl)/Zf0L/ZQL 
li^(l^.EQ*n)  GO  TO  b 

-i  / ( T ) z i / < ( I ) 

LX  ( I ) ::0  . 

J X ( 1 ) = ( -<  ( I M ) / '^  ( I ) ) ♦ ♦ T 

4 L:>MTiO.)E 
/\  z A X ( L ) 

d -'-^X  ( I ) 

L-CX(i) 

Jz^X(l) 

1F(\IL.LT.2)  GP  TO  h 
c all  ''^ult  { ax  » » cx  » ox » a » 0 » G f ,0 » NL) 

h KlTUPO 
E'MO 


aJl  - OP  » ♦ L » E 

C ROUTINE  TO  PERFORM  MATRIX  MULTIPLICATION 

GJPPOJTI  4E  '^JLK  A,.P»CO»  ATrRT»CT»DT»N) 
J I mens  n>N  A ( N ) » B ( N ) » C { N ) » D ( M ) 

ATT=A ( 1 ) 

3T  T=B( 1 ) 

CT  rzC  ( 1 ) 

OTTrO( 1 ) 

i" (N.uT.2)  GO  TO  3 
00  1 Jz?»N 

ATZA  r r + A ( J)  +R7  T + C ( J) 

SI  zATT%R(  J)  PBTT’^D(  J) 

CT  = CT  T* A ( J) PDTT  + C ( J) 

JT-CT  r^B( J)+DTT^O( J) 

ATT  = A I 
Sr T=BT 
LI T =C  T 


Ifid 


1 JTT  = [)T 
30  TO  'I 

3 AT  = /VTr 
3T  = BT  r 
CT=CTT 
■JT  = DTT 

4 4ETU^'>I 
END 


oil  ^ t ^ 

bJB^^OjTl  ME  S0LT0(  rKT^r  A3r  iVIrH-  ,HPP) 

C COMPUTES  RESPONSE  FACTORS  FOR  SOLID  MATERIAL 

JO  ible  Precision  orej»zqo 

^0=SQ4T(Z/AG) 

Z01=Z0t^l 

z:l)=z  ji 

= ♦Rl  /AG 
CO  J = ><G/'U 
1F(Z)  P»  l»Z 

P IF ( I M.mE. 1 ) GO  TO  100 
J()l=JBrj(ZOO»0) 

TX=ABS( JOl ) 
li^lTX-n. 1)0001  ) 4r4»5 
b J1 1=DJFJ(Z0D» 1) 

HF=COM*Z01^J11 /JOl 
HFl=Jll/J01/Z01 

HF2=( Jni ♦jni+jil^ ji i-Jol *J1 l/ZOl ) /JOl /JOl 
Hf^  P = -C0M*0 . 5»Z''^  ( H- 1 +H£^P  ) 

GO  TO  300 
100  C=C0S(/01) 

b=SlM(ZOl ) /ZQl 
I X=ABb(SlM(Z01 ) ) 

1 ( TX-0 . 00001  ) 4»4»3 

3 HF  = -CO\'*  ( C/S-1  ) 

HFPz-coM+o .b^ZA^ ( 1+c^ { c-s) /s/s/zoi /zon 

GO  TO  300 
1 HFzO. 

IF  ( i '>l.i^0.2  ) H~P  = -C0M*ZA/3. 

IF(lvi.EO.l)  HPnz-o  , b + COM#ZA 
Go  TO  3 [)  0 

4 HFZO. 

Hcrorn  , 

300  r^ETJ^M 
LMO 


ft)i  rOr(,*  JfJ 

bJBROJTIME  GPF  ( L)»  ZLr  1 vi»  Z) 

C COMPUTES  RESPONSE  FACTORS  FOR  GROUND  HEAT  TRANSFER 
JmEMSTOM  Z(  100)  »ZT(  5000  ) rZS(SOOO) 

JOJBLl  PRECISION  OBEJrOOEYrZQ 
PI=3. 1415927 
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T PI  ( PI  ) 

= n . un  1 
J:3  = n.  1 

I )=2^7LtS0PF(U)/S0TPl 
I ) 

) * { SQPT  { 2 . ) -?  . ) 

J :)  2 K = 3 r ^ 0 

2 2(‘<.)-2(l  >’*'(SQPT(2K)-?,*c;Qr^1'(2K~l  )+S0RT(7t<-P,  ) ) 
i F ( I V! , r , 0 ) 0 f)  7 Q 70 

1 P ( I ''1.  Ft),  j ) GO  TO  1 
{ 1 ) =7  ( 1 ) +ZL 

bO  TO  70 

1 XZFI2  ♦LOOCO.S  + EB  ) •»r) . '^674^691 
b J 'nI  r p I + n , S ♦ ( A T A N ( X ) (1  . F ♦ P I ) 

i X rn 
ri-LP-J- 

JO  17  L = l»OUOi) 

P = J + JU 
B 

1P( IX.ro. 10)  GO  TO  70 
^J0  = 0:3EJ(Z0»n) 

ZYOrOjcTY  ( ZOr  0 ) 

I Fsrx^zjotz JO*  zYO=*^zro 
lESTYzPIZ/P 

1 ESTZ-ARbC  TEST  X -TESTY) 
iP < TESTZ-0 . OOOPl ) 30r30»7l 
31  ^Z-H  + rJ  + B^TESTx 
GO  TO  X? 

30  ZZ-0^!3*P12 
i X = 1 fJ 

3?  ZT (L)=l ./Z7 
LT-L 

FEST-APS(ZT(L) ) *10 

IF  ( TEST-f) . onni ) n > n r 1 7 

17  OOOTIOUE 
n LTY:=LT/? 

LT  X-L  T Y + Z-1 
y ViAXrEf^+  ( LTX-1  ) +00 
B p r 1 . / p V;  A X 
Z J — 1 • / i,  j 
SUT  = SUN]^ZJ 
'B=:EP-JR 
JO  2B  LpWLT 
y=B+DB 

Zy-B*y+ZJ 
rS  ZP~EXP{-ZB) 

ZB  ZS(L)::(  1 . -ZP  ) ♦ -^  T ( L ) 

CALL  sr^s{zs*r)p  »sj^^»ltx) 

b<- ( SJv+bjT ) +PT 2 +BB 
bb  = CA^Pl2 
Z ( 1 ) — GG^ZL^'U 
7()  CO'jriJ.jF 
RE  T JR  g 
EMJ 
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Appendix  E 


Input  and  Output  for  the  Response  Factor  Program 


The  Response  Factor  program  (Appendix  D)  analyses  the  thermal  per- 
formance of  the  inside  space  of  the  prototype  building  under  a prescribed 
outdoor  air  temperature  cycle.  When  the  inside  air  temperature  is  thermo- 
stated,  this  program  calculates  the  rate  of  heat  loss  from  the  building 
at  prescribed  time  intervals.  If  the  inside  room  air  temperature  is  not 
controlled  and  floats  in  response  to  the  outdoor  air  temperature  cycle, 
the  program  then  calculates  time  dependent  variations  of  the  inside  room 
air  temperature.  Following  this  discussion  is  a sample  set  of  data  input 
and  the  print-out  of  corresponding  computer  results. 

A description  of  a sample  set  of  data  input  is  given  below: 

Card  Sequence 

1 Time  increment  of  the  temperature  data  in  hours 

2 Number  of  roof  layers  (includes  the  thermal  re- 


sistances of  the  inside  and  outside  surfaces) 


3 


Thermal  resistance  at  inside  surface  of  roof 


4 


Thickness,  thermal  conductivity,  density,  and 


specific  heat  of  roof 


5 


Thermal  resistance  at  outside  surface  of  roof 


6 


Description  of  the  inside  surface  of  the  roof 


7 


Description  of  roof 


le 


8 

9 

10 

11 

12 

13 

14 

15 

16 

17 

18 

19 

20 

21 

22 

23 

24 

25 

26 


Description  of  the  outside  surface  of  roof 
Number  of  wall  layers 

Thermal  resistance  at  inside  surface  of  wall 
Thickness,  thermal  conductivity,  density,  and  specific 
heat  of  wall 

Thermal  resistance  of  the  outside  surface  of  wall 
Description  of  the  inside  surface  of  wall 
Description  of  wall 

Description  of  the  outside  surface  of  wall 
Number  of  layers  of  the  floor  and  the  semi-infinite 
layer  index  (if  basement  floor) 

Thermal  resistance  at  inside  surface  of  floor 
Thickness,  thermal  conductivity,  density,  and  specific 
heat  of  the  first  solid  layer  of  the  floor  counted  from 
the  inside  surface 

Thickness,  thermal  conductivity,  density,  and  specific 
heat  of  the  second  solid  layer 

Thermal  conductivity,  density,  and  specific  heat  of  the 
earth 

Description  of  the  inside  surface  of  floor 
Description  of  the  first  solid  layer 
Description  of  the  second  solid  layer 
Description  of  the  semi-infinite  earth  layer 
Number  of  layers  for  inner  mass 

Thermal  resistance  at  the  outside  surface  of  the  in- 
terior mass 
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27  Thickness,  density,  specific  heat,  and  thermal  con- 
ductivity of  the  internal  mass 

28  Thermal  resistance  at  the  other  outside  surface  of  in- 
ternal mass 

29  Description  of  the  outside  surface  of  interior  mass 

30  Description  of  the  interior  mass 

31  Description  of  the  other  outside  surface  of  interior 

mass 

32  Blank  card  (necessary  to  show  end  of  above  data) 

33  Run  no.  card 

34-37  Outside  air  temperature 

38-41  Inside  air  temperature 

42  Roof  area,  wall  area,  floor  area,  inner  mass  surface 

area,  air  flow  (ventilation  on  air  leakage),  conductance 
of  door,  door  area,  conductance  of  window,  window  area, 
ground  temperature,  average  inside  air  temperature 


Components  possessing  significant  heat  capacity  (such  as  walls,  roof, 
etc.)  are  described  in  cards  2 through  31,  while  components  having 
negligible  heat  capacity  are  described  on  card  42.  Some  various  options 
available  in  the  Response  Factor  program  are  discussed  below. 
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Additional  layers  can  be  readily  handled.  For  example,  if  the 
roof  contains  a second  layer,  then  the  number  of  layers  given  in  card  2 
would  be  increased  to  four.  Also,  a card  giving  the  thermal  and  physical 
properties  of  this  additional  layer  and  a card  giving  a description  of 
this  layer  would  be  inserted  in  proper  sequence  (from  inside  to  outside). 
Additional  layers  in  any  component  would  be  handled  in  a similar  manner. 
If  the  additional  layer  is  an  air  insulating  layer,  then  only  an  average 
value  of  the  thermal  resistance  of  the  air  layer  would  be  specified  on 
the  card  giving  layer  properties.  Another  option  is  a floor  which  has 
no  semi-infinite  layer  (such  as  the  floor  of  a room  of  a multi-story 
building).  This  case  may  be  handled  by  omitting  the  semi-infinite  layer 
index  on  card  16  and  removing  the  card  giving  the  properties  of  the  earth 
and  the  card  giving  the  layer  description  of  the  earth.  Anotlier  option 
is  the  case  of  a building  without  a component  (such  as  a room  without 
windows).  This  case  is  handled  by  setting  the  area  of  that  component 
(given  on  card  42)  equal  to  zero.  And  finally,  the  option  for  determina- 
tion of  the  inside  room  air  temperature  (floating  test)  is  handled  by 
inserting  four  blank  cards  for  the  inside  air  temperature  (cards  38 
through  41). 
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A print-out  of  the  computer  results  follows  the  sample  set  of 
input  data.  The  first  page  of  the  print-out  of  results  gives  the  de- 
scription and  composition  of  each  building  component  having  significant 
heat  capacity.  The  second  page  gives  the  run  number,  outside  air  tem- 
perature cycle,  the  inside  air  temperature  cycle  (if  thermostated) , and 
the  data  input  given  on  card  42.  And  finally,  the  third  page  of  the 
print-out  of  the  computer  results  gives  the  inside  and  outside  air  tem- 
peratures, the  heat  fluxes  from  the  room  at  the  inside  surfaces  of  the 

-1  -2 

building  components  in  Btu  hr  ft  , the  air  infiltration  loss,  and 
the  net  heat  loss  from  the  room  at  prescribed  time  intervals,  in  Btu  hr 
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Sample  Set  of  Data  Input 
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